Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT45.9% 78.3 / 33.4 61N156K (AAC→AAA) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (24/9);  new base T (0/28);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 7.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCC  >  NC_000913/4111597‑4111892
                                                                                                                                                        |                                                                                                                                                  
caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACg‑tttg                                                                                                                                                    >  2:432625/1‑151 (MQ=255)
 aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACg‑tttgt                                                                                                                                                   <  2:508271/151‑1 (MQ=255)
        ttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑c                                                                                                                                             <  2:491265/149‑1 (MQ=255)
         tgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑cccc                                                                                                                                          >  1:735993/1‑151 (MQ=255)
          gctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCccg                                                                                                                                         <  1:543447/151‑1 (MQ=255)
             ttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCccgcc                                                                                                                                       <  1:768565/150‑1 (MQ=255)
              tttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGt                                                                                                                                     >  1:569027/1‑151 (MQ=255)
                tGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGtt                                                                                                                                    >  1:496674/1‑150 (MQ=255)
                   cATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTtagt                                                                                                                                 >  1:425982/1‑150 (MQ=255)
                    aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTtagt                                                                                                                                 >  2:580297/1‑149 (MQ=255)
                     tttATTCAACCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGTTTCAACAGCGAAGGTCACCCCCTTTT‑CGCCTTTGTCCGCGACAACATACCGCAGTGCACCTTCTGTCTGGGTGTAATAACT‑TTTGTTTTC‑CCCCGCCGTtagtag                                                                                                                               <  1:120946/150‑1 (MQ=255)
                         ttCAACCCTTCAAGCGATTCAGATAGCCCCAGCTAAATCGGTTCAACAGCGAGGGTCCCCCCTTTTC‑CGCCGTTGTCCGCAAAAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACG‑TTTTTTTTC‑CCCCGCCGTTAGTAGCg                                                                                                                             <  1:769248/148‑1 (MQ=255)
                            aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGc                                                                                                                        >  1:728622/1‑150 (MQ=255)
                            aaaCCTACAAGCGATTCAGAGAGCGCCAGTTTAATCGGTCCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGAAAACAAAACGCAGTCCCCCTTCTGTCTGGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCa                                                                                                                       <  2:230553/151‑1 (MQ=255)
                             aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGc                                                                                                                        >  1:334410/1‑149 (MQ=255)
                                ctccAAGCGATCCAGATAGCGCCAGCTTAATCGGTTCAACGGCGAAGGTCACCCCTTTTT‑CGCCGTTGCCCGCGACAACATAACGCGGTGCCCTTTCTGTCTCGGTTTAATATCG‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                                    <  2:82210/147‑1 (MQ=255)
                                     aGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                                                                                              >  2:251724/1‑151 (MQ=255)
                                        atttCAGATGGCACCAGTTAAATCTCTGCAACAGCGAAGGCCGCCCCCTTTC‑CGCTGTTGTCCGGGAAAATATAACGCAGTGCGCTTTATGTATCGGTGTAAAATCT‑TTTGTTTTA‑CCCCGCCGTTAGTAGCGACTGCAGttt                                                                                                                   <  1:14856/141‑1 (MQ=255)
                                          ttCAGATAGCGCCAGTTAAAGGGGTTAAACAGGGAGGGTAACCCCCTTTT‑CCCCTTTGTCCGCAAAAAAATACCGCTTTCCCCCTTGTGTCCCGGTGTAATACCT‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                                                                                           <  2:255716/149‑1 (MQ=255)
                                             aGATAGCGCCATTTTAATCGTTAAAACAGCGAAGGCCCCCCCTTTTG‑CCCTTTTCTCCGCCAAAACATACCGCAGTCCTCCTTCTGTCGGGGTTTAATATTT‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                                      <  1:458492/151‑1 (MQ=255)
                                             aGATAGCGCCAGCTTAATCGTTCCAACAGCGAAGGTCACCCCCTTTC‑CGCCGTTGTCCGCGACAAAATAACGCAGTGCCCCTTCTGTCTGGGTGTAATAATG‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                                                                                        <  1:263980/149‑1 (MQ=255)
                                               atggcgcCAGTTAAACCGGTTAAACAGCGAAGGCCACCCCTTTTT‑CGCCGTTTTCCGCAACAAAAAACCGCAGTCCCCCTTCTGTCTCGGTTTAAT‑ACG‑TTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTgtgt                                                                                                    <  2:702709/148‑1 (MQ=255)
                                               aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTgtgt                                                                                                    >  1:680234/1‑151 (MQ=255)
                                                 aGCACCAGCTTAACCGGTACAACAGGGAAGCCCGCCCCCTTTT‑CGCGGTTGCCCGCAACAAAATAGCGCATTGCGCCTTCTGTCGGGGTGAAATCTCT‑TTTTTTTTC‑CTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                                      <  2:595572/147‑1 (MQ=255)
                                                   cgcCAGTTTAATCGGTTAAACAGCGAAGGCCACCCCCTTTC‑CGCCGTTGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCGGGGTGTAATAACT‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                                     <  2:590512/146‑1 (MQ=255)
                                                   cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                                >  2:322116/1‑151 (MQ=255)
                                                               cgttaaaaCAGCGAAGGACACCCCCTTTT‑CGCCGTGGCCCGCAACAACATAACGCAGTGCCCCTTTTGTCTGGGTTTAATAATG‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATg                                                                                      <  1:342508/143‑1 (MQ=255)
                                                               cGGTACAACAGCG‑AGGCCCCCCCCTTTTGCGCCTTTGTCCGCGAAAACATAACGCAGTCCCCCTTCTGTCGCGGTTTAATAGTT‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg                                                                                     <  1:732325/150‑1 (MQ=255)
                                                                ggTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg                                                                                   >  1:595059/1‑151 (MQ=255)
                                                                ggTCCAACAGGGAAGGCCACCCCCTTTC‑CGCCGTTGCCCGCAAAAACATAACGCAGTCCCCCTTCTGTCTCGGTTTAATAATG‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg                                                                                     <  2:343349/149‑1 (MQ=255)
                                                                     aaaaGCGAAGGTCACCCCCTTTT‑CTCTTTTTTCCGCGACAACAAAACGCAGTGCCCTTTCTTTCTGGGTTTATTACTG‑TTTTTTTTT‑CTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt                                                                              <  1:319425/148‑1 (MQ=255)
                                                                     aaCAGCGAAGGTCCCCCCCTTTC‑CGCCTTTGCCCGCAACAACATACCGCAGTCCCCCTTCTGTCTCGGTTTAATATCT‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgg                                                                              <  2:153013/151‑2 (MQ=255)
                                                                      aCAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtg                                                                             >  1:114976/1‑151 (MQ=255)
                                                                             aggccccccccTTTTCCGCGTTTGCCGGCGAAAACAAAGCGCGTTGCCCCTTTTGTCT‑GGTTAAAAAACGTTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg                                                                         <  2:547065/142‑1 (MQ=255)
                                                                                 cacccccTTTC‑CTCCGTTGTCCCGAACAACAAACCGCGGTGCCCTTTTTGTTTCGTTTTAATAGCT‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCa                                                                    <  1:198232/146‑1 (MQ=255)
                                                                                    ccctTTTG‑CGCGTTTGCCCGCGAAAACATACCGCGGTGCCCCTTTTGTCTCGTTGTAATAGTG‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                               <  1:113902/151‑1 (MQ=255)
                                                                                    ccccTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                               >  2:313794/1‑151 (MQ=255)
                                                                                    ccccTTTC‑CGCCGTTGTCCGCAACAACATAACGCGTTGCCCCTTCTTTCGGGGTGTAATAACG‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                               <  2:736067/151‑1 (MQ=255)
                                                                                             cgccttggcccGCAAAAACATAACGCGGTCCCCTTTCTGTCTGGGTGTAATATTT‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTa                                                       <  1:314493/142‑1 (MQ=255)
                                                                                                  ttGCCCGCAAAAACAAACCGCATTCCACTTTTTGCCTGTGTTTAATACCT‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGt                                                        <  2:194392/145‑1 (MQ=255)
                                                                                                   tGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTg                                                 >  1:766593/1‑151 (MQ=255)
                                                                                                    gTCCGCAAAAACATAACGCGGTGCACCTTCTGTCTCGGTTTAATATCT‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTg                                                 <  1:712936/150‑1 (MQ=255)
                                                                                                     tCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTg                                                 >  1:762023/1‑149 (MQ=255)
                                                                                                      ccGCGACAACATACCGCAGTCCCCCTTCTGTCTGGGTGTAATAGCT‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa                                                <  2:342440/149‑1 (MQ=255)
                                                                                                           acaacaTACCGCTGTCCACCTTCTGTCTCGGTGTAATATTT‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttc                                          <  2:139038/150‑1 (MQ=255)
                                                                                                           acaacaTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcg                                         >  1:455023/1‑151 (MQ=255)
                                                                                                               caTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                                     >  1:761433/1‑151 (MQ=255)
                                                                                                                 tAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGctttcgctttcg                                   >  2:348269/1‑151 (MQ=255)
                                                                                                                    cGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCgg                                >  1:511605/1‑151 (MQ=255)
                                                                                                                     gagttgccctttCTGTCTGGTTGTATTTCGG‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                                 <  1:498200/139‑1 (MQ=255)
                                                                                                                     gCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                               >  2:169552/1‑151 (MQ=255)
                                                                                                                         tGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCtttt                           >  2:250125/1‑151 (MQ=255)
                                                                                                                           cACCTTCTGTCTCGTTGTAATACCT‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTg                         <  1:63139/151‑1 (MQ=255)
                                                                                                                            cctttCTGTCTGTGTTTAATATTG‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttttt                          <  1:395699/146‑1 (MQ=255)
                                                                                                                              cTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc                      >  2:479532/1‑151 (MQ=255)
                                                                                                                                tctgtctgGGTGTAATATCT‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct                     <  2:558546/150‑1 (MQ=255)
                                                                                                                                     ctcGGTATAAAATTG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc                 <  2:27968/149‑1 (MQ=255)
                                                                                                                                       cGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg             >  1:156766/1‑151 (MQ=255)
                                                                                                                                            taataaCG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg        >  1:122987/1‑151 (MQ=255)
                                                                                                                                             aatatCG‑GTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTGGCTCTGGCCCCTGGAt          <  1:489718/148‑1 (MQ=255)
                                                                                                                                             aataaCT‑TTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGAt          <  2:728696/148‑1 (MQ=255)
                                                                                                                                              atatCG‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGAt          <  2:687906/147‑1 (MQ=255)
                                                                                                                                                aaCT‑TTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa           <  2:115008/144‑1 (MQ=255)
                                                                                                                                                     ttttttttT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGcc  <  2:458864/149‑1 (MQ=255)
                                                                                                                                                     ttttttttC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTa     <  1:545774/146‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
CAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACG‑TTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCC  >  NC_000913/4111597‑4111892

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