Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 922,316:1 +C 15.8% intergenic (+23/+50) macB → / ← cspD macrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913922,3141.C15.8% 480.4 / 90.9 139intergenic (+21/+52)macB/cspDmacrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base . (55/62);  new base C (7/14);  total (62/77)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 148 total aligned reads displayed.

CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTG  >  NC_000913/922215‑922413
                                                                                                     |                                                                                                   
cACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑g                                                                                                      >  1:2799551/1‑100 (MQ=255)
cACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑g                                                                                                      >  2:162451/1‑100 (MQ=255)
  ccGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑cc                                                                                                   >  2:1506186/1‑100 (MQ=255)
   cGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCg                                                                                                  >  2:1249507/1‑100 (MQ=255)
      gATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATc                                                                                               <  2:2647691/100‑1 (MQ=255)
       aTTTTATTTGGCTTGTTACCCGCACGAAATGCGGCACGAATGGATCCAGTAGATGCCCTGGCACGAGAGTAATTTTTGAGATAAAAATGCTA‑G‑CCg                                                                                                  >  1:2349206/1‑96 (MQ=255)
       aTTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCg                                                                                              >  2:2430405/1‑100 (MQ=255)
        ttttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCgg                                                                                             >  1:1215245/1‑100 (MQ=255)
        atttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCgg                                                                                             <  1:2857134/99‑1 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GCCCGATCggctgg                                                                                         <  1:1811852/100‑5 (MQ=255)
             tttGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTgg                                                                                        <  1:1689092/100‑1 (MQ=255)
              ttGGCCGGTTGGGCGCACGAGATGGGGCACTAGGGGATGCGGTAGATGCTCTGGCACGAGAGTAATTTCTGGGATAAAAGTGCCATC‑CCGATCggct                                                                                           <  1:1691463/97‑3 (MQ=255)
                ggctggctACCCGCACGCAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCa                                                                                      <  1:2479981/99‑1 (MQ=255)
                ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GCCCGATC‑GGCTGGCAt                                                                                     <  2:1107407/100‑1 (MQ=39)
                ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GCCCGATC‑GGCTGGCAt                                                                                     <  2:3291203/100‑1 (MQ=39)
                ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GCCCGATC‑GGCTGGCAt                                                                                     <  1:2216487/100‑1 (MQ=39)
                ggCTGATTACTCTCACGAAATGCGGCACGACTGGATTCAGTAGGTCATCTTGCACGAGAGTAAGTTTCTAGATAAAAATCACA‑C‑CCGATCGGGCTGGCa                                                                                      <  1:366145/99‑1 (MQ=255)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GCCCGATC‑GGCTGGCAtt                                                                                    <  1:764133/100‑1 (MQ=39)
                 gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GCCCGATC‑GGCTGGCAtt                                                                                    <  2:2261184/100‑1 (MQ=39)
                  cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCAttt                                                                                   >  2:3376863/1‑100 (MQ=255)
                  cTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCAttt                                                                                   <  1:2140636/100‑1 (MQ=255)
                   tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTTGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GCCCGATC‑GGCTGGCAtttt                                                                                  <  1:619610/100‑1 (MQ=38)
                   tGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCAtttt                                                                                  >  1:3377370/1‑100 (MQ=255)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTg                                                                                 >  2:2242365/1‑100 (MQ=255)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTg                                                                                 >  1:1164921/1‑100 (MQ=255)
                    ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGGAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCAtttt                                                                                  >  1:1649222/1‑99 (MQ=255)
                     gTTACCCGCACGAAATGCGGCACGACTGGATTGAGTAGATGCGCTGGCACGAGAGTAAGTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTg                                                                                 <  1:2637913/99‑1 (MQ=255)
                     gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGc                                                                                <  2:1776357/100‑1 (MQ=255)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GTCCGATC‑GGCTGGCAtttttat                                                                               <  1:1852613/100‑4 (MQ=38)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGcc                                                                               <  2:2444300/100‑1 (MQ=255)
                      ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGcc                                                                               >  1:3143107/1‑100 (MQ=255)
                       tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GCCCGATC‑GGCTGGCAtttttatc                                                                              <  2:3541994/100‑5 (MQ=38)
                        aCCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCtt                                                                             >  1:275227/1‑100 (MQ=255)
                          ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑GCCCGATC‑GGCTGGCAtttttatctca                                                                           <  2:1728115/100‑8 (MQ=37)
                          ccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTa                                                                           >  2:1397382/1‑100 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                          <  2:3121197/100‑1 (MQ=255)
                           cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                          <  1:2242613/100‑1 (MQ=255)
                            gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAgg                                                                         >  2:1954284/1‑100 (MQ=255)
                              aCGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGAt                                                                       <  1:3404073/100‑1 (MQ=255)
                              aCGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGGAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGAt                                                                       <  1:2721289/100‑1 (MQ=255)
                                gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGt                                                                     <  1:2909464/100‑1 (MQ=255)
                                 aaaTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAATAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTa                                                                    <  1:1378012/100‑1 (MQ=255)
                                     gCGGCACGACTGGATCCAGTAGATGCTCTGGAACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTac                                                                   <  2:1704898/97‑1 (MQ=255)
                                      cGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACaa                                                               >  2:2089974/1‑100 (MQ=255)
                                      cGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACaa                                                               >  2:295843/1‑100 (MQ=255)
                                       ggCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAAt                                                              >  2:3274147/1‑100 (MQ=255)
                                        gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATg                                                             <  1:3448226/100‑1 (MQ=255)
                                           cggcTGGATGCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgag                                                           <  1:3695404/96‑1 (MQ=255)
                                           cGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATgaga                                                          >  2:3597044/1‑100 (MQ=255)
                                            gACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGAc                                                         <  2:3590495/100‑1 (MQ=255)
                                             aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                        <  1:1640663/100‑1 (MQ=255)
                                             aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                        >  2:3354750/1‑100 (MQ=255)
                                             aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                        >  2:3062128/1‑100 (MQ=255)
                                             aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACa                                                        <  1:907458/100‑1 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                     >  1:1351597/1‑100 (MQ=255)
                                                ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaa                                                     >  2:3212369/1‑100 (MQ=255)
                                                 gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaag                                                    >  2:3817876/1‑100 (MQ=255)
                                                 gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaag                                                    <  2:3643327/100‑1 (MQ=255)
                                                   tCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCA‑TTTGCCTTTAGGATGTACACAATGAGACAGAAGAGc                                                 <  2:792945/100‑1 (MQ=255)
                                                     cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCt                                                >  2:1780336/1‑100 (MQ=255)
                                                       gTAGATGCTCTGGCACGAGAGTAATTTTTTAGTTATATATGCCA‑G‑CCGATCGGGCTGGCATTTTGACTTTAGGATGTACACAATGAGACAGAAGAGCTAt                                              >  1:271876/1‑100 (MQ=255)
                                                       gTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTAt                                              >  1:133231/1‑100 (MQ=255)
                                                         aGATGGTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTGAGGATGTACACAATGAGACAGAAGAGCTATga                                            <  1:1945478/100‑2 (MQ=255)
                                                         aGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGc                                            >  2:96694/1‑100 (MQ=255)
                                                            tGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGAc                                         <  2:1091662/100‑1 (MQ=255)
                                                                ctGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc                                     >  2:60367/1‑100 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                  <  1:409504/100‑1 (MQ=255)
                                                                   gCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                  <  2:62073/100‑1 (MQ=255)
                                                                    cACGAGAGTAATTTTTGAGATAATAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCtt                                 >  1:2192232/1‑100 (MQ=255)
                                                                     aCGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc                                <  1:3354730/100‑1 (MQ=255)
                                                                      cGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                               >  2:3180132/1‑100 (MQ=255)
                                                                        agagTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTAc                             >  2:2161524/1‑100 (MQ=255)
                                                                            tAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCg                         <  1:2588396/100‑1 (MQ=255)
                                                                              aTTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGAc                       <  1:1891645/100‑1 (MQ=255)
                                                                               tttttGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                      <  2:739743/100‑1 (MQ=255)
                                                                               tttttGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACg                      <  2:1644714/100‑1 (MQ=255)
                                                                                ttttGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACgg                     <  2:726533/100‑1 (MQ=255)
                                                                                ttttGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACgg                     <  2:188368/100‑1 (MQ=255)
                                                                                 tttGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGTGCTATGCGACTGCCGCTGCTACTTCGACgcg                    <  2:2570919/100‑3 (MQ=255)
                                                                                   tGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGca                  <  2:256073/100‑1 (MQ=255)
                                                                                    gagaTAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                 <  2:2032081/100‑1 (MQ=255)
                                                                                      gaTAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa               >  1:652213/1‑100 (MQ=255)
                                                                                      gaTAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa               >  2:773443/1‑100 (MQ=255)
                                                                                       aTAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAAt              <  1:216463/100‑1 (MQ=255)
                                                                                        tAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa             >  2:862005/1‑100 (MQ=255)
                                                                                        tAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATa             >  2:916181/1‑100 (MQ=255)
                                                                                        tAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACAGGCACAATa             >  1:113367/1‑100 (MQ=255)
                                                                                         aaaaaTGCCA‑G‑CCGATCGGGATGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa            >  1:2331889/1‑100 (MQ=255)
                                                                                           aaaTGCCA‑GCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGCCTGCCGCTTCTACTTCGACGGGCACAATAaca          >  1:497995/1‑100 (MQ=38)
                                                                                           aaaTGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAaca          >  2:2753300/1‑100 (MQ=255)
                                                                                           aaaTGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAaca          >  2:1381384/1‑100 (MQ=255)
                                                                                            aaTGCCA‑GCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAacac         >  2:2117234/1‑100 (MQ=39)
                                                                                             aTGCCA‑GCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACt        >  2:3545819/1‑100 (MQ=39)
                                                                                             aTGCCA‑GCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACt        >  2:1256866/1‑100 (MQ=39)
                                                                                              tGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg       >  2:2790600/1‑100 (MQ=255)
                                                                                               gCCA‑GCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg      >  2:679588/1‑100 (MQ=255)
                                                                                                 cA‑G‑CCGATCGGGCTGGCATTGTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg    <  1:1212373/100‑1 (MQ=255)
                                                                                                  a‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCAGCTTCTACTTCGACGGGCACAATAACACTGGCGt   >  1:3565478/1‑100 (MQ=255)
                                                                                                    g‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg  <  2:3116530/100‑1 (MQ=255)
                                                                                                    g‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTCCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg  <  1:2037421/100‑1 (MQ=255)
                                                                                                     |                                                                                                   
CACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCA‑G‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTG  >  NC_000913/922215‑922413

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