Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,481 C→T 6.6% intergenic (+27/+19) umuC → / ← dsbB translesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,232,4810CT 6.6% 299.3 / 7.0 105intergenic (+27/+19)umuC/dsbBtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base C (45/52);  new base T (3/4);  total (48/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

Only 100 of 111 total aligned reads displayed.

AAAGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACT  >  NC_000913/1232332‑1232630
                                                                                                                                                     |                                                                                                                                                     
aaaGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                                                                       <  2:572561/150‑1 (MQ=255)
aaaGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                                                                       <  2:564613/150‑1 (MQ=255)
aaaGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                                                                       >  1:1934124/1‑150 (MQ=255)
    agGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAg                                                                                                                                                   <  2:1554642/150‑1 (MQ=255)
     ggggAGATGAACACTTTATTTTGCCCGGCAGCTGATCCAGCCACAATGGCAGATTAAGCGACCCATTCTTTCACCACGTTATACAACGGGAAGTTCTGATTTACTGATGGTCAAATATATATAGCGGCACGAAAAAAGCGCTCCCGCAgg                                                                                                                                                  <  1:1919562/150‑1 (MQ=255)
          gagGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACTCGAAGTTCTGATTTACTGAGGGTCAAATAAATAGAGCGGCAGGATAAAAGCGCTCCCGCAGGAgcgc                                                                                                                                             >  2:323340/1‑150 (MQ=255)
            ggAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctt                                                                                                                                           <  2:1953536/150‑3 (MQ=255)
             gAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGATATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGa                                                                                                                                          >  2:940999/1‑150 (MQ=255)
              aaCACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACACCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                                                                         <  1:1884829/150‑5 (MQ=255)
              aaCACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGaa                                                                                                                                         >  1:1546202/1‑150 (MQ=255)
              aaCACTCTATTTAGCCGGGCAGGGGATCCAGCAACACTGCCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGAATTACTGAGTGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                                                                         <  2:212168/150‑5 (MQ=255)
                cacTCTATTTTGCCGGGCAGGGGATCCATCAACAATGGCAGATGACGCGAGCCATGCTTTCACCACGTAATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                                                                       <  1:1769361/150‑1 (MQ=255)
                cacTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                                                                       <  2:1569483/150‑1 (MQ=255)
                 acTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttc                                                                                                                                      <  2:1516937/150‑8 (MQ=255)
                   tctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                                                                    >  1:362863/1‑150 (MQ=255)
                   tctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                                                                    <  2:1648865/150‑1 (MQ=255)
                   tctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                                                                    >  1:1275597/1‑150 (MQ=255)
                   tctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                                                                    <  2:736691/150‑1 (MQ=255)
                       ttttGCCGGGCAGGGGATCCAGCAACAATGGCAGATGACGCGAGCCATTCTTTCACCACGTTATACAACGCGAAGTTCTGAGTTACTGAGGGTCAAATAAATATAGCGTCAGGAAAAAAGCTCTCCCTCAGGAGCGCCGAATGGAtgatc                                                                                                                                >  2:1224827/1‑146 (MQ=255)
                       ttttGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAAGGGATTAGc                                                                                                                                >  2:100675/1‑150 (MQ=255)
                       ttatGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACACCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGCCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGc                                                                                                                                >  1:779731/4‑150 (MQ=255)
                           gCCGGGCAGTGGATCCAGCAACAATCGCAGAAGCAGCGACCCATGCTATCACCACCTAATACAACGCGAAGTTCTGATTTACCGAGTGTCAAAAAAATATAGCGGCAGGAAACAACCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAcc                                                                                                                            <  1:1243362/150‑1 (MQ=255)
                            ccGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGATCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGTAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                                                                           >  1:1217045/1‑150 (MQ=255)
                              gggCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGaa                                                                                                                         <  1:1160329/150‑1 (MQ=255)
                              gggCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGAGTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGaa                                                                                                                         >  1:280276/1‑150 (MQ=255)
                                 cAGGGGATCCAGCAAGAATGTCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGCGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAg                                                                                                                      <  2:405935/150‑1 (MQ=255)
                                     ggATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCa                                                                                                                  >  2:1745076/1‑150 (MQ=255)
                                      gATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAc                                                                                                                 >  1:1642960/1‑150 (MQ=255)
                                        tCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGt                                                                                                               >  1:726034/1‑150 (MQ=255)
                                          cAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttt                                                                                                             >  2:1041310/1‑150 (MQ=255)
                                            tcaacaaTGGCAGATGAAGCGAGGCATGATTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCTGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTTGCGACCGATCAGATCACGttttt                                                                                                           <  2:450330/149‑1 (MQ=255)
                                             caacaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttttt                                                                                                          >  2:1971171/1‑150 (MQ=255)
                                             caaaaaTGGCAGATGAAGCGAGCCATGCTTTCACAACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                                                                         <  1:390214/150‑1 (MQ=255)
                                               acaaTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCg                                                                                                        >  2:1763977/1‑150 (MQ=255)
                                                   tGGCAGTTGAAGCGAGCCATGGTTTCACCACGTTTAACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCttt                                                                                                    <  2:90754/150‑1 (MQ=255)
                                                   tGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCttt                                                                                                    <  2:848968/150‑1 (MQ=255)
                                                       aGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGTAGCGCCGATTGGTTTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAAc                                                                                                >  1:1784211/1‑150 (MQ=255)
                                                       aGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAAc                                                                                                <  2:230973/150‑1 (MQ=255)
                                                            aaGCGAGCCATTCTTTCACCAAGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                                                                           <  2:1244065/150‑1 (MQ=255)
                                                            aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                                                                           <  2:1837966/150‑1 (MQ=255)
                                                             aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTgg                                                                                          <  2:492679/150‑1 (MQ=255)
                                                              gCGAGCTATGCTTTCACCACGTTATAGAACTCGAAGTTCTGATTTACTGAGGGGCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGCGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                                                                         <  1:223170/150‑1 (MQ=255)
                                                                gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTTCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaa                                                                                       <  1:189523/150‑1 (MQ=255)
                                                                gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaa                                                                                       <  1:31548/150‑1 (MQ=255)
                                                                 aGCCATGCTTTCACCACCTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaaa                                                                                      <  1:1804873/150‑1 (MQ=255)
                                                                    cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                                                                   <  1:963964/150‑1 (MQ=255)
                                                                      tGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcacc                                                                                 >  1:1693030/1‑150 (MQ=255)
                                                                          ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAccacca                                                                             >  1:354435/1‑150 (MQ=255)
                                                                            caccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGCTCACGTTTTTTCGCTTTAAACGGCTGGGAATTCACCACCAGc                                                                           >  2:1975570/1‑150 (MQ=255)
                                                                             accacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCa                                                                          >  2:1104385/1‑150 (MQ=255)
                                                                             accacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCa                                                                          >  1:112777/1‑150 (MQ=255)
                                                                              ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGTGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                                                                         >  2:287999/1‑150 (MQ=255)
                                                                                  gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                                                                     <  2:1708790/150‑1 (MQ=255)
                                                                                   ttATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGa                                                                    <  2:1340739/150‑1 (MQ=255)
                                                                                   ttATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGa                                                                    <  2:205406/150‑1 (MQ=255)
                                                                                    tataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGAGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                                                                   <  1:110594/150‑1 (MQ=255)
                                                                                    tataCAACGCGAAGATCTGATTTACTGAGGGTCAAATAAATATTGCTGCAGGAAAAACGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                                                                   <  1:862738/150‑1 (MQ=255)
                                                                                     ataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACa                                                                  <  1:1446902/150‑1 (MQ=255)
                                                                                      taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTTCGACaa                                                                 >  2:1426560/1‑150 (MQ=255)
                                                                                       aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                                                                >  1:672629/1‑150 (MQ=255)
                                                                                        cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCCCCAGCACTGCGACAATc                                                               >  1:1359951/1‑150 (MQ=255)
                                                                                        cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc                                                               >  2:911174/1‑150 (MQ=255)
                                                                                           cgcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg                                                            <  1:1519020/150‑1 (MQ=255)
                                                                                            gcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGt                                                           <  1:903239/150‑1 (MQ=255)
                                                                                                aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc                                                       <  2:1546202/150‑1 (MQ=255)
                                                                                                    cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATa                                                   >  1:956651/1‑150 (MQ=255)
                                                                                                      gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaa                                                 >  2:1232152/1‑150 (MQ=255)
                                                                                                            cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                           >  1:1919439/1‑150 (MQ=255)
                                                                                                            cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGAGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGCAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                           >  1:1120269/1‑150 (MQ=255)
                                                                                                               tgggTCAAATATATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGag                                        <  1:173787/149‑1 (MQ=255)
                                                                                                                gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagc                                       <  1:555698/150‑1 (MQ=255)
                                                                                                                gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCAGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagc                                       <  2:355905/150‑1 (MQ=255)
                                                                                                                   tCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCAGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagcagc                                    <  1:323340/150‑1 (MQ=255)
                                                                                                                    cAAATAAATATAGCGGCAGGAAAAAAGCTCTCCCGCAGGAGCGCCGAATGGATTATCGAACGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCAACAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                   <  2:748228/150‑1 (MQ=255)
                                                                                                                    cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGATTGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTTCTACAATCAGGTAAGCGATAAAAAAACCGAGCAGcc                                   >  2:1206767/1‑150 (MQ=255)
                                                                                                                    cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                   >  2:307286/1‑150 (MQ=255)
                                                                                                                    aatataaataTAGCGGCAGGAAAAAAGCTGTCTCTCAGGAGAGCCGAATGGATTAGCGACCGATCAGATCACGTTTTTTCGCTTTAAACGGCTGGGAATTCACCACCAGCAGTGCTACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                   <  2:1826928/147‑1 (MQ=255)
                                                                                                                    aaaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGTCCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                   <  1:1241912/149‑1 (MQ=255)
                                                                                                                           atatAGCGGCCGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCgg                            <  2:145615/150‑1 (MQ=255)
                                                                                                                            tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCCCCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGc                           >  2:1909312/1‑150 (MQ=255)
                                                                                                                                  ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGAACAGATCACGTTTTTTCGCATTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                     >  2:1236246/1‑150 (MQ=255)
                                                                                                                                  ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                     <  1:940999/150‑1 (MQ=255)
                                                                                                                                  ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                     <  1:1380419/150‑1 (MQ=255)
                                                                                                                                   gCAGGAAAAAAGCGCTCCCGCAGGAGAGCCGATTGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACTATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCa                    >  2:44319/1‑150 (MQ=255)
                                                                                                                                    cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGTATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAACACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGATCAGCCACTCCGGCATTTCCac                   >  2:1749864/1‑149 (MQ=255)
                                                                                                                                    cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAAGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAg                   >  2:1441691/1‑150 (MQ=255)
                                                                                                                                     aGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGa                  <  2:1693552/150‑1 (MQ=255)
                                                                                                                                     aGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGa                  >  2:1623530/1‑150 (MQ=255)
                                                                                                                                      ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCTCCACGAGGACTGCGTCAATCTGGTAAGCGATAATAATAGCGAGCAGGCAGTGCGGCATTTCCAGAc                 <  2:1745249/150‑1 (MQ=255)
                                                                                                                                      ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAc                 <  2:892067/150‑1 (MQ=255)
                                                                                                                                       gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAcc                <  1:1464156/150‑1 (MQ=255)
                                                                                                                                       gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAAGTTTTTTCGCTTTAAAAGGCTGGGAAATCACAACCAGCACTGCGCCAATCAGGTAAGCGAAAAAAATACCGAGCAGCCACTGCGCCATTTCCAGAcc                >  2:296286/1‑150 (MQ=255)
                                                                                                                                         attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTa              >  1:472182/6‑150 (MQ=255)
                                                                                                                                         aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTa              >  2:402878/1‑150 (MQ=255)
                                                                                                                                          ttcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACGGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaa             >  2:1952943/5‑150 (MQ=255)
                                                                                                                                            cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCCCTGCGGCATTTCCAGACCTaaaa           >  1:1261443/3‑150 (MQ=255)
                                                                                                                                              gcgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaaaaa         <  2:482970/150‑1 (MQ=255)
                                                                                                                                               cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAAt        <  1:1290672/150‑1 (MQ=255)
                                                                                                                                                   cccGCAGGAGAGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGGAAGCGATAAAAATACCGAGCAGCCACTCCGGCATGTCCAGACCAAAAGAATcccc    >  1:980138/1‑149 (MQ=255)
                                                                                                                                                     cGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACt  <  1:1915102/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAGAGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACT  >  NC_000913/1232332‑1232630

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