Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913345,3280CA40.5% 121.8 / 53.9 84intergenic (+337/‑76)yahL/yahMuncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (45/5);  new base A (0/34);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.56e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345180‑345477
                                                                                                                                                       |                                                                                                                                                     
aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCaa                                                                                                                                                     <  2:382236/151‑1 (MQ=255)
 aaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCaa                                                                                                                                                     >  2:112684/1‑150 (MQ=255)
   aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCaaaaa                                                                                                                                                  <  2:513817/151‑1 (MQ=255)
    aTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAAc                                                                                                                                                 <  2:1229386/151‑1 (MQ=255)
     tGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAAcc                                                                                                                                                >  2:294246/1‑151 (MQ=255)
       aTCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGc                                                                                                                                              >  2:1291268/1‑151 (MQ=255)
              gCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGc                                                                                                                                       >  2:759062/1‑151 (MQ=255)
                aGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCtc                                                                                                                                     >  2:700372/1‑151 (MQ=255)
                    gTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCtcatc                                                                                                                                  >  1:1005330/1‑150 (MQ=255)
                    gTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCtcatc                                                                                                                                  >  2:210411/1‑150 (MQ=255)
                     ttttGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCtcatc                                                                                                                                  >  1:826619/1‑149 (MQ=255)
                      tttGTGATTGTGTGGGTTGACATTAATAAAAAAAAAATCATAATCCAATCAT‑TTTGTGGGGGGTTTTTTCTTTTTTCTTAC‑GCACCACCCCCACGCCTAACAAAAAATATGTGGAAAGC‑CTAAAATAAAAAACCGCCATTAGCTCATCAg                                                                                                                                <  1:782110/150‑1 (MQ=255)
                      tttGTGAATGTGTTGGTTGACATTCATATGAAAAAAATAATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAgg                                                                                                                               >  1:490521/1‑151 (MQ=255)
                          tgAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGa                                                                                                                              >  2:1294979/1‑148 (MQ=255)
                              tgtgTTGTTTGACATTCATAAAAAAAAAATAAAAATCCCATCTTTTTTGGGGGGTTTTTTTTT‑TGTCTTTCAC‑ACATCGCCTCCAACCCTAAAAAAAAATATGCTTAAAGG‑TAAATACAAAAAACCGCCATTAGCTCATCAGGAagag                                                                                                                          <  2:328001/148‑1 (MQ=255)
                              tgtgTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                         >  2:1286567/1‑149 (MQ=255)
                                tgtTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                                                     >  2:1011941/1‑151 (MQ=255)
                                   tggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                                  >  2:339153/1‑151 (MQ=255)
                                            ttCATAAAAAAAAAACCATAAGCCCATAATGTTTGTGGGGGTTCTTGTC‑TGTTTCTGAC‑CCACCCCGCTCACGCCTGACAAAAAATACGCTTTAGTA‑CACAGATAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGA‑AACCAGGa                                                                                                          <  2:960199/149‑1 (MQ=255)
                                                    aaaaaaaTCAGAATCCTACTATGTTTGGGGGGGTTCTTTTC‑GGTATTTAACGGCACCG‑ACTCAAGCCGAACAAAATGTTGGCCGAAAGG‑CTCAGATCAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATtgttgt                                                                                                   <  2:189027/149‑1 (MQ=255)
                                                     aaaaaaACATAATCCCATCTTTTTTGTGGGGGTTTTTTTC‑TGTATCTTAC‑GCACCCCCCTCAACCCGGAAAAAAAATTGGTTGTAAGG‑CTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg                                                                                                <  2:113173/151‑1 (MQ=255)
                                                           tCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgg                                                                                          >  1:1294504/1‑151 (MQ=255)
                                                             aaacgcccATCTTTTTTGGGGGGGTTTTTTTC‑TTTATCTAAC‑CCACCGCCCAAACGCCTAAAAAAAAATATGCTGTAAGG‑CACATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACggg                                                                                         <  1:583193/144‑1 (MQ=255)
                                                             aTAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                                                        >  1:951722/1‑151 (MQ=255)
                                                             aTAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                                                        >  1:580525/1‑151 (MQ=255)
                                                               aaTTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                                                      >  2:317741/1‑151 (MQ=255)
                                                               aaTCCCATCATTTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTAAAGCCTAAAAAAAAATATGCTGAAAGG‑CTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg                                                                                           <  2:250635/146‑1 (MQ=255)
                                                                aTTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                                                      >  1:193252/1‑150 (MQ=255)
                                                                   ccATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                                    >  1:1029739/1‑149 (MQ=255)
                                                                       catGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                                              >  2:351390/1‑151 (MQ=255)
                                                                         tttttGGGTGGGTTTTTTTC‑TTTCTTTGCC‑GCACCGCCCCAAAGCCTAACAAAAAATATGCTGAAAGG‑CTAAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                             <  1:401801/148‑1 (MQ=255)
                                                                              gtgtgGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                                        >  2:701269/1‑150 (MQ=255)
                                                                              gtgtgGGGCCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGc                                                                       >  1:1161255/1‑151 (MQ=255)
                                                                                gtgGGGTCTTTTC‑TTTTTCTTAC‑TCACCGCACCAAACCCTAACAAAAAAGAGCCTGAAGGC‑CTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                                        <  1:208934/148‑1 (MQ=255)
                                                                                gggggtttttttC‑TTTATCTCAC‑CCACCCCACACACGCAAAAAAAAAAATATGCGATGATG‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                                         <  1:331712/139‑1 (MQ=255)
                                                                                  gtttttttttC‑TGTATCTTCC‑GCGTCCCACACAGGCCAGACAGAAAATATGTTGTAAGG‑TTTTTATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt                                                                    <  2:1084977/144‑1 (MQ=255)
                                                                                  gggtttttttC‑TTTATTTTAC‑CCACCCCACAAAACCCTAAAAAAAAATATG‑TTTAAGGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc                                                                   <  2:79762/145‑1 (MQ=255)
                                                                                  gggtttttttC‑TTTATCTTAC‑CCATCGCTCTCAAGCCAGACAAAAAATATGCTGAAAGC‑CCAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt                                                                    <  2:305716/144‑1 (MQ=255)
                                                                                  ggggTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc                                                                   >  1:328660/1‑151 (MQ=255)
                                                                                    ggTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTcc                                                                  >  1:318402/1‑150 (MQ=255)
                                                                                     gTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                                                 >  2:1181889/1‑150 (MQ=255)
                                                                                         tttC‑TGTATCTTAC‑GCATAGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                            >  1:299555/1‑151 (MQ=255)
                                                                                          ttC‑TGTTTCTTCC‑GCTACCCACAAAAGCCAAAAAAAAAATTTGCGGAAAGG‑CTAATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                            <  1:1191565/150‑1 (MQ=255)
                                                                                                     ttccccaccccacacaAGCCGAAAAAAAAATATGTTGAAAGG‑TAAAAAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActc                                                   <  1:92333/140‑1 (MQ=255)
                                                                                                     tAC‑GCATCCCACAAAAGCCTGAAAAAAAATGTTCTTTAAGT‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtttt                                                           <  1:29695/141‑1 (MQ=255)
                                                                                                       c‑GCATCGCCCTCAAGAGAACCAAAAAATATGTTGTAGGC‑GTCATAGAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtcttc                                                <  2:437516/150‑1 (MQ=255)
                                                                                                        cgcaccgcctcaaaGCCTAAAAAAAAATATGCTGAAAGG‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActct                                                  <  2:585511/136‑1 (MQ=255)
                                                                                                         gcgacccACTCAACCAGGAGAAAAATTATGCGGTAAGC‑CTCATATAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAc                                                     <  1:420274/138‑1 (MQ=255)
                                                                                                          cctcccccTCAAGCCTAAAAAAAAATATGCTGAAGGC‑CTAATATAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtcttct                                               <  1:1261077/142‑1 (MQ=255)
                                                                                                           accGCACTCAGGCCGAAAAAAAAATATGTTGAAAGC‑TACATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                                 <  1:294236/144‑1 (MQ=255)
                                                                                                           aTCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                             >  1:266415/1‑150 (MQ=255)
                                                                                                           aTCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                             >  2:222213/1‑150 (MQ=255)
                                                                                                             cccACCCACGCCTCACAAAAAATAGGCTGTAAGT‑CTAAAAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                             <  1:47195/146‑1 (MQ=255)
                                                                                                             cGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          >  2:645678/1‑151 (MQ=255)
                                                                                                             cGCACCCAACCCTAACAAAAAATATGCTGAAAGG‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt                                           <  2:1052245/150‑1 (MQ=255)
                                                                                                             cGCACACAACCCTAAAAAAAAATTTGTTGTAGGA‑CTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          <  1:465208/151‑1 (MQ=255)
                                                                                                                acccAAGCCTGAAAAAAAATATGCTGAAAGG‑CTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                                 <  2:92939/138‑1 (MQ=255)
                                                                                                                 cTCAAGCCTAACAGAAAATAGGCTGTAAGT‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCa                                         <  1:94885/148‑1 (MQ=255)
                                                                                                                  tCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          >  1:91411/1‑146 (MQ=255)
                                                                                                                    aaGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATccc                                   >  2:1255429/1‑151 (MQ=255)
                                                                                                                       ccTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGAACGGGGTTCGAGGCCTCGATGGCGGGCCAAATTTTCAAAACTCTTCTGAATCAGATCCCCtt                                >  1:984394/1‑151 (MQ=255)
                                                                                                                         aaaaaaaaaaTATGCTGAAGGC‑CTAATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttat                              <  1:1270244/145‑1 (MQ=255)
                                                                                                                           aCAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttatt                             >  1:1175095/1‑150 (MQ=255)
                                                                                                                            cggAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                           >  1:772447/3‑151 (MQ=255)
                                                                                                                            atgaattttgtcTGTAGGC‑CCCAAATCAAACCCCGCCATTAGCTCATCAGGAAGAGCAGACGAAAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCAGAATCAGATCCCCt                                 <  1:932718/134‑1 (MQ=255)
                                                                                                                            aaaaaaaTATGCTGAAGGC‑CTCATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttatt                             <  1:455074/146‑1 (MQ=255)
                                                                                                                             aaaaaaTAGGCTGTAAGG‑TTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                           <  2:1283268/148‑1 (MQ=255)
                                                                                                                             aGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAGAACTCTTCTGAATCAGATCCCCTTATTatc                          >  2:715709/1‑151 (MQ=255)
                                                                                                                             aGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatc                          >  2:1263119/1‑151 (MQ=255)
                                                                                                                              aaaaattggatgaaggg‑CACATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatca                         <  1:335734/136‑1 (MQ=255)
                                                                                                                              aaaaaTAGGCTGAAAGT‑CAAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatca                         <  1:781026/150‑1 (MQ=255)
                                                                                                                               aaaaTATGTTGAAAGC‑CTAATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatc                          <  2:290567/149‑1 (MQ=255)
                                                                                                                                 aaTATGTTGTAAGG‑CTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATcc                      <  1:432441/151‑1 (MQ=255)
                                                                                                                                   tatGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt                    >  1:658853/1‑151 (MQ=255)
                                                                                                                                   tatGCTGAAAGG‑CTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATccc                     <  2:241845/150‑1 (MQ=255)
                                                                                                                                     tGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGc                  >  1:587817/1‑151 (MQ=255)
                                                                                                                                        tGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa                >  2:481365/1‑150 (MQ=255)
                                                                                                                                        tGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa                >  2:1143714/1‑150 (MQ=255)
                                                                                                                                          tAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAg              >  1:247092/1‑150 (MQ=255)
                                                                                                                                          gaaagttcaatatAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa                <  2:163729/140‑1 (MQ=255)
                                                                                                                                            aGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg           >  2:197144/1‑151 (MQ=255)
                                                                                                                                             gg‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg           >  2:421099/1‑150 (MQ=255)
                                                                                                                                                   atatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACAtt     >  2:156134/1‑151 (MQ=255)
                                                                                                                                                      tCAAAAACCGCCATGAGCTCATCAGGGAGAGCAGACGACAACCAGGATTGTTGTATGGAACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGACTCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGcc  >  1:1270578/1‑151 (MQ=255)
                                                                                                                                                       cAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGcc  >  2:771322/1‑150 (MQ=255)
                                                                                                                                                       cAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGcc  >  2:1164867/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
AAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTAC‑GCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGG‑CTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345180‑345477

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