Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,694 A→G 15.1% 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,6940AG15.1% 389.6 / 59.3 169intergenic (+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 (73/69);  new base G (14/11);  total (87/82)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-02

Only 100 of 174 total aligned reads displayed.

TCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTG  >  NC_000913/681596‑681793
                                                                                                  |                                                                                                   
 cAAAATCAGTCTGGCGGTGCTGGTGATAGCTTATGCTCTCGGCACCATCTCCATGTTGTTTAGCTGACGATGTGTGTTTTGCGCCGATGCCGGATGCAg                                                                                                    <  1:1995091/99‑1 (MQ=255)
  aaaaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCg                                                                                                  <  2:660267/100‑1 (MQ=255)
   aaaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGt                                                                                                 <  2:2856522/100‑1 (MQ=255)
      tCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTa                                                                                              <  2:59025/100‑1 (MQ=255)
         gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgc                                                                                           >  2:3496596/1‑100 (MQ=255)
         gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgc                                                                                           >  2:1296817/1‑100 (MQ=255)
         gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgc                                                                                           >  2:4120055/1‑100 (MQ=255)
           ttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgccg                                                                                         <  2:18482/100‑1 (MQ=255)
            tAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgccg                                                                                         <  2:578815/99‑1 (MQ=255)
             aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCa                                                                                       >  2:753686/1‑100 (MQ=255)
             aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCa                                                                                       >  2:957565/1‑100 (MQ=255)
               cGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGTATGCGGCGTAAACGCTGCATc                                                                                     <  2:751829/100‑1 (MQ=255)
               cGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATc                                                                                     >  2:1372821/1‑100 (MQ=255)
               cGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATc                                                                                     >  1:4356589/1‑100 (MQ=255)
                 gTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCg                                                                                   <  2:430612/100‑1 (MQ=38)
                  tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  1:1530914/100‑1 (MQ=38)
                   gCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  1:133989/99‑1 (MQ=38)
                   gCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGc                                                                                 <  2:2750804/100‑1 (MQ=38)
                   gCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTGTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGc                                                                                 <  2:683962/100‑1 (MQ=38)
                      ggttaTAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGCatc                                                                              <  2:3201712/96‑3 (MQ=25)
                      ggTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGc                                                                              >  1:996673/1‑100 (MQ=255)
                        tGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCg                                                                            >  2:4364330/1‑100 (MQ=255)
                         gATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgt                                                                           >  2:3476287/1‑100 (MQ=255)
                           tAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtgt                                                                         >  1:3231499/1‑100 (MQ=255)
                           tAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCg                                                                            >  2:2685410/1‑97 (MQ=255)
                             ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTTCTCCTTTGCCGGATGCAGCGTTTACGCCCCATCCGGCAGCCGTGttt                                                                       >  1:1220195/1‑100 (MQ=255)
                             ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGttt                                                                       <  2:595239/100‑1 (MQ=255)
                             ggTTATGCTCTCGGCAAAATCGCGATGTGGTTTGGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGACGCCGTGttt                                                                       >  1:2884446/1‑100 (MQ=255)
                                tATGCTCTCGGCAAAATCGCGATGTTGTTTAGGTGCCGAGGTCTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATg                                                                    <  1:3084366/100‑1 (MQ=39)
                                tATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATg                                                                    <  2:4002771/100‑1 (MQ=255)
                                   gCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGa                                                                 <  1:3756705/100‑1 (MQ=255)
                                   gCTCTCGGCAAAATCGCGATGTTGTTCAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGa                                                                 <  2:3766258/100‑1 (MQ=255)
                                    atctcGGCAAAATCGCGATGTTGTTTAGCTGCCGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                  >  1:562854/2‑98 (MQ=255)
                                     tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtc                                                               <  1:886826/100‑1 (MQ=255)
                                      ctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtct                                                              >  2:4266946/1‑100 (MQ=255)
                                       tcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtctc                                                             >  1:618508/1‑100 (MQ=255)
                                        cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCa                                                            >  2:1190201/1‑100 (MQ=255)
                                        cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCa                                                            >  2:1457324/1‑100 (MQ=255)
                                        cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCa                                                            >  2:1945310/1‑100 (MQ=255)
                                        cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATATCa                                                            >  1:905285/1‑100 (MQ=255)
                                          gCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                          >  2:3532237/1‑100 (MQ=255)
                                          gCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                          <  1:2612253/100‑1 (MQ=255)
                                           cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                          >  1:143180/1‑99 (MQ=255)
                                           cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                         <  1:4341012/100‑1 (MQ=255)
                                             aaaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGc                                                       <  2:2113335/100‑1 (MQ=255)
                                             aaaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGc                                                       <  2:105957/100‑1 (MQ=255)
                                             aaaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCAGCATTCGGCAGCCGTGTTTATGGGATCTCaa                                                           >  1:1038375/1‑96 (MQ=255)
                                             aaaTCGCGATGTTGTTGAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGc                                                       <  1:4300773/100‑1 (MQ=255)
                                                 cgcgATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCAGGCAGCCGTGTTTATGGGATCTCAATGGCTaaa                                                   <  1:663167/100‑1 (MQ=255)
                                                    gATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataa                                                >  1:1956187/1‑100 (MQ=255)
                                                       gttgttTAGCTGACGGTGTGTTTTTGGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcg                                             <  2:1327344/100‑1 (MQ=255)
                                                       gttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcg                                             <  1:115831/100‑1 (MQ=255)
                                                        ttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgc                                            <  2:1678122/100‑1 (MQ=255)
                                                           tttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcg                                         <  1:2479919/100‑1 (MQ=255)
                                                            ttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcgc                                        <  1:474063/100‑1 (MQ=255)
                                                              aGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgccc                                      >  2:1632870/1‑100 (MQ=255)
                                                               gCTGACGATGTGTTTTTTTCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccg                                     <  1:3541769/100‑1 (MQ=255)
                                                               gCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTGATGGGATCTCAATGGCTAAATAATCGCTcgcccg                                     >  1:2319296/1‑100 (MQ=255)
                                                                cTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgc                                    >  2:3301440/1‑100 (MQ=255)
                                                                  gACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccg                                  <  1:539677/100‑1 (MQ=255)
                                                                        gtgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAc                            <  1:2613581/100‑1 (MQ=255)
                                                                        gtgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAc                            <  2:814117/100‑1 (MQ=255)
                                                                        gtgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAc                            <  1:439223/100‑1 (MQ=255)
                                                                        gtgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCCCCGCCGAAc                            >  1:1252206/1‑100 (MQ=255)
                                                                          gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGGCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                          >  2:4169536/1‑100 (MQ=255)
                                                                           ttttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTc                         <  1:2123419/100‑1 (MQ=255)
                                                                            tttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCa                        <  2:2465002/100‑1 (MQ=255)
                                                                            tttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCa                        <  2:3185499/100‑1 (MQ=255)
                                                                              tttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATc                      <  2:263554/100‑1 (MQ=255)
                                                                              tttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATc                      <  1:821887/100‑1 (MQ=255)
                                                                                tGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCa                    >  2:3838573/1‑100 (MQ=255)
                                                                                 gcgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAc                   >  2:3648444/1‑100 (MQ=255)
                                                                                 gcgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCACCGAACTTCaagatc                   >  2:4343702/1‑95 (MQ=255)
                                                                                 acacggctGCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCACTGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAc                   >  2:3238987/8‑100 (MQ=21)
                                                                                  cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcg                  <  1:4395797/100‑1 (MQ=255)
                                                                                    ccGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcgcg                >  2:154377/1‑100 (MQ=255)
                                                                                      gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg              >  1:13898/1‑100 (MQ=255)
                                                                                       ctGCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGc             >  2:2071294/2‑100 (MQ=37)
                                                                                        tGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGcc            <  2:1767503/100‑1 (MQ=255)
                                                                                         gCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa           >  2:3329462/1‑100 (MQ=38)
                                                                                         gCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa           <  2:1381867/100‑1 (MQ=255)
                                                                                          ccGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat          >  2:524929/1‑100 (MQ=38)
                                                                                          ccGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat          >  2:3230612/1‑100 (MQ=255)
                                                                                           cGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCata         >  1:3433523/1‑100 (MQ=38)
                                                                                           cGGATGCAGCGTTTACGCCGCATCCGACAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCata         >  2:2675969/1‑100 (MQ=255)
                                                                                            ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCatat        >  1:1422276/1‑100 (MQ=38)
                                                                                            ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCatat        >  1:4218452/1‑100 (MQ=38)
                                                                                            ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCatat        >  2:1126303/1‑100 (MQ=38)
                                                                                            ggATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCatat        <  2:3780531/100‑1 (MQ=255)
                                                                                             gATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGGTCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTGCATCCACGCGGGCCATATc       >  1:2394079/1‑100 (MQ=38)
                                                                                              aTGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCg      <  2:1242068/100‑1 (MQ=255)
                                                                                               tGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGt     >  2:1473425/1‑100 (MQ=255)
                                                                                                gCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGtt    <  1:3225391/100‑1 (MQ=255)
                                                                                                 cAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGttt   >  1:173468/1‑100 (MQ=255)
                                                                                                 cAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGttt   >  1:1034540/1‑100 (MQ=255)
                                                                                                 cAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGtt    >  1:3178279/1‑99 (MQ=255)
                                                                                                  aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  <  1:2842610/100‑1 (MQ=255)
                                                                                                  aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  >  1:356006/1‑100 (MQ=255)
                                                                                                  aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  >  2:3106163/1‑100 (MQ=255)
                                                                                                  aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  <  1:4004810/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTG  >  NC_000913/681596‑681793

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: