Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,694 A→G 15.4% 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.4% 282.1 / 42.3 117intergenic (+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 (51/48);  new base G (5/13);  total (56/61)
Fisher's exact test for biased strand distribution p-value = 7.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

Only 100 of 118 total aligned reads displayed.

ATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTT  >  NC_000913/681595‑681792
                                                                                                   |                                                                                                  
  cAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGc                                                                                                  <  1:821856/100‑1 (MQ=255)
   aaaaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCg                                                                                                 <  1:266740/100‑1 (MQ=255)
     aaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGtt                                                                                               >  2:1108787/1‑100 (MQ=255)
      aTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGttt                                                                                              <  1:292369/100‑1 (MQ=255)
          gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgc                                                                                          >  2:1746446/1‑100 (MQ=255)
          gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgc                                                                                          >  2:1941879/1‑100 (MQ=255)
          gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcg                                                                                           >  1:2328496/1‑99 (MQ=255)
          gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGAATGCAGCGTTTAcgc                                                                                          >  1:3064343/1‑100 (MQ=255)
           tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                         <  2:434170/100‑1 (MQ=255)
            ttAACGGTGCTGGTGATAGGTTATGCTCTCGGTAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgccg                                                                                        >  2:1478069/1‑100 (MQ=255)
            ttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgccg                                                                                        <  2:1888922/100‑1 (MQ=255)
            ttAACGGTGCTGGGGATAGGGTATGCTCACGGCTAACTCGCGATGTTGTTTAGCTGACGAGGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                         <  2:2904784/99‑1 (MQ=255)
             tAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGtaaacgctgc                                                                                       >  2:2058685/1‑91 (MQ=255)
              aaCGGGGCTGGGGATAGGTTAGGCTCGCGGCAAAATCGTGATGTTGTGTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                         <  2:540177/97‑1 (MQ=255)
                 ggTGCTGGTGATAGGTTATGCTCTCGGCAATATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATcc                                                                                   <  2:988098/100‑1 (MQ=38)
                 ggTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATcc                                                                                   <  1:45478/100‑1 (MQ=38)
                  gTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCg                                                                                  <  1:2035848/100‑1 (MQ=38)
                  gTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCg                                                                                  <  2:2630673/100‑1 (MQ=255)
                  gTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCg                                                                                  <  1:2521885/100‑1 (MQ=255)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAATCGCTGCATCCgg                                                                                 <  1:2046390/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                 <  1:1879462/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                 <  1:1975813/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                 <  1:1799329/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCgg                                                                                 >  2:632971/1‑100 (MQ=255)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCgg                                                                                 >  1:657180/1‑100 (MQ=255)
                    gCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGc                                                                                <  2:2476384/100‑1 (MQ=255)
                     cTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGCa                                                                               <  1:1928913/100‑1 (MQ=38)
                      tGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAg                                                                              <  1:1152789/100‑1 (MQ=255)
                         tGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCg                                                                           >  2:11125/1‑100 (MQ=255)
                          gATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgt                                                                          <  2:3074350/100‑1 (MQ=255)
                          gATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCGGAGGAGGGGTTTATGGGGCGGCTGCCGTCTGCAGCGTTTgcgcgtcgtccg                                                                                  >  2:2944992/1‑84 (MQ=255)
                           aTAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtg                                                                         <  2:1623470/100‑1 (MQ=255)
                           aTAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtg                                                                         <  1:2860174/100‑1 (MQ=255)
                            tAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtgt                                                                        <  1:632972/100‑1 (MQ=255)
                            tAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtgt                                                                        <  2:2659997/100‑1 (MQ=255)
                              ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGttt                                                                      >  1:2421190/1‑100 (MQ=255)
                              ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGttt                                                                      >  1:156628/1‑100 (MQ=255)
                               gTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTa                                                                     <  1:926898/100‑1 (MQ=255)
                               gTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTa                                                                     <  2:852065/100‑1 (MQ=255)
                                ttATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTAt                                                                    <  1:1161349/100‑1 (MQ=255)
                                  aTGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATgg                                                                  <  1:607127/100‑1 (MQ=255)
                                   tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                 <  2:999272/100‑1 (MQ=255)
                                   tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCACCCGTGTTTATgtg                                                                 <  1:2111676/100‑3 (MQ=255)
                                   cgatctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGCTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                 <  1:2090785/97‑1 (MQ=255)
                                    gCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGa                                                                >  1:401716/1‑100 (MQ=255)
                                     ctctCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAt                                                               >  1:2890148/1‑100 (MQ=255)
                                      tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtc                                                              <  1:1056468/100‑1 (MQ=255)
                                      tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtc                                                              <  1:1276085/100‑1 (MQ=255)
                                         cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTGATGGGAtctc                                                            >  1:1470981/1‑99 (MQ=255)
                                         cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGATTACGCCGCATCCGGCAGCCGTGTTTATGGGagatcg                                                           >  2:2090768/1‑95 (MQ=255)
                                          ggCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCaa                                                          <  1:638001/100‑1 (MQ=255)
                                          ggCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCAGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCaa                                                          >  1:923279/1‑100 (MQ=255)
                                           gCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                         <  1:2057972/100‑1 (MQ=255)
                                           gCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                         >  1:363708/1‑100 (MQ=255)
                                           gCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                         >  1:259916/1‑100 (MQ=255)
                                            cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                        >  1:1269457/1‑100 (MQ=255)
                                            cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                        >  1:1566154/1‑100 (MQ=255)
                                            cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                        >  1:1057086/1‑100 (MQ=255)
                                             aaaaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAGTgg                                                       <  2:440175/100‑1 (MQ=255)
                                             aaaaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATgg                                                       <  1:2083905/100‑1 (MQ=255)
                                               aaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTTTTTTTTTTTTGGATGCAGCGCTTACGCCGCATCCGGCATCCGTGTTTATGGGATCTCAATGGCt                                                     >  2:622340/1‑100 (MQ=255)
                                                       tgttgtTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                             <  2:2342512/100‑1 (MQ=255)
                                                        gttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcg                                            <  2:2074073/100‑1 (MQ=255)
                                                             ttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcgc                                       >  1:136309/1‑100 (MQ=255)
                                                             ttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcgc                                       <  2:1761203/100‑1 (MQ=255)
                                                               aGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgccc                                     <  1:1108792/100‑1 (MQ=255)
                                                                gCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccg                                    >  1:115175/1‑100 (MQ=255)
                                                                 cTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgc                                   >  1:1043303/1‑100 (MQ=255)
                                                                 cTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgc                                   >  2:526980/1‑100 (MQ=255)
                                                                 cTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgc                                   >  1:1917763/1‑100 (MQ=255)
                                                                  tGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgcc                                  >  2:546159/1‑100 (MQ=255)
                                                                     cGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgcc                               >  2:903448/1‑100 (MQ=255)
                                                                          tgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGCTCTCAATGGCAAAATAATCGCTCGCCCGCCGCCGAACt                          >  1:2358159/1‑100 (MQ=255)
                                                                           gtTTTTTGCGCTGTTGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAACTAATCGCTCTCCCGCCGCCGAACtt                         >  2:1982886/1‑100 (MQ=39)
                                                                           gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                         >  1:2937249/1‑100 (MQ=255)
                                                                           gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                         >  1:2513365/1‑100 (MQ=255)
                                                                           gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                         >  2:1477318/1‑100 (MQ=255)
                                                                                ttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATcc                    <  1:2068659/100‑1 (MQ=255)
                                                                                   cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcg                 >  1:2953645/1‑100 (MQ=255)
                                                                                   cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcg                 >  1:2889080/1‑100 (MQ=255)
                                                                                   cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcg                 >  2:3040297/1‑100 (MQ=255)
                                                                                   cacggctGCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGAATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcg                 >  2:2182239/7‑100 (MQ=21)
                                                                                    gcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcgc                >  1:2172729/1‑100 (MQ=255)
                                                                                    gcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAc                  >  2:901100/1‑98 (MQ=255)
                                                                                     ccGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcgcg               >  2:2436011/1‑100 (MQ=255)
                                                                                     ccGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcgcg               >  1:2842647/1‑100 (MQ=255)
                                                                                      cGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCgg              <  2:259914/100‑1 (MQ=255)
                                                                                       gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg             <  2:2264829/100‑1 (MQ=255)
                                                                                       gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg             <  2:423142/100‑1 (MQ=255)
                                                                                       gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg             <  2:622360/100‑1 (MQ=255)
                                                                                       gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg             <  1:1462870/100‑1 (MQ=255)
                                                                                          gCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa          >  1:2560544/1‑100 (MQ=255)
                                                                                           tcGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat         <  2:2659507/99‑1 (MQ=255)
                                                                                           ccGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat         >  1:360634/1‑100 (MQ=255)
                                                                                            cGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCata        >  1:1677181/1‑100 (MQ=38)
                                                                                            cGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCata        >  1:2659535/1‑100 (MQ=255)
                                                                                            cGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTAATCCACGCGGGCCata        >  2:1820779/1‑100 (MQ=255)
                                                                                              gATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATACTCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATc      >  2:1827412/1‑100 (MQ=38)
                                                                                              gATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCatat       >  2:1465650/1‑99 (MQ=38)
                                                                                                  cAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGttt  <  2:1057081/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
ATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTT  >  NC_000913/681595‑681792

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