Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,951,598:1 +T 100% pseudogene (57/663 nt) ilvG → pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,5971.T97.5% 501.2 / 4.6 158pseudogene (56/663 nt)ilvGpseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base . (0/4);  new base T (76/78);  total (76/82)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAC  >  NC_000913/3951498‑3951691
                                                                                                    |                                                                                                
 ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAt                                                                                                 >  1:309693/1‑100 (MQ=255)
 ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGACGTTACGACCAt                                                                                                 >  1:1903803/1‑100 (MQ=255)
   aaaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAttc                                                                                               >  2:2004246/1‑98 (MQ=255)
        ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGt                                                                                          >  2:2941367/1‑100 (MQ=255)
        ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGt                                                                                          <  1:1964659/100‑1 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTg                                                                                         <  1:3092206/100‑1 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTg                                                                                         <  1:2318448/100‑1 (MQ=255)
          ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGa                                                                                        >  1:3073663/1‑100 (MQ=255)
           tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGAc                                                                                       <  1:2561457/100‑1 (MQ=255)
           tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGAc                                                                                       >  1:2340596/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACg                                                                                      <  2:443528/100‑1 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACg                                                                                      >  2:1752582/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCACCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACg                                                                                      >  2:909260/1‑100 (MQ=255)
              tACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATCAACATTCATGGCGTTACGACCATTCCCGGTGACGCt                                                                                    >  2:1748649/1‑100 (MQ=255)
               aCCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTa                                                                                   >  1:3248898/1‑100 (MQ=255)
                ccAGCATTACAGCAGCCGTTAAATCAATGCCTGGCAGCAACACTGCGCGCAGCTGCGTGATAAACATTCCTAGCGTTACGACCATTCCCGGTGACGCTAt                                                                                  >  1:3530235/1‑100 (MQ=255)
                ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCCGCAACACTGCGCGCAGCTGCGTGATTAACATTCCTGGCGATAAGACCATTCCCGGTGACGCTAt                                                                                  >  1:1815421/1‑100 (MQ=255)
                  aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCt                                                                                >  1:2603609/1‑100 (MQ=255)
                  aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCt                                                                                >  1:338443/1‑100 (MQ=255)
                   gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTa                                                                               <  1:293721/100‑1 (MQ=255)
                   gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTa                                                                               >  1:1044697/1‑100 (MQ=255)
                     aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTAcg                                                                             >  2:3302710/1‑100 (MQ=255)
                     aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTAcg                                                                             >  2:3921472/1‑100 (MQ=255)
                       tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCTGGTGACGCTATCTAcgcg                                                                           <  2:1159779/100‑1 (MQ=255)
                        aCAGTAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACgcgc                                                                          <  2:2478255/100‑1 (MQ=255)
                        aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGCTGAAAATTCATGGCGTTACGACCATTACCGGTGACGATATCTACgcgc                                                                          >  1:2275820/1‑100 (MQ=255)
                        aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACgcgc                                                                          >  1:3186861/1‑100 (MQ=255)
                         cagcagCCGGTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGcc                                                                         >  2:3495512/1‑100 (MQ=255)
                          agcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTTCCGGTGACGCTATCTACGCGCCg                                                                        >  1:531706/1‑100 (MQ=255)
                           gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgt                                                                       <  2:253826/100‑1 (MQ=255)
                           gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgt                                                                       >  1:1206473/1‑100 (MQ=255)
                           gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgt                                                                       >  2:347380/1‑100 (MQ=255)
                               ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgtt                                                                   >  2:2707392/1‑100 (MQ=255)
                               ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgtt                                                                   >  1:1367975/1‑100 (MQ=255)
                                cGTTATATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgt                                                                    <  2:2257930/98‑1 (MQ=255)
                                cGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttt                                                                  >  1:366170/1‑99 (MQ=255)
                                cGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttg                                                                  >  1:2013412/1‑100 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  1:2777062/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  2:1068749/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  1:3286532/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  2:2197246/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  1:1586643/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 >  1:3775024/1‑100 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCGACGCGCCgttgttgt                                                                 >  1:3866943/1‑100 (MQ=255)
                                    aaaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                                                              <  2:1960896/100‑1 (MQ=255)
                                     aaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                                            >  2:1717821/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                                                        >  2:2333649/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                                                        <  1:2863253/100‑1 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                                                        <  1:2060962/100‑1 (MQ=255)
                                          aTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                       <  2:2084550/100‑1 (MQ=255)
                                           tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaacn                                                      <  2:3201335/100‑2 (MQ=255)
                                           tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                       >  2:86947/1‑99 (MQ=255)
                                            gACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                                                     <  1:3317136/100‑1 (MQ=255)
                                            gACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTAACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                                                     >  1:1499285/1‑100 (MQ=255)
                                             aCTGGCAGCAACACTGCGCGCAGCTGCGTGAGGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                                                     <  2:1598414/99‑1 (MQ=255)
                                              cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTc                                                   >  1:943977/1‑100 (MQ=255)
                                               tGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                                                  <  1:3262075/100‑1 (MQ=255)
                                                ggCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCACCATTGTTGTTAAAACAACTGTCgg                                                 >  1:3024132/1‑100 (MQ=255)
                                                 gcagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa                                                <  2:2369013/100‑1 (MQ=255)
                                                  cagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTTTaacaacaacggcgcgta                                               <  2:894623/100‑9 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATTAACATTCCTGGCGATACGACCATTCCCGCTGACGCCATCTACGCGCCTTTGTTGTTAAAACAACtctc                                                   >  2:1044161/1‑94 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             >  2:238942/1‑100 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             >  1:2626806/1‑100 (MQ=255)
                                                     caACACTTCGCGCAGCTGCGTGATGAACACTCCTTGCGTTACAACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             <  1:634598/99‑1 (MQ=255)
                                                     caACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGt                                            <  2:749548/100‑1 (MQ=255)
                                                       acacTGCTCGAAGCTGCGTGATGAAAATTCCTGGCAGTGCGAAAATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaa                                          <  2:362147/100‑1 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc                                       >  1:3587886/1‑100 (MQ=255)
                                                           tGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCt                                      <  1:2941393/100‑1 (MQ=255)
                                                               cgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCgg                                  <  1:3188787/100‑1 (MQ=255)
                                                                gtaGCTGCGTGATGAACATTCTTGGCGTTATGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCgg                                  <  1:2589243/97‑1 (MQ=255)
                                                                gcAGCTGCGTGATGAACAGTCCTGGGGTTAGGACCATGCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGa                                 <  2:3744359/100‑1 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGc                             >  2:3125838/1‑100 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGc                             <  2:3786909/100‑1 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                         <  1:1794041/100‑1 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                       <  1:488553/100‑1 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCCTAAACCTGCGGATTGCGTCGTg                       <  2:2268149/100‑1 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCAACGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgt                           >  2:3439114/1‑96 (MQ=255)
                                                                            tgaACATTCCTGGCGTTACGACCATTCCCGGTGACACTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAc                     >  2:3037781/1‑100 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                    <  1:1752597/100‑1 (MQ=255)
                                                                               aCATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGCTCGTAAACCTGCGGATTGCGTCGTGACcac                  >  2:1349333/1‑100 (MQ=255)
                                                                               aCATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACAGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                  >  1:887676/1‑100 (MQ=255)
                                                                                cATTCCTTGCGTTACGACCATTCCCGGTGACGCGATCTACGCGCCGTTGTTGGTAAAACCACTGTCGGATCGTAAACCTGCCGATTGCGTCGTGACcac                  <  1:1346715/99‑1 (MQ=255)
                                                                                cATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcaca                 >  2:3669614/1‑100 (MQ=255)
                                                                                 aTTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGa               <  1:1032828/100‑1 (MQ=255)
                                                                                  ttCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                  >  2:1543188/1‑97 (MQ=255)
                                                                                       ggCGTTACGACCATTCCCGGTGACGCTATCTACTCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgg          <  1:3614904/100‑1 (MQ=255)
                                                                                       ggCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCCTGACCACAGATGTgg          >  2:3851643/1‑100 (MQ=255)
                                                                                        gCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTggg         <  1:2815826/100‑1 (MQ=255)
                                                                                          gTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  2:2775728/100‑1 (MQ=255)
                                                                                          gTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  2:764720/100‑1 (MQ=255)
                                                                                          gTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       >  1:377932/1‑100 (MQ=255)
                                                                                          gTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       >  2:605301/1‑100 (MQ=255)
                                                                                          gTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  2:3626539/100‑1 (MQ=255)
                                                                                          gTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  1:1909616/100‑1 (MQ=255)
                                                                                           ttACGACTATTCCCGGTGACTCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGGCCACAGATGTGGGgca      <  1:3672065/100‑1 (MQ=255)
                                                                                           ttACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgca      >  1:403111/1‑100 (MQ=255)
                                                                                           ttACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgca      <  1:412204/100‑1 (MQ=255)
                                                                                           ttACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgca      <  2:2929667/100‑1 (MQ=255)
                                                                                            tACGACCATTCCCAGCGACGGTATCTACGGCCCGTTGTTGTTAAAACAACTATCTGATCGTAAACCTGCGGATTCCGTCGTGTCCACAGATGTggggt       <  2:962508/98‑2 (MQ=255)
                                                                                               gACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAc  >  1:1429466/1‑100 (MQ=255)
                                                                                                    |                                                                                                
TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAC  >  NC_000913/3951498‑3951691

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