Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,484 A→G 6.3% intergenic (+30/+16) 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,4840AG 6.3% 332.9 / 9.6 110intergenic (+30/+16)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 A (48/55);  new base G (3/4);  total (51/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 111 total aligned reads displayed.

AGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGAC  >  NC_000913/1232336‑1232633
                                                                                                                                                    |                                                                                                                                                     
agGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAg                                                                                                                                                      <  2:1554642/150‑1 (MQ=255)
 ggggAGATGAACACTTTATTTTGCCCGGCAGCTGATCCAGCCACAATGGCAGATTAAGCGACCCATTCTTTCACCACGTTATACAACGGGAAGTTCTGATTTACTGATGGTCAAATATATATAGCGGCACGAAAAAAGCGCTCCCGCAgg                                                                                                                                                     <  1:1919562/150‑1 (MQ=255)
   gCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAg                                                                                                                                                   <  2:138312/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:1275597/1‑150 (MQ=255)
               tctATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                                                                       >  1:362863/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)
                           ggCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                                                                           <  2:1672679/150‑1 (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)
                                                  cAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGATCAGATCACGTTTTTTCGCTTTaaa                                                                                                    >  1:619462/1‑150 (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                                                                             >  1:112777/1‑150 (MQ=255)
                                                                         accacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCa                                                                             >  2:1104385/1‑150 (MQ=255)
                                                                              gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                                                                        <  2:1708790/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)
                                                                                       cgagAAGTTCTGATTTACTGAGGGACAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg                                                               <  2:515752/147‑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)
                                                                                                      tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAAtc                                                >  1:987833/1‑149 (MQ=255)
                                                                                                        cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                              >  1:1919439/1‑150 (MQ=255)
                                                                                                        cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGAGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGCAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATAcc                                              >  1:1120269/1‑150 (MQ=255)
                                                                                                         tGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCg                                             >  2:1915102/1‑150 (MQ=255)
                                                                                                           tgggTCAAATATATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGag                                           <  1:173787/149‑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)
                                                                                                                              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:1623530/1‑150 (MQ=255)
                                                                                                                                 aGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGa                     <  2:1693552/150‑1 (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)
                                                                                                                                           cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAAt           <  1:1290672/150‑1 (MQ=255)
                                                                                                                                               cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCAGTACCAGACCTAAAAACTcccc       >  2:301241/1‑149 (MQ=255)
                                                                                                                                                 cGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACt     <  1:1915102/150‑1 (MQ=255)
                                                                                                                                                  gCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCCTTTCCAGACCTAAAAAATCCCACTg    >  1:335598/1‑150 (MQ=255)
                                                                                                                                                   cAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGa   <  2:1934124/150‑1 (MQ=255)
                                                                                                                                                    aGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGAc  <  1:1763977/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGGCAGAGGAACACTCTATTTTGCCGGGCAGGGGATCCAGCAACAATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGAC  >  NC_000913/1232336‑1232633

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: