Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,533,536 A→G 100% A83A (GCT→GCC yjhQ ← GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,533,5360AG97.4% 471.8 / ‑0.3 148A83A (GCT→GCCyjhQGNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (93/51);  total (94/54)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

Only 100 of 150 total aligned reads displayed.

CCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGC  >  NC_000913/4533439‑4533632
                                                                                                  |                                                                                                
ccGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGc                                                                                                >  2:1826461/1‑100 (MQ=255)
  gTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGa                                                                                              >  2:2604568/1‑100 (MQ=255)
    cTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAgg                                                                                            >  1:707913/1‑100 (MQ=255)
     tGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGa                                                                                           >  1:2589066/1‑100 (MQ=255)
      ggCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAt                                                                                          >  2:383613/1‑100 (MQ=255)
      ggCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAt                                                                                          >  1:1788129/1‑100 (MQ=255)
      ggCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCCCCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAt                                                                                          >  1:3412780/1‑100 (MQ=255)
      ggCAGCCCATCAACCTCAGATGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAgg                                                                                            >  2:257417/1‑98 (MQ=255)
       gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAtg                                                                                         <  2:646607/100‑1 (MQ=255)
       gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAtg                                                                                         <  2:799193/100‑1 (MQ=255)
        cAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAtgt                                                                                        >  1:3190194/1‑100 (MQ=255)
         aGCCCATCAACCTCAGGTGTTCGATGCTAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAtgtg                                                                                       <  1:595868/100‑1 (MQ=255)
           cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                     <  1:1245291/100‑1 (MQ=255)
           cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                     >  1:2532391/1‑100 (MQ=255)
           cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTACCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                     <  2:362607/100‑1 (MQ=255)
            ccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCAt                                                                                    <  2:2029266/100‑1 (MQ=255)
             cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCCCGCCCATCCCTTTGTATTCCGGCATCACCGCCAGGGGGGCGAGGATGTGCATc                                                                                   >  1:1701439/1‑100 (MQ=255)
             cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCAGCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATc                                                                                   >  2:923565/1‑100 (MQ=255)
              aTCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCa                                                                                  <  1:3792076/100‑1 (MQ=255)
                 aaCCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACg                                                                               <  1:2421262/100‑1 (MQ=255)
                    cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                            <  1:2393154/100‑1 (MQ=255)
                    cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                            <  2:3871441/100‑1 (MQ=255)
                    cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                            <  1:3442563/100‑1 (MQ=255)
                       agttGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                            <  2:3084174/94‑1 (MQ=255)
                        ggTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATc                                                                        >  1:3882380/1‑100 (MQ=255)
                        ggTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATc                                                                        >  2:445099/1‑100 (MQ=255)
                         gTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATcc                                                                       >  1:151641/1‑100 (MQ=255)
                            cTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCAtt                                                                    <  2:2710024/100‑1 (MQ=255)
                            cTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGAATCAACAGTAAATCCAtt                                                                    >  2:3613105/1‑100 (MQ=255)
                              cGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTtc                                                                  >  1:3732127/1‑100 (MQ=255)
                                      gTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGa                                                          <  2:2763770/100‑1 (MQ=255)
                                         cGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAgg                                                       >  1:1814458/1‑100 (MQ=255)
                                         cGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAgg                                                       >  2:3195607/1‑100 (MQ=255)
                                          gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                                      >  2:1276259/1‑100 (MQ=255)
                                           tATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGtt                                                     >  1:3294723/1‑100 (MQ=255)
                                             tCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGc                                                   >  2:2965516/1‑100 (MQ=255)
                                              cAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGcc                                                  >  2:1933260/1‑100 (MQ=255)
                                              cAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGcc                                                  >  2:1934795/1‑100 (MQ=255)
                                              cAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGcc                                                  >  1:563355/1‑100 (MQ=255)
                                              cAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGcc                                                  <  1:904910/100‑1 (MQ=255)
                                               aaCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCGGGGGGGCGAGGATGGGCATAAACGGAGAAAACATTACTACTTAGAAgg                                                       >  1:3551131/1‑94 (MQ=255)
                                               aaCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGccc                                                 >  2:2161797/1‑100 (MQ=255)
                                                aCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGTCTTCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCg                                                <  1:469024/100‑1 (MQ=255)
                                                aCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCg                                                >  2:939404/1‑100 (MQ=255)
                                                aCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGAGAATCCATTTCTCCTTTGAAGGTTGCCCg                                                >  2:3009682/1‑100 (MQ=255)
                                                 ccGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCgg                                               >  1:500035/1‑100 (MQ=255)
                                                 ccGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCgg                                               >  2:2504045/1‑100 (MQ=255)
                                                 ccGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCgg                                               >  2:1683625/1‑100 (MQ=255)
                                                  cGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCggg                                              >  2:3547349/1‑100 (MQ=255)
                                                   gTCCTCCCACGCCCATCCCTTGTTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCggg                                              <  2:3513563/99‑1 (MQ=255)
                                                   gTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGt                                             >  2:590455/1‑100 (MQ=255)
                                                    tcctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                                            >  2:1540599/1‑100 (MQ=255)
                                                    tcctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                                            <  2:1038964/100‑1 (MQ=255)
                                                     cctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaa                                           >  2:1652190/1‑100 (MQ=255)
                                                      ctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaa                                          >  2:3738367/1‑100 (MQ=255)
                                                      ctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaa                                          >  1:448650/1‑100 (MQ=255)
                                                      ctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTATCCTTTGAAGGTTGCCCGGGtaaa                                          >  1:441584/1‑100 (MQ=255)
                                                        cccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaata                                        >  1:2943253/1‑100 (MQ=255)
                                                        cccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGaaaata                                        >  1:1767250/1‑100 (MQ=255)
                                                         ccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaataa                                       <  2:2957072/100‑1 (MQ=255)
                                                          cACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaataaa                                      <  2:1788764/100‑1 (MQ=255)
                                                          cACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaataaa                                      >  1:2245839/1‑100 (MQ=255)
                                                          cACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGCGGGGAGAGGATGTGCATTAACGGTGAATCCATTTCTCCTTTGAAGGTTGACCGGGtcaatcaa                                      >  2:778379/1‑100 (MQ=255)
                                                             gCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaata                                   >  2:1148615/1‑100 (MQ=255)
                                                              cccATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAatat                                  >  2:1696921/1‑100 (MQ=255)
                                                                cATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGc                                >  2:1691479/1‑100 (MQ=255)
                                                                 aTCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGcc                               <  1:3369015/100‑1 (MQ=255)
                                                                 aTCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCAAGGAAAATAAAATATGc                                >  1:1619388/1‑99 (MQ=255)
                                                                   cccTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg                             >  2:1326553/1‑100 (MQ=255)
                                                                        ggTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGc                        >  2:3481369/1‑100 (MQ=255)
                                                                         gTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCt                       <  1:3051094/100‑1 (MQ=255)
                                                                          tATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCtt                      >  2:2242014/1‑100 (MQ=255)
                                                                           aTTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCttt                     <  2:1537151/100‑1 (MQ=255)
                                                                             tCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGc                   >  1:1550581/1‑100 (MQ=255)
                                                                              ccGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGcc                  >  2:3376456/1‑100 (MQ=255)
                                                                              ccGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGcc                  <  1:3493990/100‑1 (MQ=255)
                                                                               cGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCATTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCt                 <  2:3461305/100‑1 (MQ=255)
                                                                               cGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCt                 <  1:1025719/100‑1 (MQ=255)
                                                                                ggTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGTATCCATTTCTCCTTTGAAGGTTGTCCGGGTAAATAAAATATGACCGACGGCTTTACCtt                >  1:1728227/1‑100 (MQ=255)
                                                                                 gTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGTAGGTGGCCCGGGTAAATAAAATATTCCCGACGGCTTTGcat                 >  2:2444120/1‑96 (MQ=255)
                                                                                  tATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCa              <  1:2242034/100‑1 (MQ=255)
                                                                                  tATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGCGTAAATAAAATATGCCCGACGGCTTTGCCTTCa              <  2:2589042/100‑1 (MQ=255)
                                                                                    tGACCGCCAGGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCAt             <  2:2289958/100‑1 (MQ=255)
                                                                                      aCCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATAcc          >  2:3537606/1‑100 (MQ=255)
                                                                                      aCCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATAcc          >  1:3267479/1‑100 (MQ=255)
                                                                                      aCCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATAcc          >  2:12246/1‑100 (MQ=255)
                                                                                       ccGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCg         >  2:2072541/1‑100 (MQ=255)
                                                                                       ccGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTACCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCg         >  1:3137901/1‑100 (MQ=255)
                                                                                        cGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCgg        >  1:2854295/1‑100 (MQ=255)
                                                                                         gCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCggg       >  2:1318995/1‑100 (MQ=255)
                                                                                         gCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCggg       <  1:883916/100‑1 (MQ=255)
                                                                                         gCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCggg       >  1:2616760/1‑100 (MQ=255)
                                                                                          ccAGGGGGGCGAGGATGTGCATCAACGTTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGc      <  1:1901323/100‑1 (MQ=255)
                                                                                           cAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCt     <  2:3294689/100‑1 (MQ=255)
                                                                                           cAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCt     <  1:3835951/100‑1 (MQ=255)
                                                                                            aGGGGGGCTAGTATGTGCATCTACTGTGAATCCATTTCCCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCTACGGCTTTGCCTTCATACCggg       <  2:3351283/97‑1 (MQ=255)
                                                                                            aGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTa    >  1:675512/1‑100 (MQ=255)
                                                                                            aGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTa    >  2:3375664/1‑100 (MQ=255)
                                                                                              ggggggCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTag   >  1:255998/1‑100 (MQ=255)
                                                                                               ggggGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTagc  <  1:1156660/100‑1 (MQ=255)
                                                                                                  |                                                                                                
CCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGC  >  NC_000913/4533439‑4533632

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