Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 17.4% intergenic (+112/+261) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,1770AC17.4% 257.9 / 5.9 148intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (28/75);  new base C (16/9);  total (52/96)
Fisher's exact test for biased strand distribution p-value = 8.70e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑2229326
                                                                                                                                                      |                                                                                                                                                        
     tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑G‑AGGGggc                                                                                                                                                    >  1:1030792/1‑150 (MQ=255)
      ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGggcg                                                                                                                                                   <  1:2606975/150‑1 (MQ=255)
       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGggcgg                                                                                                                                                  <  2:147120/150‑1 (MQ=255)
        tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑G‑AGGGgggggc                                                                                                                                                 >  2:1571500/1‑150 (MQ=255)
         gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGcc                                                                                                                                                <  2:3100198/150‑1 (MQ=255)
         gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGcc                                                                                                                                                <  2:879044/150‑1 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCa                                                                                                                                               <  1:1503176/150‑1 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCG‑C‑CCGGGGCGGcccac                                                                                                                                             >  1:3532158/1‑147 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑g‑ggggggggggcgag                                                                                                                                             >  1:3308449/1‑145 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑CGGGGGCGGcccaa                                                                                                                                             >  2:2995193/1‑147 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCCAt                                                                                                                                             >  1:941069/1‑150 (MQ=255)
                gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATcgcc                                                                                                                                         <  1:1653182/150‑1 (MQ=255)
                   tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCG‑C‑GGGGGGCGGCCCCAcgccgcc                                                                                                                                      >  1:3540306/1‑150 (MQ=255)
                   tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCC‑CGGGGGGGGCCAACcgcgcc                                                                                                                                       >  2:4172775/1‑147 (MQ=255)
                     tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATcgccgccgc                                                                                                                                    <  1:4103283/150‑1 (MQ=255)
                     tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATcgccgccgc                                                                                                                                    <  1:3941637/150‑1 (MQ=255)
                       aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑CGGGGGGGGCCAAAAGCcccccccc                                                                                                                                  >  2:2416831/1‑150 (MQ=255)
                       aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCCATAGCCCCCGccc                                                                                                                                  >  2:2417561/1‑150 (MQ=255)
                       aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGccc                                                                                                                                  <  2:1296906/150‑1 (MQ=255)
                          tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑G‑GGGGGGGGGGCCAAAGCccccccccccc                                                                                                                               >  2:4076834/1‑149 (MQ=255)
                          tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑G‑ACGGGGGGGGCAATAGCCGccccccccc                                                                                                                               >  1:3851388/1‑149 (MQ=255)
                              atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑ACGGGGCGGCCAAAAGCCGcccccccccccac                                                                                                                           >  2:568811/1‑145 (MQ=38)
                               ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑CGGGGGGGGGGCCAAACccccccccccccccccc                                                                                                                           >  2:865035/1‑145 (MQ=255)
                                 ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑CAGGGGCGGCCCATAGCCGCCCCCCCCCGcccccc                                                                                                                        >  2:2580036/1‑150 (MQ=38)
                                     aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                    <  2:1945919/150‑1 (MQ=255)
                                     aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGATTCATATAATCACTGGAGCAAGGAGCGCG‑C‑CTGGAGAGGGCAATGGCCGCCCACCCCCGCACACCCcgg                                                                                                                    >  1:2571271/1‑150 (MQ=255)
                                           ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑G‑AGGGGGGCGCCAAACGCCCCCCCCCCCAGCAACCCCCGGCACTgg                                                                                                              >  1:2447014/1‑150 (MQ=255)
                                              atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGa                                                                                                          <  2:1029045/150‑1 (MQ=255)
                                                   atcaGGACCTAAAATAACGATAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGATTAATATAATCACTGGCGAATGGAGCGCG‑C‑CGGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑aaa                                                                                                    <  1:2389808/150‑1 (MQ=25)
                                                    tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑G‑GGGGGGGGGGCAAAACCCCCCGCCCCCACCACCCCCCCGCGCTGG‑CGGAC‑aaaa                                                                                                   >  1:2168875/1‑150 (MQ=25)
                                                    tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑CGGGGGGGGCCCAAAGCCCCCCCCCCCCGCCCCCCCCCGCGCGGG‑CGCAC‑aaaa                                                                                                   >  2:4178616/1‑150 (MQ=25)
                                                    tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑CGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCGGGCTCTGG‑CGAAC‑aaaa                                                                                                   >  2:2840801/1‑150 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑CCGGGGCGGGCCCACCCCCCCCCCCCCCCCCCCCCCCGCCCCGGG‑CGAaa‑aaaaaa                                                                                                 >  2:564815/1‑148 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCCGGCTCTGG‑CGAAC‑AAAATc                                                                                                 >  2:2877097/1‑150 (MQ=255)
                                                      aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTGCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAAACGCCGCCGCCCCCGGCCCCCCCCGGCGCGGG‑CGCAC‑AAAATc                                                                                                 >  2:2479207/1‑150 (MQ=25)
                                                        gTCCGAAAATACCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATcgc                                                                                               <  2:1238936/150‑1 (MQ=255)
                                                         tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATcgcc                                                                                              <  2:3922399/150‑1 (MQ=255)
                                                          ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCAGG‑CGAAC‑AAAATcgccg                                                                                             >  2:2915174/1‑150 (MQ=255)
                                                          ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑ACGGGGCGGCCAAAACCCCCCCCCCCCAGCCCCCCCCCGCGCAGG‑CGCAC‑acaaaacccc                                                                                             >  2:4214856/1‑144 (MQ=25)
                                                                aTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTgcgc                                                                                       <  2:2253379/150‑1 (MQ=255)
                                                                aTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTgcgc                                                                                       <  2:2253353/150‑1 (MQ=255)
                                                                   aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑G‑GGGGGGGGGCCCCACAGCGCCGGCCCCCGCCCCCCCCCCGGCGGG‑GGGCCAAAAATCGCCGCTGCGCgg                                                                                     >  2:1858132/1‑150 (MQ=255)
                                                                      agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGcc                                                                                 <  2:2299499/150‑1 (MQ=255)
                                                                                  tcgctctTCCTGTTGCGCTCCTGTCATGGGAATTGCTACATATAACCCCTGGCGCACGGCGCGCG‑C‑GGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                     <  1:1710032/147‑1 (MQ=14)
                                                                                    tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCCATCCCCCCCCCCCCCACCACCCCCCGGCTCAGG‑CCAAC‑AAAAAACCCCCAGCGCGGTGCCCACCGCTTATacca                                                                   >  2:3588962/1‑149 (MQ=18)
                                                                                      tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCCATACCCGCCCCCCCCACCAACCCCCGGCTCATG‑CGAAA‑AAAAAACCCCCAGCCCGGTGCCCACCGCAGATCCCATa                                                                 >  2:2133803/1‑150 (MQ=12)
                                                                                      tcACCCTGTGGCGCTCCTGACATGTGCTTTGCTTCATATAATCACTGGCGCACGGCGCGCG‑C‑CGGGGGCGCCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                 <  1:1428158/150‑1 (MQ=17)
                                                                                      ccATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑CGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                 <  1:2555703/149‑1 (MQ=255)
                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                               <  2:1344555/150‑1 (MQ=255)
                                                                                            tgtgGCGCTCGTGTCTTGTGCTTTGCTTCATATAATCAGTGGCGCATGGTGCGCG‑C‑CGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑a                                                          <  1:934745/150‑1 (MQ=17)
                                                                                            tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCG‑C‑TGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑a                                                          <  1:3147106/150‑1 (MQ=255)
                                                                                            tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑a                                                          <  2:2601029/150‑1 (MQ=255)
                                                                                            tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑a                                                          <  2:2954871/150‑1 (MQ=255)
                                                                                                 ggcTCCTGTCTTGTGCATTGCTTCATATAATCCCTGGCGCACGGTGCGGG‑C‑CGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCggg                                                     <  1:776710/149‑1 (MQ=17)
                                                                                                 cgcTCCTGTCTTGTGCTTTGCTTCATATAATCACTGGCGCATGGAGCGCG‑C‑TGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCggg                                                     <  1:776623/150‑1 (MQ=255)
                                                                                                   cTCCTGTCTGGTGCTTTGCTTCATATAATCAGTGGCGCACGGAGCGCG‑C‑TGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTc                                                   <  1:1853015/150‑1 (MQ=17)
                                                                                                   cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCATGGAGCGCG‑C‑TGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTc                                                   <  1:3469719/150‑1 (MQ=255)
                                                                                                    tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑CCGGGGCGGCCAAAAGCCGCCCCCCCCCGCACCCCCCGGCTCGGGCCGAAA‑AAAAAC‑CCCCAGCGCGGTGCCCACCGCGTATACCATACGGCTAACCgggcg                                                   >  2:1752358/1‑148 (MQ=12)
                                                                                                    tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATAGCCGCCGCCCCCTGCACCCCCCGGCTCTGG‑CGAAC‑AAAAACGCCCCGGCGCGGTGCCCTCCGCGTATCCCATACGGCT‑ACCCGGTCg                                                  >  2:2912175/1‑150 (MQ=17)
                                                                                                    tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATAGCCGCCGCCCCCAGCACCCCCGGGCTCTGG‑CGAAA‑AAAATAGCCGCAGCGCGGTGCCCACGGCTTATCCCATAAGGCT‑AACCGGTCg                                                  >  2:735810/1‑150 (MQ=17)
                                                                                                     ccTGTCTTGTGCTTTGCTTAATATAATCACTGGCGCACGGAGCGCG‑C‑TGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCgg                                                 <  1:1906497/150‑1 (MQ=14)
                                                                                                     ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCTCCGCTGCGCGGTGCCCTCGGATTATCCCTTACGGCT‑ACCGGGTCgg                                                 <  2:2365088/150‑1 (MQ=255)
                                                                                                     ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCgg                                                 <  2:431194/150‑1 (MQ=255)
                                                                                                      cTGTCTTGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCggg                                                <  1:3033713/150‑1 (MQ=255)
                                                                                                      cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCggg                                                <  2:3258075/150‑1 (MQ=255)
                                                                                                      cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCggg                                                <  2:3532158/150‑1 (MQ=255)
                                                                                                      cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCggg                                                <  2:941069/150‑1 (MQ=255)
                                                                                                      cTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCggg                                                <  2:3428879/150‑1 (MQ=255)
                                                                                                         tCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCCATCGCCCCCGCCCCCAGCACCCCCCGGCGCTGG‑CCAAC‑AAAATCCCCGCAGCCCGGGGCCCACGGCGTATACCTTACCGCT‑ACCCGGTCGGgcgc                                             >  2:1041871/1‑150 (MQ=17)
                                                                                                         tCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCCAAAGCCCCCCCCCCCAGCACCCCCCGGCGCTGG‑CGGAC‑AAAAACGCCCCGGCGCGGTGGCCACCGCTTATCCCGTAAGGCG‑AACCGGTCGGgcgc                                             >  2:357048/1‑150 (MQ=12)
                                                                                                          cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGcgcg                                            <  1:345378/150‑1 (MQ=255)
                                                                                                           aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGa                                           <  2:3055392/150‑1 (MQ=255)
                                                                                                            tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAg                                          <  1:2417561/150‑1 (MQ=255)
                                                                                                            tgtgCATTGCTTCATATAATCACGGGCGCAAGGAGCGCG‑C‑CGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAg                                          <  1:4178616/150‑1 (MQ=255)
                                                                                                                  ttGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCATTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACa                                    <  2:4434/150‑1 (MQ=255)
                                                                                                                    gCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATc                                  <  2:2830478/150‑1 (MQ=255)
                                                                                                                      ttCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATccc                                <  2:3443525/150‑1 (MQ=40)
                                                                                                                       tCATATACTCCCTGGCGCACGGAGCGGG‑C‑CGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:4201573/150‑1 (MQ=11)
                                                                                                                         atatAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGt                             <  2:3366036/150‑1 (MQ=38)
                                                                                                                         atatAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGt                             <  2:3366098/150‑1 (MQ=38)
                                                                                                                             aaTCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:3974748/150‑1 (MQ=34)
                                                                                                                                cACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:4207488/150‑1 (MQ=32)
                                                                                                                                 aCTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  2:3130511/150‑1 (MQ=32)
                                                                                                                                  cTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:3167296/150‑1 (MQ=32)
                                                                                                                                  cTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:2571913/150‑1 (MQ=32)
                                                                                                                                  cTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:3994032/150‑1 (MQ=32)
                                                                                                                                       gcAAGGAGCGCG‑C‑AGGGGGCGGCCAATCCCCCCCCCCCCCAGCACCCCCCGGCTCTGG‑CGAAC‑AAAATAGCCGCGGCGCGGTGCCCTCGGCGTATCCCATACCGCT‑AACCGGTCGGGCGCGAGGTAACAACCCAGTAAAACCCGCCCGCa               >  2:3960275/1‑150 (MQ=11)
                                                                                                                                        cAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:2627440/150‑1 (MQ=32)
                                                                                                                                        cAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  2:3499643/150‑1 (MQ=32)
                                                                                                                                          aGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  2:1455886/150‑1 (MQ=32)
                                                                                                                                          aGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:4214856/150‑1 (MQ=32)
                                                                                                                                            gAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca          <  1:1201089/150‑1 (MQ=32)
                                                                                                                                               cgcg‑c‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcacac       <  1:480486/150‑2 (MQ=17)
                                                                                                                                                gcg‑c‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc      >  2:663986/1‑150 (MQ=32)
                                                                                                                                                gcg‑c‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc      >  2:1977749/1‑150 (MQ=32)
                                                                                                                                                gcg‑c‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc      >  2:2719489/1‑150 (MQ=32)
                                                                                                                                                 cg‑c‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  1:3537400/150‑1 (MQ=32)
                                                                                                                                                      aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  <  2:3989722/150‑1 (MQ=32)
                                                                                                                                                      aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  <  1:3497083/150‑1 (MQ=32)
                                                                                                                                                      aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  <  2:3989600/150‑1 (MQ=32)
                                                                                                                                                      |                                                                                                                                                        
ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCG‑C‑AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAC‑AAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCT‑ACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑2229326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑