Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,694 A→G 18.0% intergenic (+31/+29) djlC → / ← hscC J domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913681,6940AG18.0% 242.5 / 48.4 118intergenic (+31/+29)djlC/hscCJ domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor
Reads supporting (aligned to +/- strand):  ref base A (40/56);  new base G (11/10);  total (51/67)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-03

Only 100 of 124 total aligned reads displayed.

ATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTG  >  NC_000913/681595‑681793
                                                                                                   |                                                                                                   
aTCAAAATCAGTTTAACGGTGCTGGTTATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCa                                                                                                     <  2:509103/100‑1 (MQ=255)
 tCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCTGATGCAg                                                                                                    <  1:1583059/100‑1 (MQ=255)
 tCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAg                                                                                                    <  1:35606/100‑1 (MQ=255)
  cAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGc                                                                                                   <  1:1527002/100‑1 (MQ=255)
   aaaaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCg                                                                                                  >  1:2177924/1‑100 (MQ=255)
    aaaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATTCCGGATGCAGcct                                                                                                 <  1:2041992/100‑3 (MQ=255)
     aaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGtt                                                                                                >  2:321960/1‑100 (MQ=255)
       tCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGtaaa                                                                                              <  2:1971728/100‑4 (MQ=255)
         aGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcg                                                                                            >  2:1107210/1‑100 (MQ=255)
          gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGtt                                                                                                >  2:143629/1‑95 (MQ=255)
          gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgc                                                                                           >  1:1062831/1‑100 (MQ=255)
          gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCATTTAcgc                                                                                           >  1:1496771/1‑100 (MQ=255)
           tttAACGGTGCTGGTGATTGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACTATGTGTTTTTTGCGCCGATGCTGGATGCAGCGTTTAcgcc                                                                                          <  2:198158/100‑1 (MQ=255)
           tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                          >  1:2264943/1‑100 (MQ=255)
           tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                          >  1:2264171/1‑100 (MQ=255)
               aCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCAt                                                                                      >  1:282610/1‑100 (MQ=255)
                 ggTGCTTGTGATAGGTTATTCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACTATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATcc                                                                                    <  1:2488153/100‑1 (MQ=255)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGTGCCGATGCCGGATGCAGCGTTTACGCCGCATCCgg                                                                                  <  2:2235840/100‑1 (MQ=255)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  2:1768826/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  2:240127/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  2:58948/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  1:1841106/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  1:1812307/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCg                                                                                   <  2:1262349/99‑1 (MQ=255)
                    gCTGGTGATAGGTTATGCTCTCGGCAAATTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGc                                                                                 <  2:766804/100‑1 (MQ=255)
                    gCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGc                                                                                 <  1:78333/100‑1 (MQ=38)
                    gCTGGTGATAGGTGATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGTATGCGGCGTAAACTCTGCATCCgg                                                                                  <  2:466901/99‑1 (MQ=255)
                     cTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGCa                                                                                <  1:1890559/100‑1 (MQ=38)
                     cTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCa                                                                                <  2:1040266/100‑1 (MQ=255)
                      tGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGc                                                                                 >  2:2142189/1‑98 (MQ=255)
                       ggTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCGGGCAg                                                                               <  2:1305924/99‑1 (MQ=255)
                        gTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCTATGCCGTATGCAGCGTTTACGCAGCATCCGGCATcc                                                                             >  1:536845/1‑100 (MQ=255)
                        gTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGcc                                                                             <  1:391548/100‑1 (MQ=255)
                        gTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGcc                                                                             >  1:1219539/1‑100 (MQ=255)
                          gATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgt                                                                           >  2:1905528/1‑100 (MQ=255)
                             aGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACTCCGCATCCGGCAGCCGTGtt                                                                        <  2:1182152/100‑1 (MQ=255)
                              ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGttt                                                                       <  1:1577821/100‑1 (MQ=255)
                               gTTAGGCTCTCGGCAAAATCGCGATTTTGTTTAGCTTTCGATGTGTGTTTGGCGCCTATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTa                                                                      <  1:1765048/100‑1 (MQ=38)
                                 tATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATg                                                                    <  2:2182723/100‑1 (MQ=255)
                                 tATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATg                                                                    <  1:1709176/100‑1 (MQ=255)
                                   tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                  <  1:835727/100‑1 (MQ=255)
                                   tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                  <  1:2172997/100‑1 (MQ=255)
                                   tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                  <  2:1355921/100‑1 (MQ=255)
                                    gCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGa                                                                 <  2:148919/100‑1 (MQ=255)
                                     ctctCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAt                                                                >  1:2913601/1‑100 (MQ=255)
                                     ctctCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAt                                                                >  1:1919778/1‑100 (MQ=255)
                                     ctctCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCAGTGTTTATGGGa                                                                 >  1:1337053/1‑99 (MQ=255)
                                      tctcGGCAAAATCGCGATGTTGTTTAGCTGACGTTGTGTTTTTTGTGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtc                                                               <  1:1194825/100‑1 (MQ=255)
                                        tcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtc                                                               <  1:2180998/98‑1 (MQ=255)
                                         cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCa                                                            >  2:726633/1‑100 (MQ=255)
                                           gCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                          >  1:2302919/1‑100 (MQ=255)
                                            cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                         <  1:910773/100‑1 (MQ=255)
                                            cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGACAGCCGTGTTTATGGGATCTCAAt                                                          >  1:2819836/1‑99 (MQ=255)
                                             aaaaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATgg                                                        >  2:57055/1‑100 (MQ=255)
                                                aTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTa                                                     <  2:1452034/100‑1 (MQ=255)
                                                 tcgttaTGTTTTTTATTTGACTATTTGTTTTTTGCGCCGAGTCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGc                                                       <  1:1341139/92‑1 (MQ=37)
                                                 tCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTATACGCAGCATCCGGCAGCCGTGTTTATTGAAACTCAATGGCTa                                                     >  1:502954/1‑99 (MQ=255)
                                                    cgATGGTGTTTCGGTGACGATGTGTTTTTTGCGCCGATGGCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaat                                                  <  2:625774/99‑1 (MQ=39)
                                                      aTGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataat                                               <  1:1905551/100‑1 (MQ=255)
                                                      aTGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACACCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataat                                               >  1:2421653/1‑100 (MQ=255)
                                                       tgttgtTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGTGATCTCAATGGCTAAATAAtc                                              <  1:1265807/100‑1 (MQ=255)
                                                       tgttgtTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                              <  1:2506353/100‑1 (MQ=255)
                                                        gttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcg                                             <  1:692760/100‑1 (MQ=255)
                                                         ttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgt                                            <  2:1358167/100‑2 (MQ=255)
                                                           gttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctc                                          >  1:1227829/1‑100 (MQ=255)
                                                            tttAGCTGGCGATGTCTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCTGTGTTTATGGGATCTGTATGGCTAAATAAtcg                                             <  2:2389790/96‑1 (MQ=38)
                                                             ttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcgc                                        <  2:2512453/100‑1 (MQ=255)
                                                                gCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccg                                     >  2:407455/1‑100 (MQ=255)
                                                                gCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccg                                     >  1:430496/1‑100 (MQ=255)
                                                                  tGACGATTTTTTTTTTGCGTCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgcc                                   <  1:1323041/100‑1 (MQ=255)
                                                                  tGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgcc                                   >  1:1084891/1‑100 (MQ=255)
                                                                  tGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgcc                                   >  1:269591/1‑100 (MQ=255)
                                                                     cGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgcc                                >  1:799391/1‑100 (MQ=255)
                                                                     cGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGAATCTCAATGGCTAAATAATCGCTCGCccgccgcc                                >  1:1001521/1‑100 (MQ=255)
                                                                      gATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgccg                               <  2:879141/100‑1 (MQ=255)
                                                                      gATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgccg                               >  2:1052042/1‑100 (MQ=255)
                                                                       aTGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGa                              <  1:1100202/100‑1 (MQ=255)
                                                                        tgtgTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGaa                             >  1:944038/1‑100 (MQ=255)
                                                                         gtgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAc                            <  2:510252/100‑1 (MQ=255)
                                                                         gtgtTTGTTGGGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAc                            <  2:2853121/100‑1 (MQ=255)
                                                                           gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                          >  1:958609/1‑100 (MQ=255)
                                                                             tttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCa                        <  2:2245906/100‑1 (MQ=255)
                                                                              ttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCAt                       <  1:634829/100‑1 (MQ=255)
                                                                                   cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTCTGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAc                   >  2:445747/1‑99 (MQ=255)
                                                                                   cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcg                  >  2:758354/1‑100 (MQ=255)
                                                                                       gctGCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGACTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg              >  1:382132/3‑100 (MQ=37)
                                                                                       gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg              <  1:1704667/100‑1 (MQ=255)
                                                                                       gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCgg               <  1:2170182/99‑1 (MQ=255)
                                                                                       gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGAATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg              >  1:2829638/1‑100 (MQ=255)
                                                                                          gCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa           >  2:2661359/1‑100 (MQ=38)
                                                                                          gCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa           <  2:2962190/100‑1 (MQ=255)
                                                                                           ccGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat          >  1:2670222/1‑100 (MQ=38)
                                                                                           ccGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCACTCAACCGCAGCCGAACTTCATCCACGCGGGCCat          >  1:857975/1‑100 (MQ=25)
                                                                                            cGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCata         >  1:936230/1‑100 (MQ=38)
                                                                                            cGGATGCGGCGTAAACGCTGCATCCAGCAGCCGTGTTTATGGGATCGCAATGGCTAAATCATCGCTCGCCCGCCGCAGAACTTCATCCa                    >  2:1580197/1‑89 (MQ=25)
                                                                                             ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCatat        >  1:1135278/1‑100 (MQ=38)
                                                                                                   aGCGTTTACGCCGCATCCGGCAGTCGTGTTTATGGGATTTTAATGGCTAAATAATCGCTCGCGCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  <  1:1337250/100‑1 (MQ=39)
                                                                                                   aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  <  2:2264146/100‑1 (MQ=255)
                                                                                                   aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  <  1:438793/100‑1 (MQ=255)
                                                                                                   aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  <  2:944032/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTG  >  NC_000913/681595‑681793

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