Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,026,9410AC40.5% 126.0 / 16.2 79E392D (GAA→GAC) guaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (37/10);  new base C (0/32);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.19e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAA  >  NC_000913/3026799‑3027088
                                                                                                                                                     |                                                                                                                                                   
aaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGc                                                                                                                                             <  2:296266/151‑1 (MQ=255)
 aGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGcc                                                                                                                                            >  2:295173/1‑151 (MQ=255)
  gTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGcc                                                                                                                                            >  2:419859/1‑150 (MQ=255)
          aTTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCCGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCt                                                                                                                                    >  2:392377/1‑150 (MQ=255)
                aggctccccccTCTCGGCATATGAAGGGTTTTCCCGGGCC‑CCGCGGGGCGGAGCAAAATCTCTGGGCTTTGA‑CGTTTTGATGGCCAAT‑TTTTTCCCGGGA‑AAA‑GAGT‑C‑GATTTGCGTGGTGAGGGACCCCACCGCCACTCCGCTAca                                                                                                                                 <  1:762665/138‑1 (MQ=255)
                aGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGCGGTGATGGAACCCACCGCCACTCCGCTAcagc                                                                                                                               >  2:798374/1‑149 (MQ=255)
                  gCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCt                                                                                                                           >  1:331381/1‑151 (MQ=255)
                     accGCCTCTCGGCATATGAAGCTTTTCCCCTGGCC‑CCGCCCGGCGGAGCAAAATCTCTGGGCCTTGA‑CTATTTGATGGGCAAT‑TTTTTCCCGGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTgcg                                                                                                                        <  2:670012/149‑1 (MQ=255)
                       cGCCTCTCGGCATATGAAGCTTTTTCCCTGGCC‑ACGCTCGGCGGAGCAAAATCTCGGGGCTTTGA‑CGTTTTGATTGGCACT‑TTTTTCCCGGGA‑AAA‑GAGG‑CTGTTTT‑CGGGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCt                                                                                                                           <  2:531119/146‑1 (MQ=255)
                       cGCCTCTCGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCt                                                                                                                      >  1:541733/1‑151 (MQ=255)
                       cGCCCCTCGGCATATGAAGCTTTTTCCCTGCCC‑ACGCTCGGCGGAGCAAATTCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCGGGA‑AAA‑GAGG‑TTTATTT‑GGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCt                                                                                                                      <  2:853946/151‑1 (MQ=255)
                         ccccccGGCATAGGAAGCGTTTTACCGGGCC‑CCGCTCGGCGGAGCGAAAGCTCTGGGCCTTGA‑CGTTTTGATTGGCAAC‑TTTTTCCCTGGA‑AAA‑GAGG‑CGGATTT‑CTGGGTGAGGAAACCCCCCGCCACTCCGCTACAGCAGCTGCGCTAt                                                                                                                    <  2:1189859/146‑1 (MQ=255)
                         ccTCTCGGCATATGAAGCTTTTTACCTGGCC‑ACGCTGGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTCCCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTAt                                                                                                                    <  2:289624/151‑1 (MQ=255)
                             tcGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGAc                                                                                                                 >  1:1180462/1‑150 (MQ=255)
                               ggCATAGGAGGGTTTTTCCCTGGCC‑CCCCGGGGGGAAGCAAAATCTCGGGGCTTTGA‑CGATTTGATGGGCAAC‑TTTTCCCCTGGA‑AAG‑GAGG‑CTGTTTG‑GGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                              <  1:787527/151‑1 (MQ=255)
                               gctatatGAGGCG‑TTTACCCGGCCGCCGCGGGGGGGAGCAAAATCTCTGGGCCTTGG‑ATATTTAATT‑GCAACTTTTTTTCCTGAA‑AAA‑GAGT‑GTGATTT‑GGTGGGGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGa                                                                                                                  <  1:804945/144‑1 (MQ=255)
                                 cATATGAACGTTTTTACCGCGCC‑CCGCTCGGGGGAGAGAAATCTCTGGACTTGAA‑CATTTTGATGAGAATC‑TTTTTCCCAGGC‑AAA‑GAGG‑TTGTTTT‑CGTGTTGATGAAACCGGCCGCCACTCCGCTACAGCAGCTGCGCTATGACAAct                                                                                                             <  2:156871/150‑1 (MQ=255)
                                  atatGAGGCTTTTTCCCTGCCC‑CCGCTCGGGGGGGCAAAATCTCGGGCCCTTGA‑CGATTTGATGGGCATT‑TTTTTCCCTGGA‑AAA‑GAGG‑TTGTTTT‑GGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActct                                                                                                           <  1:446662/151‑1 (MQ=255)
                                  atatGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGACTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActct                                                                                                           >  2:635237/1‑151 (MQ=255)
                                            ttttCCCTGCCC‑ACGCCCGGGGGAGCAAAATCCCGGGGCTTTGA‑CGATTTGATTGGCATC‑TTTTTACCGGAC‑AAA‑GAGG‑TTGTTTT‑CGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCt                                                                                                     <  1:628971/147‑1 (MQ=255)
                                            ttttACCTGGCC‑ACGCTCGGCGGGGGGAATTCTCGGGGCTTTGA‑CGATTTTATTGGAATC‑TTTTTACCTGGC‑AAA‑GAGT‑TTGATTA‑GGGGGTAATGAAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTc                                                                                                      <  2:689856/146‑1 (MQ=255)
                                            ttttACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                                   >  2:175339/1‑149 (MQ=255)
                                            ttttACCGGCCC‑ACGCTGGGCGGGGTGAATTCTCTGGTCATTGA‑GTATTTAATGGGCATT‑TTTTTGCCTGGA‑AAG‑GGGG‑CTTATTT‑GGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTc                                                                                                      <  2:651411/146‑1 (MQ=255)
                                             tttCCCTGGCC‑CCGCTCGGCGGAGCGAAATCTCGGGGCCTTGA‑CGTTTTGATGGGCAAC‑TTTTTCCCTGGC‑AAA‑GAGG‑CTGTTTT‑CGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAg                                                                                                 <  2:301412/150‑1 (MQ=255)
                                             ttacccTGCCC‑ACGCTCGGGGGAGAAAATTTTTGGGGCCTGGA‑ATATTTGATTGGATAT‑TTTTTCCTTGAA‑AAA‑GAGT‑TTGTTTT‑GGTGGTGAAGGAACCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAg                                                                                                 <  2:517746/146‑1 (MQ=255)
                                               tACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTc                                                                                               >  2:1183786/1‑150 (MQ=255)
                                                 ccTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGAc                                                                                            >  2:84343/1‑151 (MQ=255)
                                                           gCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc                                                                                   >  2:302181/1‑151 (MQ=255)
                                                            ctggggggAGCAAAATCTCTGGCCCTTGA‑CGATTTGATTGGCACC‑TTTTTCCCGGGC‑AAA‑GGGG‑TTGATTT‑GGTGGTAATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt                                                                                    <  1:302835/143‑1 (MQ=255)
                                                               gggggAGCAAAATCTCGGGCCCTTTG‑TGATTTGATTGCACAC‑TTTTTCCCTGGC‑AAA‑GAGG‑CTGTTTT‑TGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGt                                                                                     <  1:776569/142‑1 (MQ=255)
                                                               ggcggAGCAAATTCTCTGGCCTTTGA‑CTTTTTGATGGACATT‑TTTTTCCCGGAA‑AAA‑GGGG‑TTGTTTG‑CGGGGTGAGGCACCCAACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtga                                                                               <  1:795679/151‑1 (MQ=255)
                                                                  ggAGCAAACTCTC‑GGGCCTTTAGAGATTTGTTGGTCAAC‑TTTTTACCGGGA‑AAA‑GAGC‑GTAATTT‑TGTTGTGAGGAAACCCGCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatga                                                                            <  2:883702/151‑1 (MQ=255)
                                                                    aGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg                                                                          >  2:460125/1‑151 (MQ=255)
                                                                        aaaTCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGt                                                                      >  1:116287/1‑151 (MQ=255)
                                                                        aaaTCTCGGGGCTTTGA‑CGATTTGATTGGCAAT‑TTTTTCCCTGGA‑AAG‑GGGG‑TTGTTTT‑GGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg                                                                          <  1:480430/147‑1 (MQ=255)
                                                                         aaTCTCTGGGCCTTGA‑CGTTTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGTTTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGtt                                                                     <  2:271907/151‑1 (MQ=255)
                                                                         aaTCTCGGGGCTTTGA‑CGATTTGATGGGCATC‑TTTTTCCCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGAc                                                                        <  2:27772/148‑1 (MQ=255)
                                                                          aTCTCTGGGCCTTGA‑CGATTTGATTGGCACT‑TTTTTCCCTGGA‑AAA‑GGGG‑TTGATTT‑CGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGtt                                                                     <  2:161433/150‑1 (MQ=255)
                                                                            ctctGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTggg                                                                  >  1:569922/1‑151 (MQ=255)
                                                                               tGGTCCTTGG‑CGATGTGTTGGCCACC‑TTTTTACTGGCC‑AAA‑GCGT‑GTGATTT‑CTGGTTGATGGACCCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTgg                                                                   <  1:116918/147‑1 (MQ=255)
                                                                               gggcccTGGA‑CGATTTGTTTGGC‑AC‑TTTTTCACCGGC‑AAAGGAGG‑TTGATTT‑GGTGGTGGTGCACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGa                                                               <  1:870132/147‑1 (MQ=255)
                                                                                gggCTTTGA‑CTTTTTGATGGGCATC‑TTTTCACCTGAA‑AAA‑GGGG‑CTGTTTT‑CGGGGTGTGGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGAt                                                              <  1:1238311/151‑1 (MQ=255)
                                                                                 gggcctgaccattttGATTGCCATT‑TTTTCCCCTGCA‑AAG‑GTGG‑TTGTTTT‑GGGGGTGATGGCCCCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTg                                                                    <  1:1061757/133‑1 (MQ=255)
                                                                                     ttGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCg                                                         >  1:127261/1‑151 (MQ=255)
                                                                                       gA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGtt                                                       >  2:1084059/1‑151 (MQ=255)
                                                                                            aTTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGAt                                                   >  2:65051/1‑151 (MQ=255)
                                                                                                  ttGGCAAT‑TTTTTCCCGGGA‑AAA‑GGGG‑CTAATTG‑CGGGGTAAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTAcc                                              <  2:833183/150‑1 (MQ=255)
                                                                                                   tGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGc                                            >  2:633413/1‑151 (MQ=255)
                                                                                                    gacAAC‑TTTTTACCTGGC‑AAA‑GAGC‑CTGTTTT‑CGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGc                                            <  1:356537/148‑1 (MQ=255)
                                                                                                       aaC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGTCAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCt                                        >  1:112394/1‑151 (MQ=255)
                                                                                                         accttttTCCGGGCAAAA‑GGGG‑CTGTTTT‑GGTGGGGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCa                                           <  1:1146470/145‑1 (MQ=255)
                                                                                                         aatttttCCCCGGC‑AAA‑AAGGTCTGATTT‑GGGGGTTTTGGACCCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCt                                        <  2:236271/149‑1 (MQ=255)
                                                                                                            ttttACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTTGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGt                                    >  1:567380/1‑151 (MQ=255)
                                                                                                              ttACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGtt                                   >  1:1030836/1‑150 (MQ=255)
                                                                                                               tACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTg                                  >  2:352256/1‑150 (MQ=255)
                                                                                                                 ccTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGAt                                >  1:239406/1‑150 (MQ=255)
                                                                                                                 ccTGGA‑AAA‑GGGT‑GTGATTT‑GGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATg                               <  1:1210638/151‑1 (MQ=255)
                                                                                                                 ccTGAC‑AAA‑GAGG‑CTGTTTT‑CGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATg                               <  2:513766/151‑1 (MQ=255)
                                                                                                                   tGGC‑AAA‑GAGT‑CTTTTTT‑CGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATgg                              <  1:1112251/150‑1 (MQ=255)
                                                                                                                   tGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATg                               >  1:219495/1‑149 (MQ=255)
                                                                                                                    ggC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACACATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCCCCTACGTTGATGgt                             >  1:290511/1‑150 (MQ=255)
                                                                                                                    ggC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtc                            >  2:881108/1‑151 (MQ=255)
                                                                                                                     ggaaaa‑GAGG‑CTGTTTT‑CGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATgg                              <  1:457613/146‑1 (MQ=255)
                                                                                                                     gC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcg                           >  2:1106958/1‑151 (MQ=255)
                                                                                                                      c‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcg                           >  1:1101618/1‑150 (MQ=255)
                                                                                                                        aaa‑GGGG‑TTGTTTT‑GGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgt                          <  2:293493/150‑1 (MQ=255)
                                                                                                                        aaa‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgtc                         >  2:331580/1‑151 (MQ=255)
                                                                                                                         aa‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCt                        >  1:1025923/1‑151 (MQ=255)
                                                                                                                          a‑gagG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTg                       >  2:1111447/1‑151 (MQ=255)
                                                                                                                                 cttaggt‑CGTGGTGGGGAACCCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTTCCGCACCTACGTTGATGTTCGTCTGgt                     <  1:632919/142‑1 (MQ=255)
                                                                                                                                 cTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtgt                   >  1:1000317/1‑150 (MQ=255)
                                                                                                                                  tGATTT‑CGTGGTGATGAAACCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGgtg                    <  2:850426/148‑1 (MQ=255)
                                                                                                                                          gtggtgATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAAc             <  2:1227389/148‑1 (MQ=255)
                                                                                                                                              tgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACt       >  2:362379/1‑150 (MQ=255)
                                                                                                                                              tgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACt       >  1:279330/1‑150 (MQ=255)
                                                                                                                                               gATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACta      >  2:1040609/1‑150 (MQ=255)
                                                                                                                                                aTGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaa     <  2:1187063/150‑1 (MQ=255)
                                                                                                                                                  ggAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaataa  <  2:1058978/151‑1 (MQ=255)
                                                                                                                                                   gAAACCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACtaata   >  1:262876/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCC‑ACGCTCGGCGGAGCGAAATCTCTGGGCCTTGA‑CGATTTGATTGGCAAC‑TTTTTACCTGGC‑AAA‑GAGG‑CTGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAA  >  NC_000913/3026799‑3027088

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