Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,949,2080AC21.4% 400.9 / 127.6 207A254A (GCT→GCG) yifBmagnesium chelatase family protein and putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (88/74);  new base C (24/20);  total (112/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 209 total aligned reads displayed.

GGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGT  >  NC_000913/3949059‑3949354
                                                                                                                                                      |                                                                                                                                                   
ggTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCa                                                                                                                                                     >  1:1570409/1‑150 (MQ=255)
  tAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGc                                                                                                                                                   <  3:518247/150‑1 (MQ=255)
      tcatcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCAttt                                                                                                                                               >  4:87171/1‑150 (MQ=255)
      tcatcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCAttt                                                                                                                                               >  3:810330/1‑150 (MQ=255)
       catcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTa                                                                                                                                              >  4:346569/1‑150 (MQ=255)
        atcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTGTTTTTTGTACTGATTCAGCATTTa                                                                                                                                              >  2:64584/1‑150 (MQ=255)
        atcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTAc                                                                                                                                             >  2:1435178/1‑150 (MQ=255)
        atcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTAc                                                                                                                                             <  3:640683/150‑1 (MQ=255)
        atcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTAc                                                                                                                                             >  3:412644/1‑150 (MQ=255)
        atcaAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTAc                                                                                                                                             >  4:162789/1‑150 (MQ=255)
          caAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCa                                                                                                                                           <  3:403547/150‑1 (MQ=255)
           aaGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAg                                                                                                                                          >  2:2100383/1‑150 (MQ=255)
           aaGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAg                                                                                                                                          >  4:661599/1‑150 (MQ=255)
           aaGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAg                                                                                                                                          >  2:297670/1‑150 (MQ=255)
              aaaaaGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACt                                                                                                                                       >  1:744439/1‑150 (MQ=255)
                  aGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAat                                                                                                                                   >  2:1400793/1‑150 (MQ=255)
                  aGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAat                                                                                                                                   >  2:1835293/1‑150 (MQ=255)
                    cACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAatat                                                                                                                                 >  3:493665/1‑150 (MQ=255)
                       gCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCgc                                                                                                                              >  2:1651159/1‑150 (MQ=255)
                        ctgTTATGCGCCGGGCAAATTTCCCCTGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTTAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCgca                                                                                                                             <  2:2126352/148‑1 (MQ=255)
                        ccGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCgca                                                                                                                             <  2:1118458/150‑1 (MQ=255)
                        ccGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCgca                                                                                                                             <  1:870018/150‑1 (MQ=255)
                          gaaaTGTGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCgcagc                                                                                                                           <  2:1799900/147‑1 (MQ=255)
                                cgcCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACtctc                                                                                                                     <  2:1756222/150‑1 (MQ=255)
                                    cgCGAAATTTCCCCGGGCCCTGGAATTGCGCCCCCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCTTTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGAAGCACTCTCCAGt                                                                                                                 <  4:885379/149‑1 (MQ=255)
                                      cGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGc                                                                                                               >  1:1874736/1‑150 (MQ=255)
                                       gAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGcc                                                                                                              >  1:835442/1‑150 (MQ=255)
                                       gAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCACCAGTGcc                                                                                                              >  2:214520/1‑150 (MQ=255)
                                        aaaTTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCct                                                                                                             <  4:60796/150‑1 (MQ=255)
                                             tCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCa                                                                                                        >  1:103768/1‑150 (MQ=255)
                                                   ggCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCtt                                                                                                  <  2:1899593/150‑1 (MQ=255)
                                                   ggCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCtt                                                                                                  >  1:1331729/1‑150 (MQ=255)
                                                   ggCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCtt                                                                                                  <  2:1127661/150‑1 (MQ=255)
                                                     cccTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCGTTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTaa                                                                                                <  4:184761/150‑1 (MQ=255)
                                                      ccTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACCCTCCAGTGCCTCTTCATTGCTTaaa                                                                                               >  2:1703324/1‑150 (MQ=255)
                                                      ccTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTGCCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTaaa                                                                                               >  4:471020/1‑150 (MQ=255)
                                                      ccTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTaaa                                                                                               >  3:745672/1‑150 (MQ=255)
                                                      ccTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTaaa                                                                                               >  4:139609/1‑150 (MQ=255)
                                                       cTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAAt                                                                                              <  4:717912/150‑1 (MQ=255)
                                                       cTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGAGGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAAt                                                                                              >  2:430944/1‑150 (MQ=255)
                                                        tGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATc                                                                                             <  1:2197930/150‑1 (MQ=255)
                                                          gAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCACACTGAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTg                                                                                           <  1:641511/150‑1 (MQ=255)
                                                           aaTTGCGCCACTGCTTACCTTCGCGGTTAGCGGTGCACTGTGATGTGGTGGGCGGAACGGGCGCTGCCGCTATTG‑TTTTTGTTCTGATTCTGCTTTTTTTTGACTTAATATCGCGGCACCGTCCAGTGTCTCTTTATTGCTTAATTCTgg                                                                                          <  2:972615/150‑1 (MQ=255)
                                                            aTTGCGCCACCGCCTACCATCGCAGTTAACGCTGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGc                                                                                         >  1:1689490/1‑150 (MQ=255)
                                                                 gcCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAg                                                                                    <  2:706813/150‑1 (MQ=255)
                                                                   cACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGc                                                                                  >  2:2073312/1‑150 (MQ=255)
                                                                     ccGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACTGTCGCTGACGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCa                                                                                <  4:231724/150‑1 (MQ=255)
                                                                               tCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACAGc                                                                      <  2:1499506/150‑1 (MQ=255)
                                                                                  cAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTgg                                                                   <  1:1468388/150‑1 (MQ=255)
                                                                                   aGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGc                                                                  >  1:2110101/1‑150 (MQ=255)
                                                                                   aGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGc                                                                  >  2:2036132/1‑150 (MQ=255)
                                                                                     ttAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGa                                                                >  1:2181059/1‑150 (MQ=255)
                                                                                     ttAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGa                                                                <  4:262849/150‑1 (MQ=255)
                                                                                          gATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCAtt                                                           >  3:252056/1‑150 (MQ=255)
                                                                                           aTGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTg                                                          >  1:2216993/1‑150 (MQ=255)
                                                                                           aTGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTg                                                          >  2:1861255/1‑150 (MQ=255)
                                                                                             gCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACTGCTGGCGAGCATTGtt                                                        <  1:1805058/150‑1 (MQ=255)
                                                                                              cACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGttt                                                       <  2:675229/150‑1 (MQ=255)
                                                                                                cTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTa                                                     <  3:539143/150‑1 (MQ=255)
                                                                                                 tgtgATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTAc                                                    <  1:470848/150‑1 (MQ=255)
                                                                                                 tgtgATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTAc                                                    <  1:1870027/150‑1 (MQ=255)
                                                                                                          gTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCcg                                           >  1:266616/1‑150 (MQ=255)
                                                                                                          gTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCcg                                           >  1:1260820/1‑150 (MQ=255)
                                                                                                          gTGAGCGGAACGGGCGCTGACGCCATTG‑ATTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCGGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCcg                                           >  2:1920611/1‑150 (MQ=255)
                                                                                                             aGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCcggcg                                        <  1:1621941/150‑1 (MQ=255)
                                                                                                             aGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCcggcg                                        <  2:524789/150‑1 (MQ=255)
                                                                                                             aGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCcggcg                                        <  1:901253/150‑1 (MQ=255)
                                                                                                              gCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCcggcgg                                       >  2:1291639/1‑150 (MQ=255)
                                                                                                               cGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCggcggc                                      <  1:241837/150‑1 (MQ=255)
                                                                                                                ggAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCTATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGcc                                     <  4:821575/150‑1 (MQ=255)
                                                                                                                    cGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAAt                                <  1:681185/150‑1 (MQ=255)
                                                                                                                    cGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAAt                                >  1:924587/1‑150 (MQ=255)
                                                                                                                      ggagcTGCCGCCATTG‑TTTTTGTACTGTTTCAGCATTTACCAGACTTAATATCGCAGCACGCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCa                              <  2:1910991/147‑1 (MQ=255)
                                                                                                                      ggCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCa                              <  1:1324343/150‑1 (MQ=255)
                                                                                                                       gcgcTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAg                             >  3:851814/1‑150 (MQ=255)
                                                                                                                           tGCCGTCATTG‑TTTTTGTACTGATTCAGCATTTACTAGACTTAATATATTTGCACTCTTTAGTGCCTCTTCATAGTTTAAATCTTACAACAGGCCATTAATACGGCTGGTGAGAGTTGTTTTTGCTGTTTCCGGCTGCCTCAAT‑AGTaaa                         >  4:365003/1‑150 (MQ=255)
                                                                                                                           tGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTaaa                         >  2:114269/1‑150 (MQ=255)
                                                                                                                           tGCCGCCATAG‑TGTTTGTACTGAGTCAGGATTGACCAGACTTAACATAGCAGCAAACTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCCAGCCTTGTTTTACCTGTTCCCGGCGGCC‑CCATCactgaa                         >  4:749673/1‑145 (MQ=255)
                                                                                                                            gccgccATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTaaaa                        <  2:2214244/150‑1 (MQ=255)
                                                                                                                             ccgccATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAg                       >  1:966869/1‑150 (MQ=255)
                                                                                                                              cgccATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAgg                      >  1:977541/1‑150 (MQ=255)
                                                                                                                              cgccATTG‑TTTTTGTACTGATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAgg                      >  2:2001381/1‑150 (MQ=255)
                                                                                                                              cgccAGTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCCGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAgg                      <  1:314928/150‑1 (MQ=255)
                                                                                                                               gccATTG‑TTTTTTTATTTATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGt                     >  4:289194/1‑150 (MQ=255)
                                                                                                                               gccATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGt                     >  3:87417/1‑150 (MQ=255)
                                                                                                                                ccATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGtt                    <  4:689293/150‑1 (MQ=255)
                                                                                                                                ccATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGtt                    >  1:2273453/1‑150 (MQ=255)
                                                                                                                                  aTTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTtgt                  <  4:526557/150‑1 (MQ=255)
                                                                                                                                   ttG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTtgtg                 <  2:1007435/150‑1 (MQ=255)
                                                                                                                                       tttttGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGccc              <  3:905951/150‑1 (MQ=255)
                                                                                                                                         tttGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGc            <  1:1979507/150‑1 (MQ=255)
                                                                                                                                          ttGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGTCCATTAATACGGCTGGCGAGCATTGTTTTACCTGATCCCGGCGTCC‑CAATCAGTAAAAGGTAGTGCCCGcc           >  1:647943/1‑150 (MQ=255)
                                                                                                                                             tACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGCcagc        >  3:222570/1‑150 (MQ=255)
                                                                                                                                              aCTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGCcagca       >  3:160466/1‑150 (MQ=255)
                                                                                                                                               ccccctCAGCATTTACCAGACTTAATATCGCAGCTCTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGCcagcag      <  2:1133563/145‑1 (MQ=255)
                                                                                                                                                 gATTCCGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGCCAGCAGCg    <  1:1173299/150‑1 (MQ=255)
                                                                                                                                                 gATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGCCAGCAGCg    >  4:446963/1‑150 (MQ=255)
                                                                                                                                                  aTTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGTCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGCCAGCAGCgg   >  2:862314/1‑150 (MQ=255)
                                                                                                                                                  aTTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGCCAGCAGCgg   >  2:810041/1‑150 (MQ=255)
                                                                                                                                                   tccAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTGGTGCCCGCCAGCAGCGGt  <  1:2008599/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGGGCCCTGGAATTGCGCCACCGCCTACCATCGCAGTTAACGATGCACTGTGATGAGGTGAGCGGAACGGGCGCTGCCGCCATTG‑TTTTTGTACTGATTCAGCATTTACCAGACTTAATATCGCAGCACTCTCCAGTGCCTCTTCATTGCTTAAATCTGGCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCC‑CAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGT  >  NC_000913/3949059‑3949354

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