Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,512 G→A 6.6% M1M (GTG→ATG) †‡ nrfE → heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,291,5120GA 6.6% 379.5 / 10.2 122M1M (GTG→ATG) †‡nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (58/55);  new base A (6/2);  total (65/57)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 123 total aligned reads displayed.

GTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTC  >  NC_000913/4291368‑4291660
                                                                                                                                                 |                                                                                                                                                    
gTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCCTTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                  <  1:107805/150‑1 (MQ=255)
gTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTg                                                                                                                                                 <  2:558672/150‑1 (MQ=255)
   tGGTCGTCGTGTTGATGATGCGCTTTTTCATTCTCTACGACGGACAGTTGACGGTGGCGTAAGCCAGAATAGAGGTGGTTTCTGGACGTATTCCTTCCTTAAGTCGG‑GTTTCTGGCGGTTTTGTTAAGTCTCGGGGTCAACGTGTTGAcc                                                                                                                                              <  2:834224/150‑1 (MQ=255)
   tGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCATAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAcc                                                                                                                                              <  2:1193658/150‑1 (MQ=255)
     gTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAcccc                                                                                                                                            <  2:1662930/150‑1 (MQ=255)
      tCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCg                                                                                                                                           >  1:1333874/1‑150 (MQ=255)
       cGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGt                                                                                                                                          >  2:1541577/1‑150 (MQ=255)
       cGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                          <  2:1688485/150‑2 (MQ=255)
         gCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGGGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTg                                                                                                                                        <  2:1520133/150‑1 (MQ=255)
            tgtTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGGTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACg                                                                                                                                     <  1:16599/150‑1 (MQ=255)
             gtTGATGCGGCGCTTTCTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTGGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACgg                                                                                                                                    <  2:962436/150‑1 (MQ=255)
              ttGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGGCAACGTGTTGACCCCGTTGACGGc                                                                                                                                   <  2:180477/150‑1 (MQ=255)
               tGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGcc                                                                                                                                  >  2:1134857/1‑150 (MQ=255)
                   gctgcGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTcg                                                                                                                              <  2:43232/150‑1 (MQ=255)
                      gcgcTTTTTCCTTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCCGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                           <  1:488582/150‑1 (MQ=255)
                      gcgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                           >  2:588346/1‑150 (MQ=255)
                      gcgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCgg                                                                                                                           >  1:1576079/1‑150 (MQ=255)
                       cgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                                                                          <  2:92207/150‑1 (MQ=255)
                       cgcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGGTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                                                                          <  1:1126640/150‑1 (MQ=255)
                            tttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                                                                     >  2:903882/1‑150 (MQ=255)
                            tttCATTCTCTACGCCGGACAGTTGACGGGGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGc                                                                                                                     <  1:104119/150‑1 (MQ=255)
                               cATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                                                                                                                  >  2:144320/1‑150 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTTACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                              >  1:1401685/1‑150 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                              >  1:828222/1‑150 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcc                                                                                                              >  2:539643/1‑149 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGGCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                              <  2:325452/150‑1 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TGTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGCCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                              >  1:319463/1‑150 (MQ=255)
                                   ctctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCATTCCTGAAGACGT‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGAAGTGCGGTTgcg                                                                                                              >  1:729386/1‑150 (MQ=255)
                                    tctACGCCGGACAGTTGACGGTGGCGTAAGCCGGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                             >  1:1064387/1‑150 (MQ=255)
                                    tctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                             <  1:509751/150‑1 (MQ=255)
                                     ctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCt                                                                                                            >  2:429022/1‑150 (MQ=255)
                                       aCGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTgg                                                                                                          <  2:391494/150‑1 (MQ=255)
                                         gCCGGACAGGTGACGGGGGCGTAAGCCAGAAAAGAGGGGGGTGCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                                                                        <  1:116684/150‑1 (MQ=255)
                                         gCCGAACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                                                                        <  1:1691200/150‑1 (MQ=255)
                                                 gTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACTTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatga                                                                                                >  1:860275/1‑150 (MQ=255)
                                                   tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                                                                              >  1:996246/1‑150 (MQ=255)
                                                   tGACGGGGTCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                                                                              <  2:297815/150‑1 (MQ=255)
                                                    gccGGGGGCGTAAGCCAGAAAAGAGGTGGGTTCTGGACGGATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGc                                                                                             <  2:834922/148‑1 (MQ=255)
                                                      cGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCgc                                                                                           >  2:946378/1‑149 (MQ=255)
                                                       ggTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                                                                          <  2:241146/150‑1 (MQ=255)
                                                         tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTc                                                                                        >  2:31196/1‑150 (MQ=255)
                                                         tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTc                                                                                        <  1:1381506/150‑1 (MQ=255)
                                                         tGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTCGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTc                                                                                        >  1:146185/1‑150 (MQ=255)
                                                                 gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                <  1:801075/150‑1 (MQ=255)
                                                                 gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                <  1:360237/150‑1 (MQ=255)
                                                                 gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                >  2:964811/1‑150 (MQ=255)
                                                                 gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCAt                                                                                <  1:223051/150‑1 (MQ=255)
                                                                  ccAGAAAAGAGGGGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTCTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATc                                                                               <  1:75825/150‑1 (MQ=255)
                                                                   cAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCg                                                                              >  2:206935/1‑150 (MQ=255)
                                                                    aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGCTGCGACTCACTTGCATCgg                                                                             >  1:507996/1‑150 (MQ=255)
                                                                    aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                                                                             >  2:524261/1‑150 (MQ=255)
                                                                    aGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                                                                             <  1:964811/150‑1 (MQ=255)
                                                                     gAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGc                                                                            >  1:997793/1‑150 (MQ=255)
                                                                            ggtggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                                                                     >  2:212938/1‑150 (MQ=255)
                                                                            ggtggtTTATGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGAGTTGTTGTTAAGTCTCGGGGTCAACGTGTTTACCCCGATTACCGACTTCGCGGGAGAGCGGTTGCGCTGCCCTGCCATGATGCGACTCACTTGCATCGGCATTCTgg                                                                     >  1:1280589/1‑150 (MQ=255)
                                                                              tggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                                                                   <  2:1576079/150‑1 (MQ=255)
                                                                              gggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcg                                                                   <  2:884011/149‑1 (MQ=255)
                                                                               ggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCGACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGAAGCGACTCACTTGCATCGGCATTCTGgcgc                                                                  >  1:934702/1‑150 (MQ=255)
                                                                               ggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                                                                  >  2:1005665/1‑150 (MQ=255)
                                                                               ggtTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                                                                  <  2:539210/150‑1 (MQ=255)
                                                                                    cTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc                                                             >  1:751654/1‑150 (MQ=255)
                                                                                      ggACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc                                                           <  1:213601/150‑1 (MQ=255)
                                                                                            aTTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct                                                     >  1:1250470/1‑150 (MQ=255)
                                                                                            aTTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct                                                     >  1:857340/1‑150 (MQ=255)
                                                                                            aTTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAACCCGTTGACGGCCTTCGCGAGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct                                                     >  1:1028097/1‑150 (MQ=255)
                                                                                               ccttcctGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgct                                                  >  2:1026669/1‑150 (MQ=255)
                                                                                                cttcctGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgctg                                                 >  2:1203107/1‑150 (MQ=255)
                                                                                                   cctGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGTGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTc                                              <  2:1514385/150‑1 (MQ=255)
                                                                                                   cctGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTc                                              >  1:563301/1‑150 (MQ=255)
                                                                                                    ctGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                                                                        <  2:1499318/91‑1 (MQ=255)
                                                                                                    ctGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGcc                                                                                                        >  1:1499318/1‑91 (MQ=255)
                                                                                                      gAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGGGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGcc                                           <  1:1360736/150‑1 (MQ=255)
                                                                                                        aGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATAGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCtt                                         >  1:1310779/1‑150 (MQ=255)
                                                                                                         gTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCttt                                        >  1:983990/1‑150 (MQ=255)
                                                                                                            cgtttttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAGCCCGTTGACGGCCTTCGCGAGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTgg                                      <  2:1028097/147‑1 (MQ=255)
                                                                                                                  tCTGGCGTTGTTGTTTAGTCTCGGGGTCAACGGGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACt                                <  2:1064387/150‑1 (MQ=255)
                                                                                                                   cTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTg                               >  2:331786/1‑150 (MQ=255)
                                                                                                                           gttgttAAGTCTCGGGGTCTACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGGTGCGCTGGCCTGCCATGATGCGACTCTCTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTg                       <  1:239431/150‑1 (MQ=255)
                                                                                                                           gttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTg                       <  2:1250470/150‑1 (MQ=255)
                                                                                                                           gttgttAAGTCTCGGGGGCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTg                       <  2:374795/150‑1 (MQ=255)
                                                                                                                             tgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGtt                     >  2:1394535/1‑150 (MQ=255)
                                                                                                                              gttAAGTCTCGGGGTCAACGTGTTGACGCCGTTGACGGCCTTCGCGGGTGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt                    >  1:226963/1‑150 (MQ=255)
                                                                                                                              gttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttt                    >  1:1022907/1‑150 (MQ=255)
                                                                                                                                tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGttttg                  >  1:1138610/1‑149 (MQ=255)
                                                                                                                                tAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtc                  >  2:1216256/1‑150 (MQ=255)
                                                                                                                                  gtcaacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtctc                >  1:1285820/6‑150 (MQ=255)
                                                                                                                                   gTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGACTGCGGTAGAGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtca               <  2:1067664/150‑1 (MQ=255)
                                                                                                                                    caacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcat              >  1:753676/4‑150 (MQ=255)
                                                                                                                                     ctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatc             >  2:1488551/1‑150 (MQ=255)
                                                                                                                                     ctcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatc             <  1:1583569/150‑1 (MQ=255)
                                                                                                                                      acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatca            >  2:1004012/2‑150 (MQ=255)
                                                                                                                                      acGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGAGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatca            >  1:1179462/2‑150 (MQ=255)
                                                                                                                                       cGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAg           >  1:1256812/1‑150 (MQ=255)
                                                                                                                                        ggggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCGTTGGCGTACTGACGTATTGTTTTCTCATCAGc          >  1:1347154/1‑150 (MQ=255)
                                                                                                                                          ggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGa        >  1:188051/1‑150 (MQ=255)
                                                                                                                                           gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGAt       <  1:429022/150‑1 (MQ=255)
                                                                                                                                           gTCAACGTCTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCTCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTACGCGCTCCTGCTGCACGCCATTGGCCTACTGACGTATTGTTTTCTCATCAGCGAt       >  1:338549/1‑150 (MQ=255)
                                                                                                                                              aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCTTCAGCGGTTtc    >  2:25045/1‑150 (MQ=255)
                                                                                                                                               tggtgtTGACCCCGTTGACGGCCTTCGGGGGAGTGCGGGTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTTTTGTTTTCTCATCAGCGATTtct   <  1:1675288/148‑1 (MQ=255)
                                                                                                                                                cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc  <  1:1461090/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGG‑TTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTC  >  NC_000913/4291368‑4291660

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