Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,322 C→G 5.6% intergenic (+257/+116) 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,3220CG 5.6% 267.2 / 8.3 105intergenic (+257/+116)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (56/40);  new base G (2/4);  total (61/44)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 107 total aligned reads displayed.

GCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAGCC  >  NC_000913/2229173‑2229471
                                                                                                                                                     |                                                                                                                                                     
gcgcAGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCGCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc                                                                                                                                                       <  2:1410612/150‑1 (MQ=17)
  gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa                                                                                                                                                     <  2:489199/150‑1 (MQ=32)
     gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACGTCCCTGTACACCGCGCCCTCAGCCCACATCCAtgt                                                                                                                                                  >  2:1464518/1‑150 (MQ=12)
      ggggCGGCCAACCGCGGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg                                                                                                                                                 <  1:1130525/150‑1 (MQ=12)
        ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCCGCCCCCCTCCATGTggg                                                                                                                                               >  1:725726/1‑150 (MQ=12)
         gcggcCAATCGCCGCCGCCCCCTGCCCCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGGAACCTCCCTGGAAAACGCGCCCTCCGCCCCCATCCCTGTGGGc                                                                                                                                              >  1:1302455/1‑150 (MQ=11)
         gcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCGCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCGCAGCCCACATCCATGTGGGc                                                                                                                                              <  2:293342/150‑1 (MQ=17)
         gcggcCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCGCAGCCCACATCCATGTGGGc                                                                                                                                              <  2:549565/150‑1 (MQ=17)
           gtccAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGGATGTGGGCTg                                                                                                                                            <  1:1115876/148‑1 (MQ=12)
           ggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCGCTGGCGACCAAAATCGCCGCTGCGCGGTGCCCTCGGCTTCTCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGGATGTGGGCTg                                                                                                                                            <  1:267372/150‑1 (MQ=11)
                aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCg                                                                                                                                       >  1:1224635/1‑150 (MQ=32)
                 tCGCCGCCGCCCCCGGCACCCCCGGGCTCGGGCGCCCACCAGCGCCGCGGCGCGGTGCCCTCGGCTTATCCCGTACGGCTACCGGGTCGGGCGCGAGGGAACATCCCTGTACAACGCGCCCGCAGCCCACATGGATGTGGGCTgagggcg                                                                                                                                      <  1:100017/150‑7 (MQ=11)
                   gccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAACACGCGCCCTCCGCCCCCATCCATGTGGGCTGCCCCGGcc                                                                                                                                    >  1:1185834/1‑150 (MQ=12)
                   gccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGcc                                                                                                                                    <  2:779821/150‑1 (MQ=32)
                        ggccCCCTGCACCCCCGGGCTCTGGCGAACAAACTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAg                                                                                                                               <  2:1856852/149‑1 (MQ=17)
                         gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAgg                                                                                                                              <  2:72878/150‑1 (MQ=32)
                              cTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACCTCCCTGTAATACGCGCCCTAAGCCCACATCCATGTGGGCTGCCCCGGCCTGAAGGGaccc                                                                                                                         >  1:1785499/1‑147 (MQ=11)
                               tGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAGGTGGGCTGCCCCGGCCGTCAGGGAACGc                                                                                                                        >  1:936930/1‑150 (MQ=17)
                               tGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACAGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATACATGTGGGCTGCCCCGGCCTTCAGGGAACGc                                                                                                                        >  2:1255594/1‑150 (MQ=17)
                                gCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCCCAGCCATGTGGGCTGCCCCGCCCGGCAGGGCACGcc                                                                                                                       >  1:950818/1‑150 (MQ=11)
                                gCACCCCCGGGCGCTGGCGAACAAATTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGcc                                                                                                                       <  2:365803/150‑1 (MQ=17)
                                   cccccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTTCGGCTACCGGGTCGGGCGCGGGGGAACATCCCTGGAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGGCGTCAGGGGACGCCTCg                                                                                                                    >  1:873533/1‑150 (MQ=11)
                                   cccccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCCCCCCACATCCATGTGGGCTGACCCGGCCTTCAGCGAACCCcgcg                                                                                                                    >  2:1103557/1‑150 (MQ=11)
                                         ggCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCCCATCCATGTGGGCTGCCCCGGCCGTCAGGGAACGCCTCGGCAAtt                                                                                                              >  1:268775/1‑150 (MQ=21)
                                          gCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAAttt                                                                                                             <  2:263762/150‑1 (MQ=34)
                                           ctctGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAAtttt                                                                                                            >  1:17862/1‑150 (MQ=34)
                                           ctctGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAAGTGGGCTGCCCCGGCCTTCAGGGACCGCCTCGGCAAtttt                                                                                                            >  2:440024/1‑150 (MQ=21)
                                            tctGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTg                                                                                                           <  1:1889920/150‑1 (MQ=34)
                                            tctGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTg                                                                                                           <  2:1839252/150‑1 (MQ=34)
                                            tctGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGGCGGGCGCGAGGTAACATCCCTGTATAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTg                                                                                                           <  1:2194529/150‑1 (MQ=21)
                                             ctGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACCTCCCTGTAAAACGCGCCCTCAGCCCCCATCCATGTGGGCTGCCCCGGCCTTCAGGGGACGCCTCGGCAATTTTgg                                                                                                          >  1:2087432/1‑149 (MQ=255)
                                             ctGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCCGCCGGCATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTgg                                                                                                          >  1:1467556/1‑149 (MQ=255)
                                              tGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGAc                                                                                                         <  1:998021/150‑1 (MQ=34)
                                               ggCGCACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGACCGCCTCGGCAATTTTGCCg                                                                                                        >  1:1207984/1‑150 (MQ=255)
                                               ggCGAACAAACTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACg                                                                                                        <  1:512397/150‑1 (MQ=21)
                                               ggCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGTGAACGCCTCGGCAATTTTGACg                                                                                                        <  1:2055815/150‑1 (MQ=21)
                                               ggCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGGACGCCTCGGCAATTTTGACg                                                                                                        >  1:542938/1‑150 (MQ=21)
                                               ggCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACg                                                                                                        <  1:1735476/150‑1 (MQ=34)
                                                 cGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCATACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACCCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGcc                                                                                                      >  2:2208598/1‑150 (MQ=21)
                                                  gAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCCCATACATGTGGGCTGCCCCCGCCATCAGGGCACCCCTCGGGCATTTTGACGccc                                                                                                     >  2:1840706/1‑149 (MQ=11)
                                                  gAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGGAACGTCCCTGTAAGACGCGCCCTCCGCCCCCAGACATGGTGGCTGCCCCGGCCGGCAGGGAACGCCTGGGCCATTTTGGCggcc                                                                                                     >  1:1080174/1‑147 (MQ=11)
                                                   aaCAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGccac                                                                                                    >  2:1149711/1‑150 (MQ=34)
                                                     cAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAACCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTGCAGGGAACGCCGCGGCAATTTTGACGccaccc                                                                                                  >  2:378471/1‑149 (MQ=255)
                                                     cAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGccacca                                                                                                  >  2:1245542/1‑150 (MQ=34)
                                                        aaTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAaac                                                                                               >  2:475970/1‑150 (MQ=35)
                                                           cgccgcTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGCGGGAACCTCCCTGTAAACCGCGCCCTCAGCCCACCTCCATGTGGGCTGCCCCGGCCGTCAGGGCCCGCCGCGGCCATTTTGCCGCCCCCCcacaac                                                                                            >  1:1881993/1‑150 (MQ=11)
                                                            gccgcTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAAcc                                                                                           >  2:162735/1‑150 (MQ=39)
                                                             ccgcTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTACCATCCCTGTAAAACGCGCCCTCAGCCCCCCTCCATGTGGGCTGCCCCGGCCTTCAGGGACCGCCTCGGCAAGTTTGACGCCGCCCATCACCCg                                                                                          >  1:99370/1‑150 (MQ=11)
                                                             ccgcTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTCAAACGCGCCCTCAGCCCACAGCCGTGTGGGCTGCCCCGGCCTTCAGGGCACGCCGCGGCAATTTTGACGCCACCAACCAACCg                                                                                          >  1:99373/1‑150 (MQ=11)
                                                             ccgcTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCg                                                                                          <  1:1250588/150‑1 (MQ=40)
                                                               gcTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTg                                                                                        <  2:726239/150‑1 (MQ=255)
                                                                  gcgcGGTGCCCTCGGCTTATCCCGTAAGGCTACCCGGTCGGGCGCGAGGTAACCTCCCTGTAAAACGCGCCCTCAGCCCACAACCCAGTGGGCTGCCCCCGCCGTGAGGGGACGCCACGGCAATTTTGACGCCACCCAACCACCGGGCgg                                                                                     >  2:739142/1‑150 (MQ=11)
                                                                   cgcgGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGc                                                                                    <  1:1173627/150‑1 (MQ=255)
                                                                   cgcgGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGc                                                                                    <  2:1792664/150‑1 (MQ=255)
                                                                       gTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTAt                                                                                <  1:1744120/150‑1 (MQ=255)
                                                                         gCCCTCGGCTTATCCCTTACGGCTACCGGGCCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTg                                                                              <  1:366804/150‑1 (MQ=255)
                                                                          cccTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAAATTTGACGCCCCCAAACAACCGTGCGGCCTAATGa                                                                             >  2:1957760/1‑150 (MQ=255)
                                                                           ccTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGAt                                                                            <  2:653106/150‑1 (MQ=255)
                                                                            cTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTTAACCGCGCCCTCAGCCCCCATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCGCGGCAATTTTGACGCCACCAAACAACCGGGCGGCCGATTGATa                                                                           >  1:921009/1‑150 (MQ=14)
                                                                             tCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATaa                                                                          <  1:2128843/150‑1 (MQ=255)
                                                                                gCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAaga                                                                       <  2:49089/150‑1 (MQ=255)
                                                                                  ttATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGc                                                                     <  1:1386987/150‑1 (MQ=255)
                                                                                    aTCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACATCCGTGCGGCCTATTGATAAAGAGCTa                                                                   <  1:1983091/150‑1 (MQ=255)
                                                                                       ccTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAaca                                                                <  2:817899/150‑1 (MQ=255)
                                                                                        cTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGACCGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTTacac                                                               >  1:88489/1‑150 (MQ=255)
                                                                                          tACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGGACGCCTCGGCAATTTTGACGCCCCCACACAACCGTGCGGCCTATTGATAATGAGCTGACACAt                                                             >  1:175782/1‑150 (MQ=17)
                                                                                              gCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCTATTTTGACGCCACCAAACACCCGTGCGGCCTATTGATAAAGAGCTAACACATTGTc                                                         >  1:1111482/1‑150 (MQ=255)
                                                                                              gCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTc                                                         <  2:1647874/150‑1 (MQ=255)
                                                                                               cTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAGCCGTGCGGCGTAGTGAGAAAGAGCTAACACATGGTCa                                                        >  1:169097/1‑150 (MQ=17)
                                                                                                 aCCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCaaa                                                      >  2:797876/1‑150 (MQ=255)
                                                                                                 aCCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCaaa                                                      >  2:233434/1‑150 (MQ=255)
                                                                                                    gggTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCaaaaaa                                                   >  2:1130757/1‑150 (MQ=255)
                                                                                                        cGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATc                                               <  1:397503/150‑1 (MQ=255)
                                                                                                          ggCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCAc                                             >  2:639004/1‑150 (MQ=255)
                                                                                                              cgAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGGCGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAACCATCACTATg                                         >  1:1418080/1‑150 (MQ=255)
                                                                                                              cgAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATg                                         >  1:1681406/1‑150 (MQ=255)
                                                                                                              cgAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTCCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATg                                         >  1:829242/1‑150 (MQ=255)
                                                                                                                 ggTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtt                                      >  1:1603727/1‑150 (MQ=255)
                                                                                                                  gTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGttt                                     >  2:1264673/1‑150 (MQ=255)
                                                                                                                   tAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCAATATGGtttt                                    >  2:305474/1‑150 (MQ=255)
                                                                                                                    aaCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGGGCGGCCGATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGttttt                                   >  2:143993/1‑150 (MQ=255)
                                                                                                                     aCATCCCTTTAAAACGCGCCCTTAGACCACATCCATGTCGGCTGCCCCGGCCTTCAGGGAACGCCTCGTCAATTTTCACGCCACCAAACAACCGTGCTGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtttttt                                  <  1:1710685/150‑1 (MQ=25)
                                                                                                                     aCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtttttt                                  <  1:1245542/150‑1 (MQ=255)
                                                                                                                       aTCCCTGTAAATCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCATTTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACTCATTGTCAAAAAACATCACTATGGTTTTTTag                                <  2:2208818/150‑1 (MQ=255)
                                                                                                                          ccTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGGGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGt                             >  1:1279192/1‑150 (MQ=255)
                                                                                                                           cTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGtt                            >  2:1898318/1‑150 (MQ=255)
                                                                                                                                aaaCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGTGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCg                       <  1:1103557/150‑1 (MQ=255)
                                                                                                                                aaaCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGAAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCg                       <  1:969044/150‑1 (MQ=255)
                                                                                                                                   cgcgCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTTGTGATATAGAGGGAGCACATTGTCAAAAAACATCACGATGGTTTTTTAGAGTTTCTCGata                    >  1:1825582/1‑150 (MQ=25)
                                                                                                                                   cgcgCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGata                    <  2:1728715/150‑1 (MQ=255)
                                                                                                                                   cgcgCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGata                    >  2:1432730/1‑150 (MQ=255)
                                                                                                                                     cgcCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATc                  <  2:17862/150‑1 (MQ=255)
                                                                                                                                        ccTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAt               >  2:451004/1‑150 (MQ=255)
                                                                                                                                           cAGCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCGATTGATAAACAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGc            >  2:1519208/1‑150 (MQ=255)
                                                                                                                                             gCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCCCCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCt          >  1:1640778/1‑150 (MQ=255)
                                                                                                                                              cccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTCACACATTGTCAAAAAACATCACTATGGTGTTTTAGAGTTTCTCGATATAAATTGCCTg         >  1:372831/1‑150 (MQ=255)
                                                                                                                                               caacaTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGa        >  1:1579473/3‑150 (MQ=255)
                                                                                                                                                  caTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCGTATTGATAAAGAGCTAACACATTGTCAAAAAACATCAGTATGGTTTTTTAGAGGTTCTCGATATCAATTGCCTGAATa     >  2:364490/1‑150 (MQ=255)
                                                                                                                                                  caTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAg    >  2:1018214/1‑150 (MQ=255)
                                                                                                                                                     ccATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAGcc  >  2:1299519/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAGCC  >  NC_000913/2229173‑2229471

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