Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,996,734 C→T 5.1% Q33Q (CAG→CAA) ‡ 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,7340CT 5.1% 310.7 / 3.3 98Q33Q (CAG→CAA) ‡sdiADNA‑binding transcriptional dual regulator SdiA
Reads supporting (aligned to +/- strand):  ref base C (49/42);  new base T (2/3);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 102 total aligned reads displayed.

AGCACCGGATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAATTCAA  >  NC_000913/1996585‑1996883
                                                                                                                                                     |                                                                                                                                                     
aGCACCGGATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGcc                                                                                                                                                       <  2:1418097/150‑1 (MQ=255)
 gCACCGGATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCt                                                                                                                                                      >  2:1485688/1‑150 (MQ=255)
 gCACCGGATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCt                                                                                                                                                      >  1:1296328/1‑150 (MQ=255)
 gCACCGGATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCt                                                                                                                                                      >  2:687430/1‑150 (MQ=255)
     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:1407992/150‑5 (MQ=255)
             ttGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCTTCAGGCTCagcag                                                                                                                                          <  1:1190726/150‑5 (MQ=255)
             ttGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGGGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAAttt                                                                                                                                          <  1:1108289/150‑1 (MQ=255)
               gCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCa                                                                                                                                        <  2:1352111/150‑1 (MQ=255)
                       gTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGAc                                                                                                                                >  2:727234/1‑150 (MQ=255)
                       gTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGAc                                                                                                                                >  1:505795/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:546276/1‑150 (MQ=255)
                                       taactaacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTcc                                                                                                                >  2:122067/1‑150 (MQ=255)
                                       taactaacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTcc                                                                                                                <  2:454874/150‑1 (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                                                                                                           >  2:21303/1‑150 (MQ=255)
                                            aacCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCct                                                                                                           >  1:823798/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)
                                                                    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:1628004/1‑150 (MQ=255)
                                                                           gCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAg                                                                            >  1:1622545/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                                                    >  1:1048090/1‑150 (MQ=255)
                                                                                                   aCCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATa                                                    <  2:1051997/150‑1 (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                             >  2:408537/1‑150 (MQ=255)
                                                                                                                          aTAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGagag                             >  1:935244/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)
                                                                                                                                                     |                                                                                                                                                     
AGCACCGGATCAATTGCGAGAAAGTTTTTTGCCTGATAATAACTAACCCACGCCTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCTCTTCTGCGGTCTCCATCCTCTGAAAACGCAACAGCATCGTGCGACGCCAGCTGAAAAAATCCTTATCCTGCATAGTAAACCGCAACGCCCCTGAGAGTGAGTATCATTTATAATGATAATTCAA  >  NC_000913/1996585‑1996883

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: