Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,996,736 G→C 5.1% Q33E (CAG→GAG) ‡ sdiA ← DNA‑binding transcriptional dual regulator SdiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,996,7360GC 5.1% 336.5 / 6.4 99Q33E (CAG→GAG) ‡sdiADNA‑binding transcriptional dual regulator SdiA
Reads supporting (aligned to +/- strand):  ref base G (50/44);  new base C (2/3);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 101 total aligned reads displayed.

CGGATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAATTCAAAC  >  NC_000913/1996590‑1996885
                                                                                                                                                  |                                                                                                                                                     
cGGATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAg                                                                                                                                                    <  1:337027/150‑1 (MQ=255)
 ggATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGc                                                                                                                                                   >  2:1360427/1‑150 (MQ=255)
  gATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCt                                                                                                                                                  >  1:248748/1‑150 (MQ=255)
      aaTTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAAt                                                                                                                                              <  1:1215145/150‑1 (MQ=255)
        ttGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCTTCAGGCTCagcag                                                                                                                                            <  1:1190726/150‑5 (MQ=255)
        ttGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCTTCAGGCTCagcag                                                                                                                                            <  1:1407992/150‑5 (MQ=255)
        ttGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGGGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAAttt                                                                                                                                            <  1:1108289/150‑1 (MQ=255)
          gCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCa                                                                                                                                          <  2:1352111/150‑1 (MQ=255)
                  gTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGAc                                                                                                                                  >  1:505795/1‑150 (MQ=255)
                  gTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGAc                                                                                                                                  >  2:727234/1‑150 (MQ=255)
                   ttttttGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGAcc                                                                                                                                 <  2:568786/150‑1 (MQ=255)
                        tGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCtcttc                                                                                                                            <  2:1348073/150‑1 (MQ=255)
                          ccTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTg                                                                                                                          >  2:1202048/1‑150 (MQ=255)
                          ccTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTCCTGAGCCTGAAGCTCAATTTCATGTTAGACCtcttcta                                                                                                                          >  2:1387092/1‑149 (MQ=255)
                           cTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTAATCCAGCTGCTGAGCATGAATCTCAATTTCATGGTAGAACTATTCTGc                                                                                                                         <  2:1079533/150‑1 (MQ=255)
                             gATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCgg                                                                                                                       <  1:1485688/150‑1 (MQ=255)
                              ataataACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGt                                                                                                                      >  1:1512294/1‑150 (MQ=255)
                               taataaCTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGtc                                                                                                                     >  1:472552/1‑150 (MQ=255)
                               taataaCTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGtc                                                                                                                     >  1:1502699/1‑150 (MQ=255)
                                aataaCTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGtct                                                                                                                    <  2:1169329/150‑1 (MQ=255)
                                 ataaCTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTGGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGtctc                                                                                                                   <  1:1101689/150‑1 (MQ=255)
                                 ataaCTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGtctc                                                                                                                   <  1:933923/150‑1 (MQ=255)
                                 ataaCTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGtctc                                                                                                                   <  2:1438443/150‑1 (MQ=255)
                                  taactaacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTcc                                                                                                                  <  2:1207205/150‑1 (MQ=255)
                                  taactaacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTcc                                                                                                                  <  2:454874/150‑1 (MQ=255)
                                  taactaacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTcc                                                                                                                  >  2:546276/1‑150 (MQ=255)
                                  taactaacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTcc                                                                                                                  >  2:122067/1‑150 (MQ=255)
                                  aaactaacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTCCGGTCTcc                                                                                                                  <  2:130993/149‑1 (MQ=255)
                                    actaacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCAt                                                                                                                <  2:497959/150‑1 (MQ=255)
                                     ctaacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATc                                                                                                               >  1:1099432/1‑150 (MQ=255)
                                       aacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCct                                                                                                             >  1:304620/1‑150 (MQ=255)
                                       aacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCct                                                                                                             >  1:823798/1‑150 (MQ=255)
                                       aacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCct                                                                                                             >  2:21303/1‑150 (MQ=255)
                                        acCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCctc                                                                                                            >  2:146449/1‑150 (MQ=255)
                                         cccACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGTAGCTCAATTTCATGGTAGACCTCTTCTGCTGTCTCCATCctct                                                                                                           <  2:662716/150‑1 (MQ=255)
                                         cccACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCTGTCTCCATCctct                                                                                                           >  1:143086/1‑150 (MQ=255)
                                         cccACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCctct                                                                                                           >  1:1036234/1‑150 (MQ=255)
                                         cccACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCctct                                                                                                           <  1:137109/150‑1 (MQ=255)
                                         cccACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCctct                                                                                                           >  1:26104/1‑150 (MQ=255)
                                           cACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGa                                                                                                         >  1:526000/1‑150 (MQ=255)
                                             cGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGaaa                                                                                                       <  1:721994/150‑1 (MQ=255)
                                               ccTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAAc                                                                                                     <  1:21303/150‑1 (MQ=255)
                                                cTCAGGGTAATTGGTGTAAAAAGCCACTTAAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACg                                                                                                    <  2:1093815/150‑1 (MQ=255)
                                                    gggTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAAc                                                                                                <  2:771823/150‑1 (MQ=255)
                                                    gggTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAAc                                                                                                <  2:325377/150‑1 (MQ=255)
                                                             gtgtAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTg                                                                                       >  2:1030749/1‑150 (MQ=255)
                                                               gtAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCg                                                                                     >  2:1341511/1‑150 (MQ=255)
                                                                 aaaaaGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGAc                                                                                   >  1:842915/1‑150 (MQ=255)
                                                                  aaaaGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACg                                                                                  <  1:950678/150‑1 (MQ=255)
                                                                  aaaaGCCAATTTAGGACGAGTGAATGGTACCGGGTGCAGGACACATAACGAATAGTAATCGTAATCAAGCTGCTGAGCCTGAAGCTAAATTTCATGGTAGACCTCTTCTGAGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACg                                                                                  >  2:492745/1‑150 (MQ=255)
                                                                    aaGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGcc                                                                                <  2:1440465/150‑1 (MQ=255)
                                                                      gCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAg                                                                              >  1:400241/1‑150 (MQ=255)
                                                                      gCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAg                                                                              >  1:1622545/1‑150 (MQ=255)
                                                                      gCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAg                                                                              >  1:1628004/1‑150 (MQ=255)
                                                                       ccAATTTATGTCGAGTGAATTGTACCGGGTGGCGGACAAATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGc                                                                             <  1:1387092/150‑1 (MQ=255)
                                                                        cACTTTAGGTCGAGTGAATGGTACCGGGTGTCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCt                                                                            <  1:951116/150‑1 (MQ=255)
                                                                         aCTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTg                                                                           >  2:829864/1‑150 (MQ=255)
                                                                         aCTTTAGGTCGAGTGAAGGGTACCGGGGGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTg                                                                           <  2:373361/150‑1 (MQ=255)
                                                                                gTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAAt                                                                    <  1:122067/150‑1 (MQ=255)
                                                                                  cGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATcc                                                                  <  2:505795/150‑1 (MQ=255)
                                                                                   taGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCt                                                                 <  1:1088582/149‑1 (MQ=255)
                                                                                   gAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCt                                                                 <  2:1512294/150‑1 (MQ=255)
                                                                                   gAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCt                                                                 <  2:923274/150‑1 (MQ=255)
                                                                                             tACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCAt                                                       <  1:230781/150‑1 (MQ=255)
                                                                                              aCCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATa                                                      <  2:1051997/150‑1 (MQ=255)
                                                                                              aCCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATa                                                      >  1:1048090/1‑150 (MQ=255)
                                                                                                 gggTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTa                                                   <  1:1341511/150‑1 (MQ=255)
                                                                                                  ggTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTaa                                                  >  1:86727/1‑150 (MQ=255)
                                                                                                   gTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTaaa                                                 >  1:1223808/1‑150 (MQ=255)
                                                                                                       cGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGc                                             >  2:630727/1‑150 (MQ=255)
                                                                                                           cacaTAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACg                                         >  1:812479/1‑150 (MQ=255)
                                                                                                             caTAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGcc                                       >  1:1629214/1‑150 (MQ=255)
                                                                                                             caTAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGcc                                       >  2:1481502/1‑150 (MQ=255)
                                                                                                             caTAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGcc                                       >  2:447017/1‑150 (MQ=255)
                                                                                                               tAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGcccc                                     <  2:200248/150‑1 (MQ=255)
                                                                                                                aaCGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCt                                    >  2:1424625/1‑150 (MQ=255)
                                                                                                                     aTAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGagag                               >  1:935244/1‑150 (MQ=255)
                                                                                                                     aTAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGagag                               >  2:408537/1‑150 (MQ=255)
                                                                                                                         tAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAgagtgag                           <  1:1476678/150‑1 (MQ=255)
                                                                                                                            tCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTAt                        <  2:400241/150‑1 (MQ=255)
                                                                                                                               tACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCAt                     <  2:1048090/150‑1 (MQ=255)
                                                                                                                                 ctacAGCTCCTGATCATGAACATCAATTTCATGGTATACATCTTATAAAATATACATCCTCTGAAAACGAAACAGCATCGTGCGACGCCAGATGAAAAAATCCTTATACTGAATAGTAAACCGCAACGCCACTGAGAGTGAGTATCAatt                   >  2:751700/4‑150 (MQ=255)
                                                                                                                                 cTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCAttt                   <  1:1360427/150‑1 (MQ=255)
                                                                                                                                   ccAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGTCATCTGAAATATTCGTTATCCTGCATAGTAAACCGCAACTCACCTGAGAGTGAGTATCATTtat                 <  2:817508/150‑1 (MQ=255)
                                                                                                                                   ccAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTtat                 >  1:1215570/1‑150 (MQ=255)
                                                                                                                                   ccAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTtat                 >  2:1514194/1‑150 (MQ=255)
                                                                                                                                   ccAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCATTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTtat                 >  2:50126/1‑150 (MQ=255)
                                                                                                                                    atgaaattgAGCTTCAGGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTtata                >  2:121126/8‑150 (MQ=255)
                                                                                                                                        aattgAGCTTCAGGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATaa               <  2:805108/144‑1 (MQ=255)
                                                                                                                                        aattgAGCTTCAGGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATaa               >  1:805108/4‑147 (MQ=255)
                                                                                                                                         gctgAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGAt           >  1:167309/1‑150 (MQ=255)
                                                                                                                                         gctgAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGAt           >  1:448178/1‑150 (MQ=255)
                                                                                                                                            gAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAAt        >  1:766759/1‑150 (MQ=255)
                                                                                                                                             aGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAAtt       <  2:302423/150‑1 (MQ=255)
                                                                                                                                               ccTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGAAGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTAGAATGATAATTCa     >  2:639613/1‑150 (MQ=255)
                                                                                                                                                cTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCACCGCCCCTGAGAGTGAGTATCATTTATAATGATAATTCaa    >  1:375776/1‑150 (MQ=255)
                                                                                                                                                cTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAATTCaa    <  1:492745/150‑1 (MQ=255)
                                                                                                                                                 tGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAATTCaaa   >  1:51379/1‑150 (MQ=255)
                                                                                                                                                  gAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAATTCAAAc  <  2:526000/150‑1 (MQ=255)
                                                                                                                                                  gAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAATTCAAAc  <  2:448178/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAATTCAAAC  >  NC_000913/1996590‑1996885

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