Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,503,169 (G)5→4 coding (263/1329 nt) insG ← IS4 transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,503,1650G.98.0% 431.5 / 1.3 99coding (267/1329 nt)insGIS4 transposase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (49/48);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 106 total aligned reads displayed.

ATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAA  >  NC_000913/4503069‑4503259
                                                                                                |                                                                                              
aTGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgggg                                                                                             <  1:2281104/100‑1 (MQ=255)
aTGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgggg                                                                                             <  2:487659/100‑1 (MQ=255)
 tGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggc                                                                                            >  2:517637/1‑99 (MQ=255)
  gcggcgTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggca                                                                                           >  2:2430347/1‑98 (MQ=255)
   cggcgTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggcaa                                                                                          >  1:356596/1‑97 (MQ=255)
    ggcgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggcaac                                                                                         <  1:2548965/100‑5 (MQ=255)
      cgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGa                                                                                      >  1:39190/1‑100 (MQ=255)
            gtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGc                                                                                <  1:1528326/100‑1 (MQ=255)
             ttATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCg                                                                               <  1:1751245/100‑1 (MQ=255)
                tGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGAtt                                                                            >  1:495812/1‑100 (MQ=255)
                 gCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTg                                                                           >  2:2800843/1‑100 (MQ=255)
                  ccAGAGCTGCGCTGTTTTCGTGAACACGTGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGc                                                                          >  1:1530542/1‑100 (MQ=255)
                  ccAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGc                                                                          >  1:1530488/1‑100 (MQ=255)
                  ccAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGc                                                                          >  2:1079209/1‑100 (MQ=255)
                   cAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGcc                                                                         >  1:108342/1‑100 (MQ=255)
                    agagCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGccc                                                                        <  2:2843584/100‑1 (MQ=255)
                    agagCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGccc                                                                        <  1:2424280/100‑1 (MQ=255)
                       gCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGgc                                                                     >  1:1963781/1‑100 (MQ=255)
                        cTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGgca                                                                    <  1:619537/100‑1 (MQ=255)
                          gcgcTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGgcagc                                                                  >  1:1654326/1‑100 (MQ=255)
                          gcgcTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGgcagc                                                                  >  2:1632869/1‑100 (MQ=255)
                          gcgcTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGgcagc                                                                  >  2:1777571/1‑100 (MQ=255)
                            gcTGTTTTCGTGAACACGCGGCGGACAGCCTCAGTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGAGCAACGAAGGGGCGATTGCCCGGCAGCat                                                                >  2:1745691/1‑100 (MQ=255)
                             cTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCatg                                                               >  1:346826/1‑100 (MQ=255)
                             cTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGCCAGCatg                                                               >  2:71110/1‑100 (MQ=255)
                              tGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGGATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCatga                                                              <  2:2450299/100‑1 (MQ=255)
                              tGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCatga                                                              <  1:2630132/100‑1 (MQ=255)
                               gTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCatgat                                                             >  1:906558/1‑100 (MQ=255)
                               gTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCatgat                                                             >  2:2655452/1‑100 (MQ=255)
                               gTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCatgat                                                             >  1:1898609/1‑100 (MQ=255)
                                ttttCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCatgatg                                                            <  1:2666600/100‑1 (MQ=255)
                                 tttCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGt                                                           >  1:373369/1‑100 (MQ=255)
                                   tCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTcc                                                         <  2:1825708/100‑1 (MQ=255)
                                     gTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAg                                                       <  1:2017518/100‑1 (MQ=255)
                                       gAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGc                                                     >  2:379635/1‑100 (MQ=255)
                                       gAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGc                                                     <  2:101513/100‑1 (MQ=255)
                                         acacGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGa                                                   >  2:2120387/1‑100 (MQ=255)
                                         acacGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGa                                                   <  1:1444704/100‑1 (MQ=255)
                                         acacGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGa                                                   >  2:2456906/1‑100 (MQ=255)
                                          cacGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGAt                                                  <  1:1401439/100‑1 (MQ=255)
                                          cacGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGAt                                                  >  1:34552/1‑100 (MQ=255)
                                          cacGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGAt                                                  >  1:786467/1‑100 (MQ=255)
                                           acGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGAtt                                                 <  1:1964720/100‑1 (MQ=255)
                                             gcggcgGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTca                                               <  2:1654313/100‑1 (MQ=255)
                                              cggcggACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTcac                                              <  1:2430363/100‑1 (MQ=255)
                                              cggcggACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTcac                                              >  1:1589891/1‑100 (MQ=255)
                                                    aCAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTg                                        <  1:2426897/100‑1 (MQ=255)
                                                      aGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATACCGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGt                                      <  1:379135/100‑1 (MQ=255)
                                                      aGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGt                                      <  1:1063910/100‑1 (MQ=255)
                                                      aGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGt                                      <  1:745715/100‑1 (MQ=255)
                                                      aGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGt                                      <  2:1293666/100‑1 (MQ=255)
                                                       gCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTg                                     >  1:2580861/1‑100 (MQ=255)
                                                       gCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTg                                     >  1:1680921/1‑100 (MQ=255)
                                                       gCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTg                                     >  1:1067515/1‑100 (MQ=255)
                                                       gCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTg                                     <  1:1484093/100‑1 (MQ=255)
                                                        ccTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGa                                    <  2:2732961/100‑1 (MQ=255)
                                                        ccTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGa                                    >  2:1382328/1‑100 (MQ=255)
                                                         cTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGaa                                   >  1:1233826/1‑100 (MQ=255)
                                                          tCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAag                                  >  1:2878466/1‑100 (MQ=255)
                                                          tCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAag                                  <  1:2644464/100‑1 (MQ=255)
                                                           cACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAaga                                 <  2:108342/100‑1 (MQ=255)
                                                           cACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAaga                                 <  2:184647/100‑1 (MQ=255)
                                                            aCTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAagag                                <  2:2382524/100‑1 (MQ=255)
                                                            aCTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAagag                                >  2:2829226/1‑100 (MQ=255)
                                                             cTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAgg                               >  2:32530/1‑100 (MQ=255)
                                                              ttCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGt                              <  2:1234559/100‑1 (MQ=255)
                                                               tCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGtt                             >  1:1443609/1‑100 (MQ=255)
                                                               acccAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAgg                               <  1:2102467/98‑1 (MQ=255)
                                                                  cAGGCGCTGGCGGGCCTGAATAACGGCACT‑‑GGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGGTTCACAATCTGGTGAAGAGGTTCttt                         >  2:1544009/1‑100 (MQ=255)
                                                                   aGGCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCttt                         <  2:2795560/100‑1 (MQ=255)
                                                                    ggCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCtt                          <  2:600341/98‑1 (MQ=255)
                                                                    ggCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTa                        <  2:34552/100‑1 (MQ=255)
                                                                    ggCGCTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTa                        <  1:517639/100‑1 (MQ=255)
                                                                      cgcTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGTCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTAc                       >  1:292163/1‑100 (MQ=255)
                                                                       gcTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGc                     <  2:1606041/100‑1 (MQ=255)
                                                                        cTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCt                    <  1:1114403/100‑1 (MQ=255)
                                                                        cTGGCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCt                    <  2:2103664/100‑1 (MQ=255)
                                                                          ggCGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCa                  <  2:65813/100‑1 (MQ=255)
                                                                            cGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAg                >  1:287975/1‑100 (MQ=255)
                                                                            cGGGCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAg                >  1:223547/1‑100 (MQ=255)
                                                                               gCCTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAgcgc             >  1:1929594/1‑100 (MQ=255)
                                                                                ccTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGTAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGcc            >  2:1891470/1‑100 (MQ=255)
                                                                                ccTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGcc            >  1:392408/1‑100 (MQ=255)
                                                                                ccTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGcc            <  2:2580842/100‑1 (MQ=255)
                                                                                 cTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCa           <  1:2196218/100‑1 (MQ=255)
                                                                                 cTGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCa           <  2:1446112/100‑1 (MQ=255)
                                                                                  tGAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCAt          >  2:1106487/1‑100 (MQ=255)
                                                                                   gAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATg         <  2:1067513/100‑1 (MQ=255)
                                                                                   gAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATg         <  1:2163801/100‑1 (MQ=255)
                                                                                   gAATAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATg         >  1:103915/1‑100 (MQ=255)
                                                                                    aaTAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGc        <  2:2632339/100‑1 (MQ=255)
                                                                                     aTAACGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGcc       >  1:2850047/1‑100 (MQ=255)
                                                                                      tagcGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGcc       <  1:2045757/96‑1 (MQ=255)
                                                                                        aCGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaac    <  2:1421820/100‑1 (MQ=255)
                                                                                        aCGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaac    >  1:2210053/1‑100 (MQ=255)
                                                                                         cGGCACTGGGTCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGTAGTGGTTCTTTACGCTCAAGCGCCATGCcaaca   <  2:1898679/100‑1 (MQ=255)
                                                                                         cGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaaca   >  1:1977579/1‑100 (MQ=255)
                                                                                         cGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaaca   >  1:2497562/1‑100 (MQ=255)
                                                                                         cGGCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaac    <  2:2878445/99‑1 (MQ=255)
                                                                                          ggCACTGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCcaacaa  <  1:938885/100‑1 (MQ=255)
                                                                                                |                                                                                              
ATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCATGATGTCCAGGCGATTCACAATCTGGTGAAGAGGTTCTTTACGCTCAAGCGCCATGCCAACAA  >  NC_000913/4503069‑4503259

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