Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,140,353 A→G 5.0% R149R (CGT→CGC hybE ← hydrogenase 2‑specific chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,140,3530AG 5.0% 433.6 / 7.4 140R149R (CGT→CGChybEhydrogenase 2‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base A (52/80);  new base G (3/4);  total (56/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

CGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACAC  >  NC_000913/3140204‑3140498
                                                                                                                                                     |                                                                                                                                                 
cGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTCCTCCGCCTGCCGTAGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTATCGCCACTCCTGCGACCAAACAGCACGGCGCGa                                                                                                                                                   >  2:1869191/1‑150 (MQ=255)
  cAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACg                                                                                                                                                 <  2:359695/150‑1 (MQ=255)
   aGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGc                                                                                                                                                <  1:1208643/150‑1 (MQ=255)
       gagCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTg                                                                                                                                            <  2:1912500/150‑1 (MQ=255)
         gcgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCg                                                                                                                                          <  2:979722/150‑1 (MQ=255)
         gcgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCg                                                                                                                                          <  2:1306180/150‑1 (MQ=255)
          cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGctg                                                                                                                                         <  2:1986882/150‑3 (MQ=255)
            cgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGctgtt                                                                                                                                       >  2:147990/1‑146 (MQ=255)
            cgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgtgt                                                                                                                                       <  2:214937/150‑1 (MQ=255)
             gcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCATAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGctgttt                                                                                                                                      <  1:1888945/150‑6 (MQ=255)
             gcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGtgtg                                                                                                                                      >  2:843351/1‑150 (MQ=255)
              ccAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTgg                                                                                                                                     >  2:698199/1‑150 (MQ=255)
               cAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGt                                                                                                                                    <  2:43623/150‑1 (MQ=255)
                aaTTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTa                                                                                                                                   <  2:2166140/150‑1 (MQ=255)
                     ccAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATcc                                                                                                                              >  1:706233/1‑150 (MQ=255)
                       aGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCgg                                                                                                                            <  1:1421645/150‑1 (MQ=255)
                        gCCACACGGCGGTGACGCTGTTAACATCGTGCTGCTCCGCCTTACGTTGGATAATTTCAACGGTGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGa                                                                                                                           <  2:2146748/150‑1 (MQ=255)
                        gCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGa                                                                                                                           >  1:991699/1‑150 (MQ=255)
                         ccACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGAt                                                                                                                          <  1:1338295/150‑1 (MQ=255)
                              cggcggTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTg                                                                                                                     <  2:2229283/150‑1 (MQ=255)
                              cggcggTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTg                                                                                                                     <  2:563968/150‑1 (MQ=255)
                                 cggTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACt                                                                                                                  <  2:8581/150‑1 (MQ=255)
                                    tGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGc                                                                                                               >  2:2349986/1‑150 (MQ=255)
                                     gACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCa                                                                                                              >  1:2010121/1‑150 (MQ=255)
                                      aCGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAg                                                                                                             <  1:1269530/150‑1 (MQ=255)
                                       cgcgCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGc                                                                                                            >  2:758979/1‑150 (MQ=255)
                                            ttAACATCGTGCTGCTCCGCCTACCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAg                                                                                                       <  2:392067/150‑1 (MQ=255)
                                             tAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAgg                                                                                                      >  1:2036807/1‑150 (MQ=255)
                                               aCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGAt                                                                                                    >  2:498387/1‑150 (MQ=255)
                                               aCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGAt                                                                                                    <  2:1805877/150‑1 (MQ=255)
                                               aCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGAt                                                                                                    <  1:994500/150‑1 (MQ=255)
                                               aCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGAt                                                                                                    >  1:1944292/1‑150 (MQ=255)
                                                cATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATc                                                                                                   <  1:2432598/150‑1 (MQ=255)
                                                   cGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCAtt                                                                                                >  1:1629352/1‑150 (MQ=255)
                                                       ctgctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTg                                                                                            <  1:1081764/150‑1 (MQ=255)
                                                       ctgctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTg                                                                                            <  1:2397860/150‑1 (MQ=255)
                                                         gctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgcg                                                                                          >  2:1758562/1‑150 (MQ=255)
                                                          ctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgcgc                                                                                         <  1:1324131/150‑1 (MQ=255)
                                                          ctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgcgc                                                                                         <  1:1605610/150‑1 (MQ=255)
                                                          ctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgcgc                                                                                         >  1:2306523/1‑150 (MQ=255)
                                                            ccGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAg                                                                                       <  1:274831/150‑1 (MQ=255)
                                                            ccGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGCATCATTCGTGCGCAg                                                                                       <  2:760291/150‑1 (MQ=255)
                                                              gCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTc                                                                                     <  1:830892/150‑1 (MQ=255)
                                                                cTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCAt                                                                                   >  1:2195646/1‑150 (MQ=255)
                                                                 tGCCGTTCGATAATTTCAACGGCGCTCTCGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGAGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATc                                                                                  >  1:1048720/1‑150 (MQ=255)
                                                                  gCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATcc                                                                                 >  2:631745/1‑150 (MQ=255)
                                                                    cGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGt                                                                               <  1:896729/150‑1 (MQ=255)
                                                                    cGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGt                                                                               <  2:1856026/150‑1 (MQ=255)
                                                                     gTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTc                                                                              <  1:1375732/150‑1 (MQ=255)
                                                                      ttGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCa                                                                             <  1:2283532/150‑1 (MQ=255)
                                                                       tGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAg                                                                            >  1:1201285/1‑150 (MQ=255)
                                                                       tGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAg                                                                            <  2:1308580/150‑1 (MQ=255)
                                                                        ggATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAgg                                                                           >  1:1657699/1‑150 (MQ=255)
                                                                          aTAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCg                                                                         <  2:191324/150‑1 (MQ=255)
                                                                                 cAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGccc                                                                  <  1:425698/150‑1 (MQ=255)
                                                                                 cAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGccc                                                                  <  1:2160991/150‑1 (MQ=255)
                                                                                  aaCGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCct                                                                 >  1:638931/1‑150 (MQ=255)
                                                                                  aaCGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCct                                                                 >  1:650769/1‑150 (MQ=255)
                                                                                         ctctGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAAt                                                          <  1:2349986/150‑1 (MQ=255)
                                                                                         ctctGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAAt                                                          <  2:1225155/150‑1 (MQ=255)
                                                                                          tctGGCAAAGCGACCACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCCGGGCGCGCCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCAGTCGTGCGCAGTCATCCGTCAGGCGTTGGCCAGCTTCaagc                                                         >  2:2482949/1‑148 (MQ=255)
                                                                                           ctGGCAAAGAGACAACTCATGCATTTTCACCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCg                                                        <  2:570147/150‑1 (MQ=255)
                                                                                            tGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGa                                                       >  2:1441058/1‑150 (MQ=255)
                                                                                             ggCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGAc                                                      >  2:1501188/1‑150 (MQ=255)
                                                                                               cAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACAt                                                    >  2:459228/1‑150 (MQ=255)
                                                                                                aaaGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATg                                                   >  2:513961/1‑150 (MQ=255)
                                                                                                aaaGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATg                                                   <  2:1211176/150‑1 (MQ=255)
                                                                                                   gagaCAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCtg                                                >  1:244730/1‑150 (MQ=255)
                                                                                                   cagaCAACTCATGCATTTTCGCCACTCCTGCGACAAAACAGCAGGGCGCGACGCCCTGGGTGTGGTACATCCGGATTCGTGACTGGCAGCTAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCtg                                                <  2:659835/149‑1 (MQ=255)
                                                                                                     gaCAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCtgtg                                              >  2:958849/1‑150 (MQ=255)
                                                                                                     gaCAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCtgtg                                              <  1:2380004/150‑1 (MQ=255)
                                                                                                       cAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCg                                            >  1:1090680/1‑150 (MQ=255)
                                                                                                           tCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAg                                        <  2:706233/150‑1 (MQ=255)
                                                                                                            cATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgc                                       >  1:909007/1‑150 (MQ=255)
                                                                                                             aTGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcg                                      <  1:2030038/150‑1 (MQ=255)
                                                                                                              tGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcgg                                     >  2:218663/1‑150 (MQ=255)
                                                                                                                cATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcggcg                                   <  1:1836390/150‑1 (MQ=255)
                                                                                                                 aTTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGa                                  <  2:1718326/150‑1 (MQ=255)
                                                                                                                  ttttCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGAc                                 <  1:876509/150‑1 (MQ=255)
                                                                                                                   tttCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACa                                <  2:1964543/150‑1 (MQ=255)
                                                                                                                     tCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATc                              <  2:1842718/150‑1 (MQ=255)
                                                                                                                      cGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCa                             >  2:1665117/1‑150 (MQ=255)
                                                                                                                       gCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAg                            >  1:2427997/1‑150 (MQ=255)
                                                                                                                         cACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCg                          <  1:914703/150‑1 (MQ=255)
                                                                                                                         cACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCg                          <  1:245979/150‑1 (MQ=255)
                                                                                                                            tCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAAc                       <  1:594055/150‑1 (MQ=255)
                                                                                                                            tCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAAc                       <  2:1629352/150‑1 (MQ=255)
                                                                                                                            tCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAAc                       <  2:1153086/150‑1 (MQ=255)
                                                                                                                               tGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAgg                    >  1:673407/1‑150 (MQ=255)
                                                                                                                                gCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGga                   <  2:244730/150‑1 (MQ=255)
                                                                                                                                   aCCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGgagag                >  1:364146/1‑150 (MQ=255)
                                                                                                                                    ccAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGgagaga               >  2:736820/1‑150 (MQ=255)
                                                                                                                                        cacgcagGGCGTCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATAtt           >  2:218556/4‑150 (MQ=255)
                                                                                                                                           gcagGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCg        <  1:591635/150‑1 (MQ=255)
                                                                                                                                           gcagGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCg        <  1:2012758/150‑1 (MQ=255)
                                                                                                                                             agGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGaa      <  1:551329/150‑1 (MQ=255)
                                                                                                                                              gggCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGaaa     >  2:628161/1‑150 (MQ=255)
                                                                                                                                              gggCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCATTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGaaa     >  2:1785835/1‑150 (MQ=255)
                                                                                                                                              gggCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCATTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGaaa     >  2:1776876/1‑150 (MQ=255)
                                                                                                                                                   cgACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAacac  <  1:731203/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAACAC  >  NC_000913/3140204‑3140498

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