Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA39.7% 91.9 / 36.8 73intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (40/4);  new base A (0/29);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 2.15e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTG  >  NC_000913/5010‑5303
                                                                                                                                                       |                                                                                                                                                 
atcatcAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGgag‑a                                                                                                                                                 <  1:945554/151‑1 (MQ=255)
  catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑aa                                                                                                                                                <  2:488408/150‑1 (MQ=255)
  catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AAt                                                                                                                                               >  2:582413/1‑151 (MQ=255)
  catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AAt                                                                                                                                               >  2:91786/1‑151 (MQ=255)
     caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATaac                                                                                                                                            >  2:665383/1‑151 (MQ=255)
     caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATaac                                                                                                                                            >  1:381416/1‑151 (MQ=255)
                aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAATTTCTCAATAAATG‑CGGTCACTTAGAGATTAGGATTGCGGAG‑AATAACAAacgccg                                                                                                                                    >  2:284228/1‑148 (MQ=255)
                 ttCATTATCTAAATCAGGCGGGGTTGGGTTTTAGGCAGCCCGGTTTTTTTATGAAGAAATGATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTAAATAAATG‑CGGTAATTTAGAGATGAGGATTGCGGAG‑AATAACAAccgccg                                                                                                                                    <  2:100570/147‑1 (MQ=255)
                     ttATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATc                                                                                                                            >  2:898890/1‑151 (MQ=255)
                              tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                   >  2:351667/1‑151 (MQ=255)
                              tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                   >  1:656529/1‑151 (MQ=255)
                               cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                  >  2:929176/1‑151 (MQ=255)
                                 ggCGGGTTTTGTTTTTATGCGCCCCGGTTTTTTTATGAAAAAATTAGGGAA‑AAAAATGACGGGAAAAAAGGAAAATTTCCAAATAAATG‑CGGAAATTAAGAGATTAGGATTGCGAAA‑AAAAAAAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                  <  1:619872/149‑1 (MQ=255)
                                 ggCGGGGTTTGTTTTTATGCACCCCGGTTTTTTTATGAAAAAATTAGGAAA‑AAAAATGACAGGAAAAAAGGAAAATTTCCAAAAAAATG‑GGGTAACTTAGAGATTAGGATGGCGGAA‑AAAAAAAACCGCCGTTCTCATCGAGTAAtct                                                                                                                   <  1:224810/148‑1 (MQ=255)
                                  gCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                               >  1:228283/1‑151 (MQ=255)
                                        tttttttttACGCACCCCGGTTTTTTTATAAAAAAATTATGAAA‑AAAAATAAAAGGAAAAAAGGAAAAATTCTCAATAAATG‑CGTTACCTAAGAGATAAGGATTGCGGAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                          <  2:904873/145‑1 (MQ=255)
                                        ttgttttttCTGCAGCCCGGTTTTTTAAAAAAAAATTTAGGAAA‑AAAAATAACAGGAAAAAAGGAGAAATTCTAAAAAAAGG‑CGGTAACTGAGAGATTAGTATTGCGGAA‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                          <  2:437378/145‑1 (MQ=255)
                                            ttttaaagCAACCCGGTTTTTTAATGAAGAAATTAAGAAG‑AAAAAGGGAGGGAAAAAAGGAAAAATT‑TCAAAAAAGGTTGGTATCATAGAGTTTGGGGTTG‑GGAGAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGa                                                                                                        <  2:687465/141‑1 (MQ=255)
                                              tttATGCACCCCGGTTTTTTTATAAAAAAATAAGGAAA‑AAAAAGGACGGGAAAAAAGAAAAAATTCCCAAAAAATG‑GGGAAATTGAGAAATTAGGATGGGGGAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                    <  2:753980/150‑1 (MQ=255)
                                               ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                    >  2:677944/1‑149 (MQ=255)
                                                       ccgggtttttttATGAAAAAATAAGGAAA‑AAAAAGGACGGGAAAAAAGGAAAAATTCTAAATAAAGG‑CGGTAACTTAGAGATTAGGATTGCGGAA‑AAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                           <  2:208873/144‑1 (MQ=255)
                                                        ccGGTTTTTTAAAAAAAAAATTAGGGAA‑AAAAATGAAGGGAAAAAAGGAAAAATTCTAAAAAAAGG‑CGGTAATTAAGAGATTAGGATGGCGGAA‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                           <  1:1012761/149‑1 (MQ=255)
                                                        ccGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                         >  2:783360/1‑151 (MQ=255)
                                                         cGGCTTTTTAATAAAAAAATTATGGAA‑AAAAATGACAGGGAAAAAGGAAAAATTCAAAATAACCG‑CGATAAATTAAATCTTAGGATGGGGCAG‑AATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                        <  1:438342/151‑1 (MQ=255)
                                                          ggtttttttAGAAAAAAATAAGGGAA‑AAAAATGACGGGAAAAAGGGAAAAATTCTAAATAAATG‑GGGTAAATTAGAGATTAGGATTGCGGAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                        <  2:461037/147‑1 (MQ=255)
                                                             ttttttATGAAAAAATTATGAAA‑AAAAAGGACAGGAAAAAAGAAAAAATCCTAAAAAAAGG‑CGTTAATTTAGAAATTAGAATTGCGAAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                      <  2:868436/149‑1 (MQ=255)
                                                                tttATGCAGAAATTTTGGAG‑AAAAATGAAAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCACCGATCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                        >  2:540103/1‑144 (MQ=255)
                                                                  tATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                >  2:974806/1‑150 (MQ=255)
                                                                   agaaaaaaaTTAGGGAG‑AAAAATGAAAGGAAAAAAGAAAAAATTACAAAAAAAGG‑CGGTAATTAAGAGATTAGGATTGGGGAA‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                <  2:303398/143‑1 (MQ=255)
                                                                     gaagaaATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                             >  1:137985/1‑150 (MQ=255)
                                                                      aagaaATTATGAAA‑AAAAAGGACAGGAAAAAAGAAGAATTTCAAAATAAATG‑CGTTACCTTAGAGATTAGGATTGCGGAA‑AATACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                           <  1:1088627/151‑1 (MQ=255)
                                                                       aaaaaTTAGGAAA‑AAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGG‑CGGTAACTAAGAGATTAGGATGGGGGAA‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc                                                                          <  1:506640/149‑1 (MQ=255)
                                                                       aaaaaTTAAGGAAAAAAAATGACAGGAAAAAAGGAAAAATTCAAAATAAAGG‑GGATAAATAAAAGATAAGGATTGCGGAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                             <  1:521327/147‑1 (MQ=255)
                                                                        aaaaTAATGAAA‑AAAAATGACAGGAAAAAGGGAAAAATTCTAAATAAAGG‑CTGTAATTTAGAGATTAGGATTGAGAAA‑AAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                         <  2:350181/150‑1 (MQ=255)
                                                                         aaaTTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                        >  2:238847/1‑151 (MQ=255)
                                                                          aaTTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                       >  1:728806/1‑151 (MQ=255)
                                                                            ttATGAAA‑AAAAATGAAAGGAAAAAAGGAAAAATTCTAAAAAAAGG‑CGTTAACTTAGAGATTAGGATTGGGGAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                      <  1:91806/150‑1 (MQ=255)
                                                                                 gag‑aAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                 >  1:1081365/1‑150 (MQ=255)
                                                                                     aaaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                 >  2:306989/1‑147 (MQ=255)
                                                                                     aaaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                 >  1:717520/1‑147 (MQ=255)
                                                                                      aaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                            >  2:987584/1‑151 (MQ=255)
                                                                                       aaaTGAAAGGAAAAAAGGAGAATTTCTAAATAAATG‑CGTTACCTAAGAGATTAGGTTGGCGGAA‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                           <  1:314271/151‑1 (MQ=255)
                                                                                         aTGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                          >  2:330859/1‑150 (MQ=255)
                                                                                            aCAGGAAAAAAGGAGAAATTCCAAAAAAAGG‑GGGAAATTTAGAGATTGGGATTGCGGAA‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                        <  1:255562/149‑1 (MQ=255)
                                                                                            aCAGGAAAAAAGGAAAATTTCTAAAAAAATG‑CGTTAACTAAGAGTTTGGGATGGCGAAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCGGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                       <  1:438951/150‑1 (MQ=255)
                                                                                            aCAGGAAAAAAGGAAAAATTCTAAATAAAGG‑CGGTAACTTAGAGATTAGGATTGCGGAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                        <  1:869757/149‑1 (MQ=255)
                                                                                               gggAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                   >  1:686955/1‑151 (MQ=255)
                                                                                               gggAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGCGAAAAAGATGCAATCTATc                                                   >  1:392231/1‑151 (MQ=255)
                                                                                               ggaaaaaaGAAAAATTTCTAAATAAATG‑GGTTATCTTAGAGGTAAGGTTGGCGAAA‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                       <  1:629423/144‑1 (MQ=255)
                                                                                                 aaaaaaGGAGAAATTCACAAAAAAGG‑CGTTAACTTAGAGATTAGGATGGCGGAA‑AAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                 <  2:409977/150‑1 (MQ=255)
                                                                                                  aaaaaGGAGAATTTTAAAAAAAGTG‑CGGTAATTAAGAGATTAGGATGGGGGAA‑AAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                  <  2:895575/149‑1 (MQ=255)
                                                                                                  aaaaaGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACTTGTGAAAAAGATGCAATCTATCGt                                                 >  2:738693/1‑150 (MQ=255)
                                                                                                   aaaaGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                               >  1:797598/1‑151 (MQ=255)
                                                                                                      agaaaaaaTTCTCAAAAAAGG‑GGGTAATTAAGAGATTAGGATTTGGGAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                             <  2:228220/145‑1 (MQ=255)
                                                                                                         agaAATTCTAAATAAAAG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACAGCAGTTCTCATCGAGTAATCTCCAGATATCGACCCATAACGAGCAGTGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                         >  2:158629/1‑151 (MQ=255)
                                                                                                          gaAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                        >  2:1065446/1‑151 (MQ=255)
                                                                                                           aaaTTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                       >  1:525943/1‑151 (MQ=255)
                                                                                                                 tcAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                 >  1:824748/1‑151 (MQ=255)
                                                                                                                  cAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                >  1:484902/1‑151 (MQ=255)
                                                                                                                       aagg‑CGGTAACTTAGAAATTGGGATTGCGAAA‑AATACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                           <  2:121400/148‑1 (MQ=255)
                                                                                                                            cGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                       >  2:110144/1‑151 (MQ=255)
                                                                                                                                aaCTTAGAAATTAGAATGGCGAAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                   <  1:331949/151‑1 (MQ=255)
                                                                                                                                aaCTAAGAGATTAGGATGGGGGAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCt                     <  1:306965/149‑1 (MQ=255)
                                                                                                                                  cTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                  >  2:421317/1‑150 (MQ=255)
                                                                                                                                    tAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg               >  2:200054/1‑151 (MQ=255)
                                                                                                                                    tAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATATCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTACCCTGTGAAAAAGATCCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg               >  2:27127/1‑151 (MQ=255)
                                                                                                                                       agaTTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc             >  1:51032/1‑150 (MQ=255)
                                                                                                                                       agaTTAGGATGGCGAAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc             <  1:823152/150‑1 (MQ=255)
                                                                                                                                        gaTTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc             >  1:648323/1‑149 (MQ=255)
                                                                                                                                          ttAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc          >  2:286034/1‑150 (MQ=255)
                                                                                                                                            aGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac        >  1:762537/1‑150 (MQ=255)
                                                                                                                                              gATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg     >  1:6272/1‑151 (MQ=255)
                                                                                                                                                  ggggAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGccg  <  1:899182/148‑3 (MQ=255)
                                                                                                                                                  ggggAA‑AAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGc    <  1:138085/146‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
ATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATG‑CGGTAACTTAGAGATTAGGATTGCGGAG‑AATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTG  >  NC_000913/5010‑5303

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