Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,440,923 C→A 100% G36V (GGC→GTC)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,440,9230CA99.4% 2480.6 / 2.2 660G36V (GGC→GTC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base A (320/335);  total (324/336)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 660 total aligned reads displayed.

GGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAG  >  NC_000913/3440777‑3441069
                                                                                                                                                  |                                                                                                                                                  
 gTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACaaa                                                                                                                                                >  2:2765558/1‑150 (MQ=255)
  ttGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAg                                                                                                                                               <  2:6708234/150‑1 (MQ=255)
       cAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACg                                                                                                                                          <  2:1414975/150‑1 (MQ=255)
       cAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACg                                                                                                                                          >  2:2929842/1‑150 (MQ=255)
       cAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACg                                                                                                                                          >  2:4192989/1‑150 (MQ=255)
        aGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTGCCCAGAGTATGGACAAAGCCACGc                                                                                                                                         >  1:1253477/1‑150 (MQ=255)
             ttCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTaa                                                                                                                                    <  2:8129779/150‑1 (MQ=255)
             ttCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTaa                                                                                                                                    >  1:290831/1‑150 (MQ=255)
              tCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTaaa                                                                                                                                   <  1:1010490/150‑1 (MQ=255)
                cAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAgg                                                                                                                                 <  1:6870352/150‑1 (MQ=255)
                cAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAgg                                                                                                                                 <  1:6056758/150‑1 (MQ=255)
                 aGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGc                                                                                                                                <  1:2995220/150‑1 (MQ=255)
                         ttCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGt                                                                                                                        >  1:3574392/1‑150 (MQ=255)
                         ttCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGt                                                                                                                        >  2:6028946/1‑150 (MQ=255)
                          tCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGTTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTc                                                                                                                       >  2:2555335/1‑150 (MQ=255)
                          tCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTc                                                                                                                       <  2:1675328/150‑1 (MQ=255)
                           ccTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCa                                                                                                                      >  1:3092827/1‑150 (MQ=255)
                               aaCGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCtt                                                                                                                  >  2:2130984/1‑150 (MQ=255)
                                aCGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTg                                                                                                                 >  2:8811539/1‑150 (MQ=255)
                                 cGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTgg                                                                                                                <  1:1097833/150‑1 (MQ=255)
                                 cGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTgg                                                                                                                <  2:7293227/150‑1 (MQ=255)
                                  gCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                               <  1:2340570/150‑1 (MQ=255)
                                  gCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                               >  1:2601011/1‑150 (MQ=255)
                                  gCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                               >  1:6457802/1‑150 (MQ=255)
                                     ttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGGTCAAGGGTCACCTTGgcgtg                                                                                                            >  2:5474344/1‑150 (MQ=255)
                                         ttGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTc                                                                                                        >  1:449718/1‑150 (MQ=255)
                                            gTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGaa                                                                                                     <  1:4083962/150‑1 (MQ=255)
                                            gTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGaa                                                                                                     <  2:337472/150‑1 (MQ=255)
                                              gCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACt                                                                                                   >  1:799662/1‑150 (MQ=255)
                                                tGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCa                                                                                                 >  1:3315121/1‑150 (MQ=255)
                                                tGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCa                                                                                                 <  2:5743697/150‑1 (MQ=255)
                                                  tACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACt                                                                                               >  1:7295763/1‑150 (MQ=255)
                                                   aCTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACtt                                                                                              <  1:1874301/150‑1 (MQ=255)
                                                     tCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGc                                                                                            >  2:3985791/1‑150 (MQ=255)
                                                      cATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                           <  1:5612938/150‑1 (MQ=255)
                                                       aTGTAGTACACCATCATTCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGTGTCACCTTGGCGTGCGTCGAACTCACTTGCTc                                                                                          <  2:8362775/150‑1 (MQ=255)
                                                        tGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCg                                                                                         >  1:1357027/1‑150 (MQ=255)
                                                        tGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCg                                                                                         >  1:6592343/1‑150 (MQ=255)
                                                            gtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                                     <  1:7771365/150‑1 (MQ=255)
                                                             taCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATc                                                                                    <  1:4196578/150‑1 (MQ=255)
                                                                acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAc                                                                                 <  1:986006/150‑1 (MQ=255)
                                                                 ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc                                                                                <  2:5331930/150‑1 (MQ=255)
                                                                  cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCa                                                                               <  2:906148/150‑1 (MQ=255)
                                                                   aTCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAg                                                                              >  2:472318/1‑150 (MQ=255)
                                                                        tctcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                         >  2:3522296/1‑150 (MQ=255)
                                                                         ctcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCgg                                                                        <  1:1771102/150‑1 (MQ=255)
                                                                         ctcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCgg                                                                        >  2:8487759/1‑150 (MQ=255)
                                                                          tcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGt                                                                       <  1:6016048/150‑1 (MQ=255)
                                                                           cAACCTCGGTCACCGCGCAACCCAGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtt                                                                      >  1:7444273/1‑150 (MQ=255)
                                                                            aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGttt                                                                     <  2:8651132/150‑1 (MQ=255)
                                                                              ccTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTa                                                                   <  1:786310/150‑1 (MQ=255)
                                                                               cTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAg                                                                  <  2:1260779/150‑1 (MQ=255)
                                                                                  ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGaaa                                                               >  1:2240659/1‑150 (MQ=255)
                                                                                  ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGaaa                                                               <  2:7652325/150‑1 (MQ=255)
                                                                                  ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGaaa                                                               >  2:3671385/1‑150 (MQ=255)
                                                                                    tCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAAct                                                             >  2:7714781/1‑150 (MQ=255)
                                                                                        cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGt                                                         <  2:6872803/150‑1 (MQ=255)
                                                                                         gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTc                                                        <  2:6363969/150‑1 (MQ=255)
                                                                                         gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTc                                                        <  2:4638037/150‑1 (MQ=255)
                                                                                            cAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                     <  2:4786949/150‑1 (MQ=255)
                                                                                             aaCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAg                                                    >  1:3408184/1‑150 (MQ=255)
                                                                                             aaCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAg                                                    >  1:1528950/1‑150 (MQ=255)
                                                                                               cccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGaa                                                  >  2:6528883/1‑150 (MQ=255)
                                                                                               cccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGaa                                                  <  2:4246283/150‑1 (MQ=255)
                                                                                                ccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAc                                                 >  1:5822525/1‑150 (MQ=255)
                                                                                                ccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAc                                                 >  1:2755913/1‑150 (MQ=255)
                                                                                                   gCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCt                                              <  2:2392061/150‑1 (MQ=255)
                                                                                                       cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAt                                          <  1:3174316/150‑1 (MQ=255)
                                                                                                        gATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAtt                                         <  1:8135561/150‑1 (MQ=255)
                                                                                                        gATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAtt                                         >  1:2052973/1‑150 (MQ=255)
                                                                                                         aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtg                                        >  2:3211926/1‑150 (MQ=255)
                                                                                                         aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtg                                        <  1:8380756/150‑1 (MQ=255)
                                                                                                             gagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                    <  2:6129722/150‑1 (MQ=255)
                                                                                                                  gAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctctc                               <  1:9027234/150‑1 (MQ=255)
                                                                                                                  gAATACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctctc                               >  2:8067670/1‑150 (MQ=255)
                                                                                                                    aTACGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGGATTGTGTCCTCTCtt                             <  1:8298881/150‑1 (MQ=255)
                                                                                                                      aCGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTg                           <  1:858119/150‑1 (MQ=255)
                                                                                                                      aCGGCGCAGTGCGTTACCCAGAGTATGGAAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTg                           <  2:7875759/150‑1 (MQ=255)
                                                                                                                       cGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTgg                          <  2:7478332/150‑1 (MQ=255)
                                                                                                                       cGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTgg                          >  1:3819433/1‑150 (MQ=255)
                                                                                                                       cGGCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTgg                          <  1:3771730/150‑1 (MQ=255)
                                                                                                                        ggCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGt                         >  2:3800977/1‑150 (MQ=255)
                                                                                                                        ggCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGt                         <  1:2508633/150‑1 (MQ=255)
                                                                                                                        ggCGCAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGt                         >  1:4385415/1‑150 (MQ=255)
                                                                                                                         gcgcAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTa                        <  2:7033087/150‑1 (MQ=255)
                                                                                                                          cgcAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTGTTTGGTAc                       >  2:1115210/1‑150 (MQ=255)
                                                                                                                           gcAGTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACt                      >  1:5118873/1‑150 (MQ=255)
                                                                                                                              gTGCGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAg                   >  1:4583344/1‑150 (MQ=255)
                                                                                                                                 cGTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCtt                <  1:8811539/150‑1 (MQ=255)
                                                                                                                                  gTTACCCAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCttt               >  1:6211109/1‑150 (MQ=255)
                                                                                                                                       ccAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg          <  2:3024105/150‑1 (MQ=255)
                                                                                                                                       ccAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg          <  1:4925180/150‑1 (MQ=255)
                                                                                                                                       ccAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg          >  1:7127368/1‑150 (MQ=255)
                                                                                                                                        cAGAGTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg         >  2:498102/1‑150 (MQ=255)
                                                                                                                                         agagTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa        <  1:8879857/150‑1 (MQ=255)
                                                                                                                                         agagTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa        <  2:1817386/150‑1 (MQ=255)
                                                                                                                                         agagTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACCCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa        >  1:6306913/1‑150 (MQ=255)
                                                                                                                                            gTATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTa     <  1:5278784/150‑1 (MQ=255)
                                                                                                                                              aTGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTaaa   <  1:7001240/150‑1 (MQ=255)
                                                                                                                                               tGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAg  >  1:76579/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAG  >  NC_000913/3440777‑3441069

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