Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913827,5070TG47.9% 143.4 / 27.3 96D492A (GAC→GCC) ybhFputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (12/38);  new base G (45/1);  total (57/39)
Fisher's exact test for biased strand distribution p-value = 5.82e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGG  >  NC_000913/827363‑827656
                                                                                                                                                |                                                                                                                                                         
cTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGc                                                                                                                                                   <  2:1593054/151‑1 (MQ=255)
      cttAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGCAg                                                                                                                                            >  2:1548353/2‑151 (MQ=255)
      aTTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAg                                                                                                                                            >  1:134627/1‑151 (MQ=255)
       ttAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGGGACGCCGGGAgg                                                                                                                                           >  2:1562367/1‑150 (MQ=255)
             tCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtc                                                                                                                                     <  2:497125/151‑1 (MQ=255)
                 gcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAAACCCCGGAGGTGGGtcgttg                                                                                                                                  >  1:165445/1‑146 (MQ=255)
                 gcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGTttgtcg                                                                                                                                 >  1:687237/1‑151 (MQ=255)
                      acacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGTGGGGGGCAACGCGGGCAGTAGGCtcctgcagaga                                                                                                                            >  2:44346/1‑142 (MQ=255)
                      acacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGaaa                                                                                                                            >  1:1467415/1‑151 (MQ=255)
                       cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCACAATTCACGGCGGGTGGGGGGGTCACCGCCGGGAGTCGGTTCGTCGAGaaaa                                                                                                                           >  2:1105786/1‑150 (MQ=255)
                       cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAAc                                                                                                                           <  1:627672/151‑1 (MQ=255)
                       cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAAc                                                                                                                           <  1:586720/151‑1 (MQ=255)
                       cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCCACGCCGGAAGGCGGTTCGCCCAGaaaa                                                                                                                           >  2:1366503/1‑150 (MQ=255)
                         ccAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAg                                                                                                                         <  1:504394/151‑1 (MQ=255)
                           aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAGTCGGTTCGTCCGGAAACAGaa                                                                                                                       >  2:280284/1‑151 (MQ=255)
                           aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGTCAACCCCGGGAGGCGGGTCGTCGGGAAACaaaa                                                                                                                       >  1:92373/1‑151 (MQ=255)
                           aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAACCCCGGAGGTGGGttgtccgggaaacgaa                                                                                                                        >  2:369215/1‑137 (MQ=255)
                            ggCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGGCGAGAAAAAGaaa                                                                                                                      >  1:1046624/1‑150 (MQ=255)
                            ggCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCCTCGAGGAACAGAAt                                                                                                                      >  1:1001095/1‑151 (MQ=255)
                            ggCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGACGGGTCGTGGGGAAAAAGaaa                                                                                                                      >  1:366934/1‑150 (MQ=255)
                               cGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                   <  2:1307894/151‑1 (MQ=255)
                                 aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                                                                                                  <  2:1577605/150‑1 (MQ=255)
                                 aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                                                                                                  <  1:1429961/150‑1 (MQ=255)
                                      gTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGTCGGGTCGTCGGGAAAAAAAAAGTCCGGGTCa                                                                                                            >  1:885103/1‑151 (MQ=255)
                                          cAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGTCAAAGCCCGAAGCTGCGCTGGCGGGAAACACAACGTCCGGGTCATGc                                                                                                         >  2:698585/1‑150 (MQ=255)
                                          cAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCa                                                                                                        <  1:995220/151‑1 (MQ=255)
                                          cAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCa                                                                                                        >  1:170230/1‑151 (MQ=255)
                                              atTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                                                                                                    <  2:178398/151‑1 (MQ=255)
                                              atTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAATCGGGTCGTCGGGAAACAGAAAGTCCCGTTCCTGCATTAg                                                                                                    >  1:992386/1‑151 (MQ=255)
                                                        aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcc                                                                                          <  2:1201729/151‑1 (MQ=255)
                                                         tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCGGAAGGCGGGTTGTCGAGAAACAGAAAGTCCGGGTCATGCCTCAGCGAAAAGgcc                                                                                          >  1:1338456/1‑150 (MQ=255)
                                                         tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAAGGGGGTCGGCGAGAAAAAAAAAGGCCGGGTCATGCATAAGCGAACAgggc                                                                                          >  1:473276/1‑148 (MQ=255)
                                                         tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGAAACAGAAAGTCCGGGTTATGCATCAGCGAACAGgcc                                                                                          >  1:250874/1‑150 (MQ=255)
                                                         tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGGCGGGTTGTCGAGAAAAAGAAAGTCCGGGTCATGGATCAGCGAACaaggca                                                                                         >  1:841647/1‑146 (MQ=255)
                                                         tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGTAACGCCGGGAGGCGGGTCGGCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACACgcc                                                                                          >  2:703422/1‑150 (MQ=255)
                                                         tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGTTTGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAAAAAgcca                                                                                         >  1:553186/1‑151 (MQ=255)
                                                         tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGTTTGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACaaggc                                                                                          >  2:1520103/1‑146 (MQ=255)
                                                         tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGGAGTCGGGTCGTCGAGAAAAAGAATGTCCGGGTTATGCATCAGCGCACAAgcc                                                                                          >  2:602964/1‑150 (MQ=255)
                                                             tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgc                                                                                     <  2:169955/151‑1 (MQ=255)
                                                             tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAAAAAgccagcg                                                                                      >  2:191915/1‑150 (MQ=255)
                                                               aaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACCCCGGAAGGCGGTTCGGCGGGAAACAAAAAGGCCGGGTCATGCCTCAGCGAACAGgcccgcgcc                                                                                    >  1:1000949/1‑150 (MQ=255)
                                                                     gtggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGAAGTCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATGCATCAGAGAAAAGGCCAGCGCCCGgcgc                                                                             >  1:940780/1‑150 (MQ=255)
                                                                         tgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGttt                                                                        <  1:905632/151‑1 (MQ=255)
                                                                              aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAc                                                                   <  1:872727/151‑1 (MQ=255)
                                                                              aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAc                                                                   <  1:1379070/151‑1 (MQ=255)
                                                                              aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAc                                                                   <  1:1213544/151‑1 (MQ=255)
                                                                              aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAc                                                                   <  2:948001/151‑1 (MQ=255)
                                                                              aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGGACAGGCCCGCGCCCGCCGCCGTTTaaaaa                                                                   >  1:209139/1‑150 (MQ=255)
                                                                               tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAGGCGGGTTGTCGAGAAAAAAAATGTGCGGGTTTTGGCTCAGGGGACAAGCCAGCGCCCGCCGCTGTTTTAAAca                                                                  >  2:1271589/1‑150 (MQ=255)
                                                                               tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGTTTGTCGGGAAACAGAAAGTCCGGGTCCTGCATCACCGAACAGGCCCGCGCCAGCCCCCGGTTTaaaa                                                                    >  1:963537/1‑150 (MQ=255)
                                                                                  ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCCGGAGGCGGGTCGTCGAGAAAAAGAATGTCCGGTTTCTGCATCAGCGAACAGGCCCGCGCCAGCCGCCGTTTTAAAACTaa                                                               >  2:709607/1‑151 (MQ=255)
                                                                                   cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCCAAGCCGGAAGTCGGGTTGTCGGGAAAAAGAATGGCCGGGTCATGCCTCAGCGGACAGGCCCGCGCCACCCGCTGTTTTAAAACTAAt                                                              >  1:1366567/1‑151 (MQ=255)
                                                                                       aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                           <  2:543336/150‑1 (MQ=255)
                                                                                       aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                           <  2:603724/150‑1 (MQ=255)
                                                                                       aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                           <  2:1083486/150‑1 (MQ=255)
                                                                                       aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                           <  2:1292530/150‑1 (MQ=255)
                                                                                        cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCCTCCAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAGGCCACCGCCAGCCGCCGGTTTAAACCTAATAGCAg                                                         >  2:738438/1‑151 (MQ=255)
                                                                                         gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGTAAAGCCGGAAGTCGGTTCGTCGAGAAAAAAAAAGTCCGGGTCATTCATCAGAGAACAGGCCAGCGCCAGCCGCCGTTTAAAAAATAATGGCAgg                                                        >  1:1255085/1‑150 (MQ=255)
                                                                                         gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACCCCGGAAATCGGTTCGTCGAGAAAAAAAAAGTCCGGGTTCTGCATCAGCGAACAGGCCAGCGCCCCCCGCCGTTTTAAAACCAATTGCCGt                                                        >  1:491852/1‑151 (MQ=255)
                                                                                          ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGTGTGAGGGGGTCAACGCCGGAAGTCTGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCGGCGCCAGCCGCTGTTTAAAGCCTAATGGCAGt                                                        <  2:1057153/150‑1 (MQ=255)
                                                                                            tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGTTGGGGGGGGTAACGCCCGGAGGCAGGTCCTCCAGAAAACGAATGTTCGGGTTCTGCATCAGCGCACAAGCCCGCGCCACCCCCCGTTTCTAACCTAATCGccagtc                                                      >  2:60487/1‑145 (MQ=255)
                                                                                            tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCATCGAACAGGCCAGCGCCCGCCGCCGTTTTAAAACCAATGGGAggt                                                       >  2:850293/1‑147 (MQ=255)
                                                                                              tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATc                                                   <  1:186282/151‑1 (MQ=255)
                                                                                               ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCg                                                  <  1:33228/151‑1 (MQ=255)
                                                                                               ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGTTCCTCGAGAAACAGAAAGTCCGGTTCATGCCTCAGCGAACAGGCCAGCGCCAGCCGCTGTTTTAAACCTTATGGGCGATCATCg                                                  >  2:1308163/1‑151 (MQ=255)
                                                                                                tctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCgg                                                 <  1:892933/151‑1 (MQ=255)
                                                                                                   aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCTACTCCGCGGAGCGGTTGCACGAGAAAAACGAAGCCCCGGTCCTGCCTCCGACCACAGCGCACCGCCAGCCGCTGTTTTTAACCTAATTGCAGTTcgacac                                                 >  2:821154/1‑143 (MQ=255)
                                                                                                         cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGGTTGTCCCGAAAAAGAAAGTCCGGGTCATGCCTCAGCGGAAAGGCCAGCGCCAACCCCCGTTTTAAAACTAAATGCGGTTTATCGgggggggggg                                        >  1:756227/1‑151 (MQ=255)
                                                                                                         cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGAAGGCGGGTCGTCGAGAAACAGAAAGTCCGGGTCCTGCATCAACGAACAGGCCCGCGCCAGCCGCCGGTTTAAACCCAATGGCAGTGTGTCGGTgggggggg                                        >  2:524766/1‑151 (MQ=255)
                                                                                                         cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGGAGACGGGTCGTCGAGAAACAGAAAGTCCGGTGCATGCATCAGCGAACAGGCCAGCGCCCACCGCCGTTTTAAAACTAAAGGCAGtttttctggtgcggggg                                        >  2:246418/1‑136 (MQ=255)
                                                                                                          tGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGa                                       <  1:369319/151‑1 (MQ=255)
                                                                                                           gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACCCCGGGAATCGGTTTTTCGTGAAACAAAATGTCCGGTTCCTGTCACAGCGAACAGGCCAGCGCCAGCCCCTGTTTTAAACCTAATGGGCGGTCCTCGgtggggtggggg                                      >  1:81504/1‑149 (MQ=255)
                                                                                                           gTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGTTTGTCGAGAAAAAGAAAGGCCGGGTTATGCATCCGCGAACAGGCCCGCGCCAGCCGCCGTTTAAAAACCAATGGCCGGTCATCGGTGGGGGGGGAg                                      >  2:832755/1‑151 (MQ=255)
                                                                                                           gTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAAAAGGCCCGCGCCCGCCGCCGTTTAAAAACTAAAGGCAGTTTATCGgggggggggggg                                      >  1:992725/1‑149 (MQ=255)
                                                                                                            ttGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAgg                                     <  2:1036565/151‑1 (MQ=255)
                                                                                                                  tgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATAc                               <  2:1000453/151‑1 (MQ=255)
                                                                                                                    cAGCCAAAACCCGCGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTGCGTCGATGAACAGAATGTCCGGTGCATGCATCAGCGAACAGGCCAGCGCCAGCCGCGGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTTGGAGGTGATACt                              <  2:751871/150‑1 (MQ=255)
                                                                                                                       ccAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTc                          <  1:353001/151‑1 (MQ=255)
                                                                                                                         aaaaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAg                        <  1:820674/151‑1 (MQ=255)
                                                                                                                           aaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGc                      <  2:288905/151‑1 (MQ=255)
                                                                                                                           aaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGc                      <  2:538665/151‑1 (MQ=255)
                                                                                                                           aaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGc                      <  2:552967/151‑1 (MQ=255)
                                                                                                                           aaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGc                      <  2:1403977/151‑1 (MQ=255)
                                                                                                                           aaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGc                      <  1:661248/151‑1 (MQ=255)
                                                                                                                                aCGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAAc                 <  2:96950/151‑1 (MQ=255)
                                                                                                                                 cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACg                <  1:833215/151‑1 (MQ=255)
                                                                                                                                   gcggGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGGAGTTCCTCGGGGGCGGGGGGGGGGATACTTTTCAGGCCGAAAGc               >  2:78929/1‑150 (MQ=255)
                                                                                                                                     gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCt             >  2:709865/1‑150 (MQ=255)
                                                                                                                                     gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCt             >  1:615726/1‑150 (MQ=255)
                                                                                                                                         gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCt         >  2:928546/1‑150 (MQ=255)
                                                                                                                                           gggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCAt      >  1:1220852/1‑151 (MQ=255)
                                                                                                                                            ggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCa       <  2:885429/149‑1 (MQ=255)
                                                                                                                                            ggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATg     >  1:982499/1‑151 (MQ=255)
                                                                                                                                              ggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCg   >  2:1568349/1‑151 (MQ=255)
                                                                                                                                              ggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCg   >  2:204876/1‑151 (MQ=255)
                                                                                                                                               gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCgg  >  1:1162482/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                         
CTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGG  >  NC_000913/827363‑827656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: