Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,650,1470GA20.5% 273.3 / 61.6 128pseudogene (70/765 nt)ydfEQin prophage; pseudogene;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (50/51);  new base A (12/14);  total (62/66)
Fisher's exact test for biased strand distribution p-value = 8.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

Only 100 of 131 total aligned reads displayed.

CCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATA  >  NC_000913/1650050‑1650246
                                                                                                 |                                                                                                   
ccTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTa                                                                                                   <  2:2106634/100‑1 (MQ=255)
 cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAg                                                                                                  >  1:784197/1‑100 (MQ=255)
 cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAg                                                                                                  <  1:1838489/100‑1 (MQ=255)
  ttAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGc                                                                                                 <  1:3231732/100‑1 (MQ=255)
  ttAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGc                                                                                                 <  1:2447619/100‑1 (MQ=255)
   tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGcc                                                                                                <  1:959067/100‑1 (MQ=255)
     aaTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCcac                                                                                              <  2:101231/100‑1 (MQ=255)
      aTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCcaca                                                                                             <  2:1376726/100‑1 (MQ=255)
        cTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGc                                                                                           >  1:648760/1‑100 (MQ=255)
        cTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGc                                                                                           <  1:879652/100‑1 (MQ=255)
        cTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGc                                                                                           <  2:1135277/100‑1 (MQ=255)
        cTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGc                                                                                           <  1:1895910/100‑1 (MQ=255)
        cTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGc                                                                                           <  2:2656591/100‑1 (MQ=255)
         tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTCTTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTTCAGTCACAGcg                                                                                          <  2:1768224/100‑2 (MQ=255)
         tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcg                                                                                          <  1:1024163/100‑2 (MQ=255)
         tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGACAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcg                                                                                          <  1:2585342/100‑2 (MQ=255)
          gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgc                                                                                         <  1:1811060/100‑3 (MQ=255)
          gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgc                                                                                         <  1:862545/100‑3 (MQ=255)
          gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgc                                                                                         <  1:3172781/100‑3 (MQ=255)
           tCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgcc                                                                                        <  2:2745332/100‑4 (MQ=255)
            cATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCAtt                                                                                       >  1:64585/1‑100 (MQ=255)
            cATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgcct                                                                                       <  1:1441009/100‑5 (MQ=255)
             aTGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTg                                                                                      <  1:1786634/100‑1 (MQ=255)
             aTGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTg                                                                                      <  2:1110755/100‑1 (MQ=255)
             aTGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTg                                                                                      <  1:2842245/100‑1 (MQ=255)
                 ggTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGc                                                                                  >  1:2559659/1‑100 (MQ=255)
                 ggTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGc                                                                                  >  2:109593/1‑100 (MQ=255)
                 ggTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGc                                                                                  >  2:516419/1‑100 (MQ=255)
                  gTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCt                                                                                 >  2:1846943/1‑100 (MQ=255)
                    tAAAACAAATGAGTAACGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTgg                                                                               <  2:2748787/100‑1 (MQ=255)
                    tAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTgg                                                                               <  1:2904829/100‑1 (MQ=255)
                     aaaaCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGc                                                                              <  1:1106616/100‑1 (MQ=255)
                     aaaaCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGc                                                                              >  1:764851/1‑100 (MQ=255)
                      aaaCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcg                                                                             <  2:2040816/100‑1 (MQ=255)
                       aaCAAATGAGTAATGTCTTTCTTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGCAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcga                                                                            >  2:2066642/1‑100 (MQ=255)
                         cAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgacg                                                                          >  2:2958415/1‑100 (MQ=255)
                          aaaTGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgacga                                                                         >  2:393365/1‑100 (MQ=255)
                             tGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACg                                                                      >  2:3147119/1‑100 (MQ=255)
                               aGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGcc                                                                    <  2:2178268/100‑1 (MQ=255)
                                gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCg                                                                   >  1:596654/1‑100 (MQ=255)
                                gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCg                                                                   >  2:2110131/1‑100 (MQ=255)
                                        tttATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCg                                                             >  1:438337/1‑98 (MQ=255)
                                         ttATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGc                                                          <  2:273415/100‑1 (MQ=255)
                                           aTTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaa                                                        >  2:887915/1‑100 (MQ=255)
                                             ttGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaaaa                                                      <  1:1328863/100‑1 (MQ=255)
                                                 gccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTc                                                  <  2:3138128/100‑1 (MQ=255)
                                                  ccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCGGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTc                                                  <  2:814057/99‑1 (MQ=255)
                                                  ccgcTAGTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCa                                                 <  2:2545680/100‑1 (MQ=255)
                                                   cgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                                <  1:1205181/100‑1 (MQ=255)
                                                     cTATTCCTGACGAACTTGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                                <  2:2301673/98‑1 (MQ=255)
                                                     cTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTg                                              >  1:1267117/1‑100 (MQ=255)
                                                       aTTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGc                                            <  2:764849/100‑1 (MQ=255)
                                                             gACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTcc                                      >  1:451933/1‑100 (MQ=255)
                                                             gACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTcc                                      <  2:1736388/100‑1 (MQ=255)
                                                               cGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTg                                    >  1:1850113/1‑100 (MQ=255)
                                                                  aCTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGaa                                 >  1:3344400/1‑100 (MQ=255)
                                                                   cTGGTAACAAGGGAATAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGTCGACGTACGCCGTGCTCGAGGAAAATTTCCCTGGCAATTCCTGGAAc                                <  2:3001599/100‑1 (MQ=255)
                                                                   cTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAAc                                <  2:2394503/100‑1 (MQ=255)
                                                                   cTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAAc                                <  1:440390/100‑1 (MQ=255)
                                                                   cTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAAc                                <  2:3001725/100‑1 (MQ=255)
                                                                    tGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACa                               <  2:3344381/100‑1 (MQ=255)
                                                                     ggCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAt                              >  2:1441305/1‑100 (MQ=255)
                                                                     ggCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAt                              >  2:2132122/1‑100 (MQ=255)
                                                                       caacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTa                            >  1:671025/1‑100 (MQ=255)
                                                                       caacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTa                            >  1:3194044/1‑100 (MQ=255)
                                                                        aacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGGGACGAACGCCGTGCTCGGGCAAAATTTCACTGGCAATTCCTGGAACAtt                             >  2:3237460/1‑98 (MQ=255)
                                                                           aaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCg                        >  2:1288968/1‑100 (MQ=255)
                                                                            aGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCgg                       >  1:2530931/1‑100 (MQ=255)
                                                                            aGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCgg                       >  2:1305210/1‑100 (MQ=255)
                                                                            aGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCgg                       <  1:1450433/100‑1 (MQ=255)
                                                                            aGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCgg                       >  1:1599030/1‑100 (MQ=255)
                                                                             gggAAGAAGGCGCTGTTGTTGTAGCCAAAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGACAATTCCTGGAACATTATCCggt                      <  2:2001589/100‑2 (MQ=255)
                                                                               gaagaaGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgctg                    <  2:3097160/100‑1 (MQ=255)
                                                                                aagacgGCGCTGTGGCTGTAGCCACAGCCATTGGAGATGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGCCAATTCCTGGAACATTATCCGgctg                    <  2:3102749/99‑1 (MQ=255)
                                                                                    agGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAg               >  1:2225840/1‑100 (MQ=255)
                                                                                    agGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAg               <  2:2468137/100‑1 (MQ=255)
                                                                                    agGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAg               >  1:457123/1‑100 (MQ=255)
                                                                                    aatggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAg               >  1:1387014/5‑100 (MQ=255)
                                                                                     ggCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAgg              >  2:725905/1‑100 (MQ=255)
                                                                                     atggcTGTGGCTACAGCCACAGCCATTGAAGCTTGCGACGGACGCCGTGCTCGAGCAAAATTTCACGGGCAAGGCCTGGAACAGTATCCGGCTGCTCAg               >  2:1957141/4‑99 (MQ=255)
                                                                                     atggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAgg              >  1:2478012/4‑100 (MQ=255)
                                                                                       ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc            >  1:2398646/2‑100 (MQ=255)
                                                                                       ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc            >  1:2831431/2‑100 (MQ=255)
                                                                                       ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATATCACTTGCAATTCCTGGAACATTATCTGGCTGCTCAGGa             >  1:2831730/2‑99 (MQ=255)
                                                                                       ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCCCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc            >  1:1126025/2‑100 (MQ=255)
                                                                                       cgcTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc            >  2:2540753/1‑100 (MQ=255)
                                                                                       cgcTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc            <  2:478590/100‑1 (MQ=255)
                                                                                       cgcTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc            <  2:332246/100‑1 (MQ=255)
                                                                                        gcTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACt           >  1:2364376/1‑100 (MQ=255)
                                                                                        gcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTTCTCAGGACt           >  1:3006072/1‑100 (MQ=255)
                                                                                        gcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACt           >  2:2466309/1‑100 (MQ=255)
                                                                                         cTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTg          >  1:190085/1‑100 (MQ=255)
                                                                                         cTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTg          >  1:2218202/1‑100 (MQ=255)
                                                                                          tgtgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgc         >  1:652768/1‑100 (MQ=255)
                                                                                           gtgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgcg        <  2:2501381/100‑1 (MQ=255)
                                                                                           gtgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgc         <  2:340193/99‑1 (MQ=255)
                                                                                               cTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTta    <  1:1056926/100‑1 (MQ=255)
                                                                                                tGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCCCGAGCAAAATTACACTGGCAAGTCCTGGAACAGTATCCGGCTGCTCAGGACCGCGCGta    >  2:3296228/1‑99 (MQ=255)
                                                                                                 gTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTtata  <  1:2709403/100‑1 (MQ=255)
                                                                                                 gTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTtata  <  1:2958432/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATA  >  NC_000913/1650050‑1650246

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: