Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,125,966 C→A T135T (ACC→ACA rimJ → ribosomal‑protein‑S5‑alanine N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,125,9660CA97.7% 308.2 / ‑0.4 88T135T (ACC→ACArimJribosomal‑protein‑S5‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (45/41);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

ATGCCTGCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACCCAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGCTATGC  >  NC_000913/1125869‑1126061
                                                                                                 |                                                                                               
aTGCCTGCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGcacaca                                                                                               <  1:677821/100‑1 (MQ=255)
   ccTGCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACacaaca                                                                                            <  1:1780951/100‑1 (MQ=255)
   ccTGCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACacaaca                                                                                            >  1:1619134/1‑100 (MQ=255)
    cTGCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACacaaca                                                                                            <  1:627711/99‑1 (MQ=255)
    cTGCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACat                                                                                           >  1:1004338/1‑100 (MQ=255)
     tGCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACata                                                                                          <  2:1940801/100‑1 (MQ=255)
     tGCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACata                                                                                          <  2:1024445/100‑1 (MQ=255)
      gCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACatat                                                                                         <  2:380999/100‑1 (MQ=255)
      gCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACatat                                                                                         >  1:1875558/1‑100 (MQ=255)
      gCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACata                                                                                          >  1:723675/1‑99 (MQ=255)
       cTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATAtt                                                                                        <  2:1468919/100‑1 (MQ=255)
       cTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATAtt                                                                                        <  2:1469432/100‑1 (MQ=255)
        tATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTc                                                                                       >  2:391912/1‑100 (MQ=255)
         aTCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCa                                                                                      <  2:277599/100‑1 (MQ=255)
           ctcGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATc                                                                                    >  2:1568397/1‑100 (MQ=255)
            tcGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCg                                                                                   <  1:649749/100‑1 (MQ=255)
             cGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGc                                                                                  >  2:1552771/1‑100 (MQ=255)
             cGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGc                                                                                  >  1:1880569/1‑100 (MQ=255)
               gTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTTATCGCAt                                                                                >  1:2282103/1‑100 (MQ=255)
               gTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCAt                                                                                >  2:1408933/1‑100 (MQ=255)
               gTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCAt                                                                                <  2:2197677/100‑1 (MQ=255)
                  aTTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTAt                                                                             <  1:1552775/100‑1 (MQ=255)
                   ttCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATg                                                                            <  1:80824/100‑1 (MQ=255)
                   ttCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATg                                                                            <  2:512331/100‑1 (MQ=255)
                   ttCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGACCTGCCGGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATg                                                                            <  2:1013225/100‑1 (MQ=255)
                     cGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGc                                                                          >  2:783384/1‑100 (MQ=255)
                      gATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCt                                                                         >  2:1528393/1‑100 (MQ=255)
                          gggCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTTATCGCATTATGGCTAAtt                                                                     >  2:545029/1‑100 (MQ=255)
                          gggCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAAtt                                                                     >  1:1074121/1‑100 (MQ=255)
                           ggCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATta                                                                    >  1:1049444/1‑100 (MQ=255)
                            gCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATtat                                                                   >  2:2324010/1‑100 (MQ=255)
                              aaaaaTGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTatat                                                                 <  1:656886/100‑1 (MQ=255)
                               aaaaTGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATg                                                                <  2:1162926/100‑1 (MQ=255)
                                taaTGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGc                                                               <  2:947692/99‑1 (MQ=255)
                                 aaTGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGcc                                                              <  2:1619129/100‑1 (MQ=255)
                                 aaTGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGca                                                              >  1:947694/1‑99 (MQ=255)
                                    ggCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGca                                                           >  2:1895893/1‑100 (MQ=255)
                                    ggCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGca                                                           >  2:721384/1‑100 (MQ=255)
                                    ggCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGca                                                           >  1:1016964/1‑100 (MQ=255)
                                    ggCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGca                                                           >  1:1336503/1‑100 (MQ=255)
                                     gCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGcac                                                          >  2:2090859/1‑100 (MQ=255)
                                     gCAGGGCAAAGCACTTATGTTTGAAGCCCTGACTGCCGCCATTAGTTATATGCAGCGCACACAACATATTCTTCGCATTATGGCTAACTATGTGCTGca                                                           <  2:1229444/99‑1 (MQ=255)
                                          gCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATaa                                                     <  2:1337446/100‑1 (MQ=255)
                                              aGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGc                                                 <  2:438680/100‑1 (MQ=255)
                                                 aCTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGc                                              >  1:1557833/1‑100 (MQ=255)
                                                  cTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCg                                             >  2:907731/1‑100 (MQ=255)
                                                  cTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCg                                             >  2:2112424/1‑100 (MQ=255)
                                                    cATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGt                                           >  1:2308294/1‑100 (MQ=255)
                                                      tGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGa                                         <  1:2324019/100‑1 (MQ=255)
                                                        tttGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGAtt                                       >  1:1058198/1‑100 (MQ=255)
                                                        tttGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGAtt                                       <  2:510521/100‑1 (MQ=255)
                                                          tGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTa                                     <  1:391914/100‑1 (MQ=255)
                                                             aGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTg                                  <  1:721385/100‑1 (MQ=255)
                                                             aGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTg                                  <  1:1895900/100‑1 (MQ=255)
                                                             aGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTg                                  <  1:1568402/100‑1 (MQ=255)
                                                             aGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTg                                  <  2:692202/100‑1 (MQ=255)
                                                              gCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTgg                                 <  1:1122531/100‑1 (MQ=255)
                                                              gCCCTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAAGTATATGCCGCACAATAAACGCAGCGGTGATTTACTgg                                 >  1:105971/1‑100 (MQ=255)
                                                                ccTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTTATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGgcg                               <  1:545031/100‑1 (MQ=255)
                                                                ccTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGgcg                               >  1:2327001/1‑100 (MQ=255)
                                                                ccTGACCGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGgcg                               <  2:650015/100‑1 (MQ=255)
                                                                     ccGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACt                          <  1:1214719/100‑1 (MQ=255)
                                                                      cGCAGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTg                         >  2:336385/1‑100 (MQ=255)
                                                                         aGCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGt                      >  1:1885585/1‑100 (MQ=255)
                                                                          gCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGtt                     >  1:2046200/1‑100 (MQ=255)
                                                                          gCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGtt                     >  1:2170031/1‑100 (MQ=255)
                                                                          gCCATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGtt                     <  2:1016962/100‑1 (MQ=255)
                                                                            cATTCGTTATATGCCGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCCGCGGTGATTTACTGGCGCGACTGGGtttt                   >  2:476820/1‑100 (MQ=255)
                                                                            cATTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGtttt                   >  1:497865/1‑100 (MQ=255)
                                                                             aTTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTg                  <  1:2090866/100‑1 (MQ=255)
                                                                             aTTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTg                  <  2:1520787/100‑1 (MQ=255)
                                                                             aTTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTg                  <  2:1193758/100‑1 (MQ=255)
                                                                             aTTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTg                  <  2:2282094/100‑1 (MQ=255)
                                                                             aTTCGTTATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTATTg                  <  1:1408936/100‑1 (MQ=255)
                                                                                  ttATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGaaaaa             <  1:1178891/100‑1 (MQ=255)
                                                                                  ttATATGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGaaaaa             >  1:562709/1‑100 (MQ=255)
                                                                                   tataTGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAaag            <  2:1875552/100‑1 (MQ=255)
                                                                                   tataTGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAaag            <  1:907733/100‑1 (MQ=255)
                                                                                   natatGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAaag            >  2:443636/2‑100 (MQ=255)
                                                                                    atatGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAaaga           >  2:720332/1‑100 (MQ=255)
                                                                                     tatGCAGCGCACCCAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAaagaa          >  2:352223/1‑100 (MQ=255)
                                                                                       tGCAGCGCACCCAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAgg        <  1:352225/100‑1 (MQ=255)
                                                                                       tGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAgg        >  1:1319925/1‑100 (MQ=255)
                                                                                       tGCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAgg        >  1:83690/1‑100 (MQ=255)
                                                                                        gCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGc       >  2:1861966/1‑100 (MQ=255)
                                                                                        gCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGc       >  2:550721/1‑100 (MQ=255)
                                                                                        gCAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGc       >  2:984100/1‑100 (MQ=255)
                                                                                         cAGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGCt      >  2:1840965/1‑100 (MQ=255)
                                                                                          aGCGCACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGCTa     >  1:737233/1‑100 (MQ=255)
                                                                                            cgcACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGCTATg   <  2:2046193/100‑1 (MQ=255)
                                                                                             gcACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGCTATGc  >  1:1274186/1‑100 (MQ=255)
                                                                                             gcACACAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGCTATGc  <  1:1022616/100‑1 (MQ=255)
                                                                                                 |                                                                                               
ATGCCTGCTATCTCGGTTATTCGATTGGGCAAAAATGGCAGGGCAAAGGACTCATGTTTGAAGCCCTGACCGCAGCCATTCGTTATATGCAGCGCACCCAACATATTCATCGCATTATGGCTAATTATATGCCGCACAATAAACGCAGCGGTGATTTACTGGCGCGACTGGGTTTTGAAAAAGAAGGCTATGC  >  NC_000913/1125869‑1126061

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