Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 13,957 A→T M599L (ATG→TTG)  dnaK → chaperone Hsp70, with co‑chaperone DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091313,9570AT96.5% 287.0 / 0.8 87M599L (ATG→TTG) dnaKchaperone Hsp70, with co‑chaperone DnaJ
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (41/43);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAAT  >  NC_000913/13860‑14053
                                                                                                 |                                                                                                
gTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTg                                                                                                <  1:342949/100‑1 (MQ=255)
gTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTg                                                                                                <  2:2355660/100‑1 (MQ=255)
 tCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGTTTTCCCAGAAATTGTTgg                                                                                               <  1:2204655/100‑1 (MQ=255)
    gcgcTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGaaa                                                                                            <  1:1033549/100‑1 (MQ=255)
     cgcTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAAt                                                                                           >  2:2263154/1‑100 (MQ=255)
      gcTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATc                                                                                          >  2:833286/1‑100 (MQ=255)
       cTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCg                                                                                         <  1:1269148/100‑1 (MQ=255)
          ccTGCACTGGAAACTGCTATGAAAGGTGAAGATAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGccc                                                                                      <  1:123228/99‑1 (MQ=255)
             gCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCcagc                                                                                   >  2:163755/1‑100 (MQ=255)
             gCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCcagc                                                                                   >  2:1431450/1‑100 (MQ=255)
             gCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCcagc                                                                                   <  1:116613/100‑1 (MQ=255)
             gCACCGGAAACTACTCTGAAAGGGGAAGACAAAACCGCCATCGAAGCGAAAATGCAGGAACCGGCCCAGGGTTCCCAGAAACCGTTGGAAAAAGCCcagc                                                                                   >  2:1121928/1‑100 (MQ=255)
              cACTGGAAACTGCTCTGAAAGGTGAAGACAAGGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCcagca                                                                                  >  1:1144011/1‑100 (MQ=255)
              cACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCcagca                                                                                  >  1:1388700/1‑100 (MQ=255)
              cACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCcagca                                                                                  <  1:733335/100‑1 (MQ=255)
              cACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCcagca                                                                                  <  1:2062080/100‑1 (MQ=255)
               aCTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCcagcag                                                                                 <  1:1602788/100‑1 (MQ=255)
                cTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCcagatcg                                                                                >  1:2307033/1‑96 (MQ=255)
                  ggAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCaa                                                                              >  1:337254/1‑100 (MQ=255)
                      aCTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATg                                                                          <  1:1162540/100‑1 (MQ=255)
                         gCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGccc                                                                       >  2:1664273/1‑100 (MQ=255)
                          ctctGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCca                                                                      <  1:2581542/100‑1 (MQ=255)
                          ctctGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTATCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATACCca                                                                      >  2:1827393/1‑100 (MQ=255)
                             tGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCcagca                                                                   >  1:1885/1‑100 (MQ=255)
                              gAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCcggcag                                                                  <  2:2196784/100‑1 (MQ=255)
                              gAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCcagcag                                                                  >  1:2066743/1‑100 (MQ=255)
                              gAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCcagcag                                                                  <  1:910972/100‑1 (MQ=255)
                                  ggTGAAGACGAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCACCAGa                                                              >  2:2519809/1‑100 (MQ=255)
                                  ggTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGa                                                              >  1:757298/1‑100 (MQ=255)
                                   gTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGAc                                                             >  1:2543547/1‑100 (MQ=255)
                                   gTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGAc                                                             >  1:1784353/1‑100 (MQ=255)
                                    tGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACt                                                            >  1:985417/1‑100 (MQ=255)
                                    tGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACt                                                            >  2:1628785/1‑100 (MQ=255)
                                       aGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGcc                                                         >  1:927642/1‑100 (MQ=255)
                                       aGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGcc                                                         <  1:617479/100‑1 (MQ=255)
                                        gACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCg                                                        <  1:1104452/100‑1 (MQ=255)
                                         aCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCgg                                                       >  2:2116500/1‑100 (MQ=255)
                                         aCAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCgg                                                       <  2:176295/100‑1 (MQ=255)
                                               ccGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGa                                                 <  2:607392/100‑1 (MQ=255)
                                               ccGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGa                                                 <  2:581809/100‑1 (MQ=255)
                                                 gCTATCGAAGCGAAACTGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATACCCAGCAGCAGACTGCCGGTGCTGATg                                               <  2:195516/100‑1 (MQ=255)
                                                  cTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGc                                              <  1:1185719/100‑1 (MQ=255)
                                                    aTCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCtt                                            >  2:864652/1‑100 (MQ=255)
                                                    aTCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCtt                                            >  2:2152566/1‑100 (MQ=255)
                                                      cGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCt                                          <  1:1595146/100‑1 (MQ=255)
                                                      cGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCt                                          <  2:1066091/100‑1 (MQ=255)
                                                       gAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTg                                         >  1:1563258/1‑100 (MQ=255)
                                                       gAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTg                                         >  1:899515/1‑100 (MQ=255)
                                                         aGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCa                                       <  2:373759/100‑1 (MQ=255)
                                                         aGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCa                                       <  2:2274850/100‑1 (MQ=255)
                                                         aGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCa                                       >  1:2186572/1‑100 (MQ=255)
                                                         aGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCa                                       <  1:1831041/100‑1 (MQ=255)
                                                          gCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCaa                                      <  1:2054332/100‑1 (MQ=255)
                                                           cGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCaaa                                     >  2:2011034/1‑100 (MQ=255)
                                                           cGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCaaa                                     >  1:396737/1‑100 (MQ=255)
                                                            gAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAaac                                    <  2:62094/100‑1 (MQ=255)
                                                            gAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAaac                                    >  1:2365128/1‑100 (MQ=255)
                                                            gAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAaac                                    >  1:2539744/1‑100 (MQ=255)
                                                             aaaaTGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAaaca                                   <  1:2233673/100‑1 (MQ=255)
                                                              aaaTGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAaacaa                                  <  1:1451626/100‑1 (MQ=255)
                                                                aTGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAAcg                                <  2:1830725/100‑1 (MQ=255)
                                                                aTGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAAcg                                <  1:2479851/100‑1 (MQ=255)
                                                                aTGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAAcg                                <  2:850403/100‑1 (MQ=255)
                                                                 tGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAAcgc                               >  1:1470436/1‑100 (MQ=255)
                                                                 tGCAGGAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAAcgc                               >  1:1990097/1‑100 (MQ=255)
                                                                   cacgAACTCGCACCGGTTTCCCAGAAACGGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGa                             <  2:2163068/97‑1 (MQ=255)
                                                                     ggAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGaaa                           >  1:812835/1‑100 (MQ=255)
                                                                      gAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAg                          <  2:1684374/100‑1 (MQ=255)
                                                                      gAACTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAg                          <  2:1369142/100‑1 (MQ=255)
                                                                         cTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATg                       >  1:2215221/1‑100 (MQ=255)
                                                                         cTGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATg                       <  1:681924/100‑1 (MQ=255)
                                                                          tGGCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGa                      >  2:915601/1‑100 (MQ=255)
                                                                           ggCACAGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGGGCTGATGCTTCTGCAAACAACGCGAAAGATGAc                     >  2:542326/1‑100 (MQ=255)
                                                                             cacaGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAATATGCCCAGCAGCAGGCTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGa                   >  1:1279297/1‑100 (MQ=255)
                                                                             cacaGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGa                   >  2:1841393/1‑100 (MQ=255)
                                                                              acaGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGAt                  >  2:124083/1‑100 (MQ=255)
                                                                              acaGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGAt                  >  1:2017704/1‑100 (MQ=255)
                                                                                aGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGAtgt                <  1:1323668/100‑1 (MQ=255)
                                                                                aGGTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGAtgt                <  1:1216616/100‑1 (MQ=255)
                                                                                 ggTTTCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGAtgtt               <  1:847986/100‑1 (MQ=255)
                                                                                    ttCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTc            <  2:491907/100‑1 (MQ=255)
                                                                                    ttCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTc            >  1:655467/1‑100 (MQ=255)
                                                                                     tCCCAGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCg           <  2:2251315/100‑1 (MQ=255)
                                                                                       ccAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGAc         <  1:879584/100‑1 (MQ=255)
                                                                                         aGAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGc       <  1:1056683/100‑1 (MQ=255)
                                                                                          gAAACTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCt      >  2:1043906/1‑100 (MQ=255)
                                                                                              cTGTTGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAAt  >  1:2247377/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAAT  >  NC_000913/13860‑14053

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