Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,621,1180AC44.4% 201.0 / 30.9 125V7G (GTG→GGG) eamAcysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (51/18);  new base C (0/55);  total (51/74)
Fisher's exact test for biased strand distribution p-value = 8.42e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

ACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620971‑1621268
                                                                                                                                                     |                                                                                                                                                         
  tttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CAcccc                                                                                                                                                       >  2:286545/1‑150 (MQ=255)
   ttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CAcccc                                                                                                                                                       <  2:530627/149‑1 (MQ=255)
   ttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCa                                                                                                                                                      <  2:60160/150‑1 (MQ=255)
     cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑Tc                                                                                                                                                   <  1:864189/150‑1 (MQ=255)
               aCAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑c                                                                                                                                        >  1:770328/1‑150 (MQ=255)
               aCAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑c                                                                                                                                        <  1:1563638/150‑1 (MQ=255)
                   aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCt                                                                                                                                   <  1:372929/151‑1 (MQ=255)
                   aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCt                                                                                                                                   <  2:376327/151‑1 (MQ=255)
                    aaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTa                                                                                                                                  <  1:1688934/151‑1 (MQ=255)
                       gATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCg                                                                                                                               <  1:1157320/151‑1 (MQ=255)
                          aGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCct                                                                                                                            <  1:1387406/151‑1 (MQ=255)
                             cgaaaaaGCGACCAGAAAAAACCGCAAGCCGGCCACCTAAAGCGGTGCCTTGTTATAAAGCCCCACTTTG‑ATGACCGCAAAATTTAGCCCCCATTCGACCACTCCCAGTGGCGCCAA‑CCCCCCA‑CCTTTTCGCGA‑CATTCTACCGCCTCTg                                                                                                                         <  1:1117629/148‑1 (MQ=255)
                              ggAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTg                                                                                                                         >  1:301417/1‑150 (MQ=255)
                                   aGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTc                                                                                                                   >  1:22939/1‑151 (MQ=255)
                                   aGCGACCAGAATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAACCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTc                                                                                                                   <  2:1437195/151‑1 (MQ=255)
                                    gCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTc                                                                                                                   <  1:836354/150‑1 (MQ=255)
                                    gCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTc                                                                                                                   >  2:260057/1‑150 (MQ=255)
                                     cGACCAGCATAAAGCGAAAACCGGCCAGCATCAGCGGTGGCATGTTAGGAACCCCCACTTTG‑ATGCCCACAAAATTTAGCCCCCATGCGCCCTCTCCCAGTAGCTCCAC‑CCCCCCA‑TCTTTCCGCGA‑CATCCTACCGCCTCTGAATTTCa                                                                                                                  <  1:628336/150‑1 (MQ=255)
                                        ccAGCATAAAGCGCAAACCGGCAAGCACCAGCGGTGGCATGTTATGACCCCCCCCTTTG‑ATGACCAAAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAA‑CACCCCA‑TCTTTCCGCGA‑CATTCTACCGCCTCTGAATTTCATCt                                                                                                               <  2:450644/150‑1 (MQ=255)
                                         cAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCttt                                                                                                             >  2:279728/1‑151 (MQ=255)
                                         cAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCttt                                                                                                             >  1:188118/1‑151 (MQ=255)
                                             aaaaaGCGCAACCCGGCCAGCAGCAGGGGGGGCATGTTAGGAAGCCCCTCTTTG‑AGGCCCAAAAAATTTAGCCCCCATCCGCCCCCTCCCAGTAGCCCCAA‑CCCCCCA‑TCTTTCCGCGA‑CTTTCTACCGCCTCTGAATTTCATCtt                                                                                                              <  1:1412026/144‑1 (MQ=255)
                                              tAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                         >  2:195721/1‑150 (MQ=255)
                                              tAAAGCCCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                         >  1:1539391/1‑150 (MQ=255)
                                               aaaGCGCAACCCGCCCAGCATCAGCGGTGGCATGTTATAAACCCCCATTTTG‑ATGCCCACAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAA‑CCCCCCA‑TCTTTCCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTg                                                                                                           <  2:108508/147‑1 (MQ=255)
                                               aaaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCCGAATTTCATCTTTTGTa                                                                                                         >  2:1503098/1‑149 (MQ=255)
                                                aaGCGAAAACCGCCCAGCATCAGCGGTGGCATTTTATGAACCCCCACTTTG‑ATGACCAAAAAATTTAGCCCCCATACGCCCCCTCCCAGTAGCGCCAC‑CCCCCCA‑TCTTTCCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                      <  1:64721/151‑1 (MQ=255)
                                                  gcgcAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCaa                                                                                                    <  1:609295/151‑1 (MQ=255)
                                                   cgcAAACCGGCCAGCATCAGCGGTGGCATTTTATAAAGCCCCTTTTTT‑ATCACCACAAAATTTACCCCCCATCCGCCCACTCCCAGTAGCGCCAA‑CCCCCCA‑TCTTTCCGCGA‑AATTCTCCCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                                     <  1:909228/149‑1 (MQ=255)
                                                     cAAACCGGTCATAATCAACGTTGATGTTTTAAGAGGCCCCCTTTTG‑AAGACAAAAAATTTAAGCCCCCATACGGCCACTCCGGGAAGCGACAA‑CCCCCCA‑TCTTTCGGAGA‑TTTTCTCCCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑Ca                                                                                                <  2:1561602/151‑1 (MQ=255)
                                                     cAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑c                                                                                                 >  1:39304/1‑150 (MQ=255)
                                                      aacccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCATTTGG‑ATGACCAAAAAATTTACCCCCCATACGCCCACTCCCAGTAGCCCCAC‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCaa                                                                                                    <  2:1032817/144‑1 (MQ=255)
                                                      aacccGGCCAACATCAGCGGTGGCATTTTATGAACCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTCCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑Ca                                                                                                <  2:1171402/147‑1 (MQ=255)
                                                      aacccGCCCACCATCAGCGGTGGTATGTTATGAAGCCCCTCTTTG‑ATCACCAAAAATTTGCCCCCCCATACGACCACCACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGC‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCaa                                                                                                    <  1:1285423/144‑1 (MQ=255)
                                                      aaaCCGGCCAGCTCCAGCGGTGGAATGTTATGAAGCCCAACTTGG‑ATGACAAAAAAATTTACCCCCCACACGCCCCCTACCAGTCGCGCAAA‑CACCCCC‑TCTTTCCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                                     <  2:742691/146‑1 (MQ=255)
                                                       acccGGCCACCATCAGCGGTGGCATGTTATGAAGCCCCATTTGG‑ATGACCACAAAATTTAGCCCCGATACGCCCACTACAAGTGGCCCCAA‑CCCCCCA‑TCTTTTCGAGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAAc                                                                                              <  2:1532804/149‑1 (MQ=255)
                                                         cgggCCAGCATCCGCGGGGGGATGTAATG‑ACCCCCCCTTTGAATAACCAAAAAATTAAGCCCCCACCCGACCCCAACCATTAGCG‑CAA‑CCCCCCATTCTTTTCGATA‑AATCCGACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACtt                                                                                            <  1:1563205/149‑1 (MQ=255)
                                                         ccGGCCAGCATCAGCGGTGGCTTGTTATGAAGCCCCACTTTG‑ATGCCCAAAAAATTTAGCCCCCATACGCCCCCTCCCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACtt                                                                                            <  2:754748/151‑1 (MQ=255)
                                                         ccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTACCCCCCATACGACCACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACtt                                                                                            <  2:1289745/151‑1 (MQ=255)
                                                          cGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACtt                                                                                            >  2:52106/1‑150 (MQ=255)
                                                           ggCCAGTATCAGCGGGGTAATGTTATGAAGCCCCACTTTG‑ATGACAAAAAAATTTACCCCCTATTCTCCCACTACCAGTCGCTCCAA‑CCCCCTA‑TTTTTTC‑CGACCATCCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑c                                                                                                 <  2:375985/144‑1 (MQ=255)
                                                           ggCCAGCATCAGCGGTGGCATGTTATAAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCAATACGCCCACTCCCAGTAGCGCCAA‑CCCCCCA‑CCTTTCCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTa                                                                                           <  1:701372/150‑1 (MQ=255)
                                                                 cATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAAt                                                                                    >  1:8315/1‑151 (MQ=255)
                                                                  aTCGGCGGGGGAATGTTATGAGGCCCCACTTTG‑ATGCCCAAAAATTTTACCCCCAAGACGCCCACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTCCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGc                                                                                         <  2:968241/145‑1 (MQ=255)
                                                                    cAGCGGTGGCATTTTATGAACCCCCCCTTTG‑ATGCCCACAAAATTTACCCCCCATCCGCCCCCTCCCAGTAGCCCCAC‑CCCCCCT‑TTTTTTCGCGA‑CTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTa                                                                                 <  1:156597/151‑1 (MQ=255)
                                                                    cAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAAttt                                                                                  >  2:212456/1‑150 (MQ=255)
                                                                      gCGGTGGCATGTTATGAACCCCCCCTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTAc                                                                                <  2:1278801/150‑1 (MQ=255)
                                                                        ggTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACttt                                                                             >  1:1182646/1‑151 (MQ=255)
                                                                         gTGGCATTTTATGAACCCCCCCTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACtttt                                                                            <  2:1132233/151‑1 (MQ=255)
                                                                          tGGCGTGTTAGGAAGCCCCACTTTG‑ATGCCCACAAAATTGACCCCCCATACGACCCCTACCAGAGGCGCCAA‑CCCCCCT‑TTTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACtttt                                                                            <  2:1415650/150‑1 (MQ=255)
                                                                          tGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGACTTTACTTTTc                                                                           >  1:1538917/1‑151 (MQ=255)
                                                                             cATTTTATGAAGCCCCTCTTTG‑ATGCCCAAAAATTTTAGCCCCCATACGCCCCCTCCCAGTAGCCCCAA‑CCCCCCA‑TTTTTCCGCGA‑CTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCtt                                                                         <  1:993492/150‑1 (MQ=255)
                                                                              aTGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCttt                                                                        >  2:1647672/1‑150 (MQ=255)
                                                                                 ttATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAC‑CACCCCA‑TCTTTTCGCGA‑CTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTaa                                                                      <  2:1360386/149‑1 (MQ=255)
                                                                                  tATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAAc                                                                     >  2:345331/1‑149 (MQ=255)
                                                                                   aTGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                                  >  1:1672313/1‑151 (MQ=255)
                                                                                   aTGAACCCCCACTTTG‑ATGACCACAAATTTTAGCCCCCATCCGACCACTCCCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAg                                                                   <  1:182600/150‑1 (MQ=255)
                                                                                    tGAACCCCCACTTTG‑ATGACCACAAAATTTACCCCCCATACGCCCCCTCCCAGTAGCGCCAA‑CCCCCCA‑TCTTTCCCCAA‑CTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                                  <  1:187409/150‑1 (MQ=255)
                                                                                     taaGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                                >  1:1194395/2‑151 (MQ=255)
                                                                                      aaGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAA‑CCCCCCA‑TCTTTCCCCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                                 <  2:670187/149‑1 (MQ=255)
                                                                                         ccccATTTTA‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGAAGCGCCAA‑CCCCCCT‑TTTTTTCGCAA‑ATTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAt                                                              <  2:1103435/149‑1 (MQ=255)
                                                                                         ccccACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                            >  2:104427/1‑151 (MQ=255)
                                                                                         ccccACTTTG‑ATGACCAAAAAATTTACCCCCAATCCGACCACTCCCAGTATCCCCAA‑CCCCCCA‑TTTTTGCCCGA‑ATTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAt                                                              <  1:693917/149‑1 (MQ=255)
                                                                                          cccACTTTG‑ATGACCACAAAATTTAGCCCCCATCCGACCACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                           <  2:168911/151‑1 (MQ=255)
                                                                                          cccACTTTG‑ATGACCACAAAATTTAGCCCCCATACGCCCCCTCCCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                           <  1:231110/151‑1 (MQ=255)
                                                                                          cccACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                           <  2:898569/151‑1 (MQ=255)
                                                                                              cTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGtt                                                         >  1:1574433/1‑149 (MQ=255)
                                                                                                ttG‑ATGCCCAAAAAATTTACCCCCCATCCGCCCACTCCCAGTAGCCCCAA‑CCCCCCA‑TCTTTTCGCGC‑CTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTa                                                        <  1:1151452/148‑1 (MQ=255)
                                                                                                ttG‑AGGCCCACAAATTTTACCCCCCATACGCCCCCTACCTGTAGCGCCAA‑CACCCCT‑TCTTTCCCCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGt                                                      <  1:405448/150‑1 (MQ=255)
                                                                                                  g‑atgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTC‑CGG‑CATTCCACCCCCTCCGGATTTCATCCTTTTGAAACAATCCACCTTAGTTG‑ATTTTCTTTTTTTTAAACGGTGGTTCGTTAAGCGcc                                                  >  1:1201357/1‑151 (MQ=255)
                                                                                                         cacaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                              >  2:720276/1‑150 (MQ=255)
                                                                                                         cacaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                              >  2:838258/1‑150 (MQ=255)
                                                                                                         cacaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                              >  2:847146/1‑150 (MQ=255)
                                                                                                            aaaaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacg                                           >  2:115915/1‑150 (MQ=255)
                                                                                                               aTTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                       >  1:220329/1‑151 (MQ=255)
                                                                                                                tttAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                                      >  1:1022785/1‑151 (MQ=255)
                                                                                                                 ttACCCCCCATACGACCACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                                      <  2:1198256/150‑1 (MQ=255)
                                                                                                                  tAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAt                                    >  1:463009/1‑151 (MQ=255)
                                                                                                                   aGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                   >  1:887838/1‑151 (MQ=255)
                                                                                                                     cccccATACGACCACTACCAGTAGCCCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTa                                  <  2:1074601/150‑1 (MQ=255)
                                                                                                                      ccccacccaaCCATTCCCAGTAG‑GCCAACCACCCCA‑TTTTTTCGCCA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                                <  1:1616201/142‑1 (MQ=255)
                                                                                                                      ccccAAACGACCTCTACCAGTAGCGCCAC‑CTCCCCA‑TCTTTTCGCTA‑CATCCTTCCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                                <  2:626095/151‑1 (MQ=255)
                                                                                                                       cccATACGACCACTCCCGTTAGCGCCAC‑CCCCCCA‑TTTTTCCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATg                               <  2:943049/151‑1 (MQ=255)
                                                                                                                        ccATACGCCCCCTCCCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCAA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                              <  1:688385/151‑1 (MQ=255)
                                                                                                                            aCGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCa                           >  1:327023/1‑150 (MQ=255)
                                                                                                                             cGACCACTCCCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcgc                            <  2:1175322/148‑1 (MQ=255)
                                                                                                                              gACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                        >  1:1687350/1‑151 (MQ=255)
                                                                                                                                ccccTACCAGTAGCCCAAC‑CCCCCCC‑TCTTTTCGCGA‑AATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                      <  2:244111/148‑1 (MQ=255)
                                                                                                                                 cACTACCAGTAGCGCCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAg                     <  1:640383/151‑1 (MQ=255)
                                                                                                                                   cTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc                   >  2:294947/1‑151 (MQ=255)
                                                                                                                                    tACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg                  >  1:331210/1‑151 (MQ=255)
                                                                                                                                     aCCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGc                 >  2:342938/1‑151 (MQ=255)
                                                                                                                                            gcgcCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACtt          >  2:845945/1‑151 (MQ=255)
                                                                                                                                             cgcCAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACt           <  1:1429179/149‑1 (MQ=255)
                                                                                                                                             cgcCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc         >  1:521632/1‑151 (MQ=255)
                                                                                                                                               ccAA‑CCCCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg       <  1:1388217/151‑1 (MQ=255)
                                                                                                                                               ccAA‑CCCCCCA‑TCTTTCCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg       <  2:842026/151‑1 (MQ=255)
                                                                                                                                                cAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg       >  1:743811/1‑150 (MQ=255)
                                                                                                                                                    cacCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAttt   >  2:472773/1‑151 (MQ=255)
                                                                                                                                                     acCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtttt  >  2:179873/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                         
ACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTG‑ATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAA‑CACCCCA‑TCTTTTCGCGA‑CATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAT‑CAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620971‑1621268

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: ‑