Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,636,300 T→G pseudogene (68/528 nt) nohQ ← pseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,636,3000TG99.1% 316.7 / NA 106pseudogene (68/528 nt)nohQpseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (65/40);  total (65/41)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 106 total aligned reads displayed.

CTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTC  >  NC_000913/1636155‑1636447
                                                                                                                                                 |                                                                                                                                                   
cTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGcc                                                                                                                                                  >  2:1543637/1‑149 (MQ=38)
 tGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGt                                                                                                                                               >  1:501523/1‑151 (MQ=38)
 tGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGt                                                                                                                                               >  1:1333853/1‑151 (MQ=38)
 tGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGt                                                                                                                                               >  1:1098480/1‑151 (MQ=38)
  gCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGtt                                                                                                                                              >  2:1177167/1‑151 (MQ=38)
      cAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTg                                                                                                                                           >  2:508112/1‑150 (MQ=25)
          tcttcAACTGCCCGGCGCAGCTTTTCGTTCTCAAGTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTGCTGAATg                                                                                                                                       <  2:661216/150‑1 (MQ=255)
            ttcAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATgg                                                                                                                                      >  2:444334/1‑149 (MQ=25)
             tcAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATgg                                                                                                                                      <  1:1116720/148‑1 (MQ=25)
             tcAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGt                                                                                                                                     <  1:558439/149‑1 (MQ=25)
             tcAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTAc                                                                                                                                   >  1:1374622/1‑151 (MQ=25)
              cAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTAc                                                                                                                                   <  1:964441/150‑1 (MQ=25)
              cAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTAc                                                                                                                                   <  1:1574864/150‑1 (MQ=25)
              cAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACg                                                                                                                                  <  2:351045/151‑1 (MQ=25)
                 cTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACTGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGAt                                                                                                                               >  1:674501/1‑151 (MQ=25)
                      cGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCg                                                                                                                          >  2:1629132/1‑151 (MQ=18)
                      cGGAGCATCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCGTACCATCTTATAACGGCGGTAGAGTCATCAAGCACCTCATTACCCTTTCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGTCAGTTCTGAATGGTACGGATCCTc                                                                                                                           >  2:1176550/1‑150 (MQ=11)
                         cgcAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCAc                                                                                                                       >  1:329933/1‑151 (MQ=18)
                           cAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATGCCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCg                                                                                                                     <  1:627108/151‑1 (MQ=18)
                           cAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCg                                                                                                                     >  2:88289/1‑151 (MQ=18)
                            aGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCg                                                                                                                     >  1:1009455/1‑150 (MQ=18)
                             gCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCTCCTCATTGCCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGaa                                                                                                                   <  2:1445991/151‑1 (MQ=16)
                                tttCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGaaaa                                                                                                                 <  1:954182/150‑1 (MQ=18)
                                 ttCGTTCTCAATTTCAGCATCCCTTTCGTCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAAt                                                                                                                <  2:1582527/150‑1 (MQ=18)
                                   cGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTc                                                                                                             >  1:758291/1‑151 (MQ=18)
                                    gTTCTCAATTTCAGCATCCCTGTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCGCCTCATTGCCCTTGCCGCCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTc                                                                                                             <  1:1596541/150‑1 (MQ=255)
                                         cAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTcagcca                                                                                                        <  2:645888/150‑1 (MQ=18)
                                         cAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTcagc                                                                                                          <  1:1617897/148‑1 (MQ=18)
                                           aTTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAgct                                                                                                     >  1:1479253/1‑151 (MQ=18)
                                           aTTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAgct                                                                                                     >  1:1292146/1‑151 (MQ=18)
                                             ttCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTATCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAgctgc                                                                                                   >  2:1142581/1‑151 (MQ=18)
                                               cAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCtt                                                                                                 <  1:687623/151‑1 (MQ=18)
                                               cAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCtt                                                                                                 <  1:1162583/151‑1 (MQ=18)
                                                  cATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCttttt                                                                                              <  2:1393337/151‑1 (MQ=18)
                                                  cATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCttttt                                                                                              >  1:342493/1‑151 (MQ=18)
                                                     cccTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCTGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTttgtt                                                                                           >  1:1007389/1‑151 (MQ=18)
                                                      ccTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTttgttg                                                                                          <  2:947249/151‑1 (MQ=18)
                                                         ttCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACt                                                                                       >  1:432539/1‑151 (MQ=18)
                                                          tCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACtt                                                                                      >  1:1570703/1‑151 (MQ=18)
                                                             gCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTcc                                                                                    >  1:270543/1‑150 (MQ=18)
                                                             gCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCa                                                                                   <  1:51009/151‑1 (MQ=18)
                                                              cATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCAt                                                                                  <  1:387762/151‑1 (MQ=18)
                                                               aTACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCCCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCAtt                                                                                 >  2:536974/1‑151 (MQ=18)
                                                                  ccATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGtt                                                                              <  1:483647/151‑1 (MQ=25)
                                                                    aTCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTc                                                                             >  2:1036563/1‑150 (MQ=25)
                                                                      cTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCAtt                                                                          >  2:1149786/1‑151 (MQ=37)
                                                                        tataACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTcc                                                                        >  1:1487940/1‑151 (MQ=25)
                                                                        tataACGGCGGCAGAATCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTcc                                                                        >  2:1020182/1‑151 (MQ=25)
                                                                         ataACGGCGGCAGAGTCATAAAGCACCTCATTACCCTCGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGGACGGATACTCGCACCGAAAATGGCAGCCAGCGGCTTTTTGTTGACTTCCATTGTTCATTcc                                                                        >  2:572211/1‑150 (MQ=255)
                                                                          taACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCAc                                                                      >  1:1277796/1‑151 (MQ=25)
                                                                          taACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCAc                                                                      >  2:139128/1‑151 (MQ=25)
                                                                          taACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCAc                                                                      >  1:1126020/1‑151 (MQ=25)
                                                                           aaCGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACg                                                                     >  1:420171/1‑151 (MQ=21)
                                                                           aaCGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACg                                                                     >  1:973711/1‑151 (MQ=21)
                                                                                 ggcAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCaaa                                                               <  2:507294/151‑1 (MQ=25)
                                                                                    agagTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAAc                                                            >  1:1394554/1‑151 (MQ=25)
                                                                                    agagTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAAc                                                            >  1:246049/1‑151 (MQ=25)
                                                                                    agagTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAAc                                                            >  1:264117/1‑151 (MQ=25)
                                                                                    agagTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAAc                                                            >  1:284445/1‑151 (MQ=25)
                                                                                     gagTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACa                                                           <  1:965935/151‑1 (MQ=25)
                                                                                      agTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAAc                                                            >  1:637480/1‑149 (MQ=25)
                                                                                        tCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACag                                                          >  2:1638170/1‑149 (MQ=25)
                                                                                           tAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGaa                                                      <  2:559724/150‑1 (MQ=25)
                                                                                            aaaGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACagag                                                        >  2:905630/1‑147 (MQ=25)
                                                                                            aaaGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGaaa                                                     <  2:693246/150‑1 (MQ=25)
                                                                                            aaaGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGaaa                                                     <  2:1473309/150‑1 (MQ=25)
                                                                                              aGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGa                                                  <  1:871253/151‑1 (MQ=25)
                                                                                               gCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGaa                                                 <  2:1528294/151‑1 (MQ=25)
                                                                                                   cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAAc                                               >  2:1028900/1‑149 (MQ=25)
                                                                                                   cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACg                                              >  1:1078778/1‑150 (MQ=25)
                                                                                                     cATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGa                                             <  1:1094921/149‑1 (MQ=25)
                                                                                                     cATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACa                                           <  2:1291202/151‑1 (MQ=25)
                                                                                                       ttACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACag                                          <  2:1623609/150‑1 (MQ=25)
                                                                                                         aCCCTTCACTCCGCCTCGCAGGACGGACACTCCCTTAGCCTGCCAGCTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTTACTTCCATTGTGCATGCCACGGCCAAAAACAGAGAAAGGAAAcagga                                           <  2:196074/147‑5 (MQ=11)
                                                                                                           ccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGcc                                     >  1:191349/1‑151 (MQ=25)
                                                                                                           ccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGcc                                     >  1:1630467/1‑151 (MQ=25)
                                                                                                            cTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGccc                                    >  2:1100953/1‑151 (MQ=37)
                                                                                                             ttGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGccc                                    >  1:95888/1‑150 (MQ=37)
                                                                                                               gCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCaa                                  <  1:982353/150‑1 (MQ=38)
                                                                                                                ccACCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCaaa                                 >  2:1313388/1‑150 (MQ=38)
                                                                                                                  aCCGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGc                              >  1:1635238/1‑151 (MQ=38)
                                                                                                                    cGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTc                            >  1:1455377/1‑151 (MQ=38)
                                                                                                                    cGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTc                            <  1:1218833/151‑1 (MQ=38)
                                                                                                                    cGCCTCGCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTc                            >  1:766479/1‑151 (MQ=38)
                                                                                                                          gCAGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTc                      <  2:1211276/151‑1 (MQ=255)
                                                                                                                            aGAACGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAg                    <  1:1121630/151‑1 (MQ=255)
                                                                                                                               aCGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCAc                 >  2:1494391/1‑151 (MQ=255)
                                                                                                                                cGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCAcc                >  2:1631255/1‑151 (MQ=255)
                                                                                                                                cGGGCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCAcc                >  2:511445/1‑151 (MQ=255)
                                                                                                                                  ggCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTg              <  1:961883/151‑1 (MQ=255)
                                                                                                                                   gCATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGt             >  1:1631478/1‑151 (MQ=255)
                                                                                                                                    cATTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTc            >  1:736089/1‑151 (MQ=255)
                                                                                                                                     aTTCCCTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCg           >  1:1516141/1‑151 (MQ=255)
                                                                                                                                         ccTGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTc       >  1:1034751/1‑151 (MQ=255)
                                                                                                                                           tGTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGCAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCct     >  2:758532/1‑151 (MQ=255)
                                                                                                                                           tGTTCCGGACAGTTTTGAATGGTACGGATACTCGCCCCGAAAATGTCTGCCAACGGCTTTTTGTTGACTTCCATTGTTCATTCAACGGACAAAAACAGAGGAAGAAAACTACAGATGGCCAAAACCTCGGCTTCAGCAACTGTCGTTTcc      >  2:713881/1‑150 (MQ=11)
                                                                                                                                            gTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCct     <  2:202454/150‑1 (MQ=255)
                                                                                                                                            gTTCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCct     <  1:371458/150‑1 (MQ=255)
                                                                                                                                              tCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttc  >  2:180011/1‑151 (MQ=255)
                                                                                                                                              tCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCcttt   <  1:772923/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTC  >  NC_000913/1636155‑1636447

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