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.

TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑T‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAC  >  NC_000913/3951498‑3951691
                                                                                                    |                                                                                                
 ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑t                                                                                                 >  1:309693/1‑100 (MQ=255)
 ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGACGTTACGACCA‑t                                                                                                 >  1:1903803/1‑100 (MQ=255)
   aaaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑ttc                                                                                               >  2:2004246/1‑98 (MQ=255)
        ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGt                                                                                          >  2:2941367/1‑100 (MQ=255)
        ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGt                                                                                          <  1:1964659/100‑1 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTg                                                                                         <  1:3092206/100‑1 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTg                                                                                         <  1:2318448/100‑1 (MQ=255)
          ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGa                                                                                        >  1:3073663/1‑100 (MQ=255)
           tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGAc                                                                                       <  1:2561457/100‑1 (MQ=255)
           tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGAc                                                                                       >  1:2340596/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACg                                                                                      <  2:443528/100‑1 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACg                                                                                      >  2:1752582/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCACCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACg                                                                                      >  2:909260/1‑100 (MQ=255)
              tACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATCAACATTCATGGCGTTACGACCATT‑CCCGGTGACGCt                                                                                    >  2:1748649/1‑100 (MQ=255)
               aCCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTa                                                                                   >  1:3248898/1‑100 (MQ=255)
                ccAGCATTACAGCAGCCGTTAAATCAATGCCTGGCAGCAACACTGCGCGCAGCTGCGTGATAAACATTCCTAGCGTTACGACCATT‑CCCGGTGACGCTAt                                                                                  >  1:3530235/1‑100 (MQ=255)
                ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCCGCAACACTGCGCGCAGCTGCGTGATTAACATTCCTGGCGATAAGACCATT‑CCCGGTGACGCTAt                                                                                  >  1:1815421/1‑100 (MQ=255)
                  aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCt                                                                                >  1:2603609/1‑100 (MQ=255)
                  aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCt                                                                                >  1:338443/1‑100 (MQ=255)
                   gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTa                                                                               <  1:293721/100‑1 (MQ=255)
                   gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTa                                                                               >  1:1044697/1‑100 (MQ=255)
                     aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTAcg                                                                             >  2:3302710/1‑100 (MQ=255)
                     aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTAcg                                                                             >  2:3921472/1‑100 (MQ=255)
                       tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCTGGTGACGCTATCTAcgcg                                                                           <  2:1159779/100‑1 (MQ=255)
                        aCAGTAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACgcgc                                                                          <  2:2478255/100‑1 (MQ=255)
                        aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGCTGAAAATTCATGGCGTTACGACCATT‑ACCGGTGACGATATCTACgcgc                                                                          >  1:2275820/1‑100 (MQ=255)
                        aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACgcgc                                                                          >  1:3186861/1‑100 (MQ=255)
                         cagcagCCGGTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGcc                                                                         >  2:3495512/1‑100 (MQ=255)
                          agcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑TCCGGTGACGCTATCTACGCGCCg                                                                        >  1:531706/1‑100 (MQ=255)
                           gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgt                                                                       <  2:253826/100‑1 (MQ=255)
                           gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgt                                                                       >  1:1206473/1‑100 (MQ=255)
                           gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgt                                                                       >  2:347380/1‑100 (MQ=255)
                               ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgtt                                                                   >  2:2707392/1‑100 (MQ=255)
                               ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgtt                                                                   >  1:1367975/1‑100 (MQ=255)
                                cGTTATATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgt                                                                    <  2:2257930/98‑1 (MQ=255)
                                cGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgttt                                                                  >  1:366170/1‑99 (MQ=255)
                                cGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgttg                                                                  >  1:2013412/1‑100 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  1:2777062/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  2:1068749/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  1:3286532/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  2:2197246/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 <  1:1586643/100‑1 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCgttgttgt                                                                 >  1:3775024/1‑100 (MQ=255)
                                 gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCGACGCGCCgttgttgt                                                                 >  1:3866943/1‑100 (MQ=255)
                                    aaaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                                                              <  2:1960896/100‑1 (MQ=255)
                                     aaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa‑a                                                            >  2:1717821/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aacaa                                                        >  2:2333649/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aacaa                                                        <  1:2863253/100‑1 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aacaa                                                        <  1:2060962/100‑1 (MQ=255)
                                          aTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aacaac                                                       <  2:2084550/100‑1 (MQ=255)
                                           tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aacaacn                                                      <  2:3201335/100‑2 (MQ=255)
                                           tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aacaac                                                       >  2:86947/1‑99 (MQ=255)
                                            gACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTg                                                     <  1:3317136/100‑1 (MQ=255)
                                            gACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTAACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTg                                                     >  1:1499285/1‑100 (MQ=255)
                                             aCTGGCAGCAACACTGCGCGCAGCTGCGTGAGGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTg                                                     <  2:1598414/99‑1 (MQ=255)
                                              cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTc                                                   >  1:943977/1‑100 (MQ=255)
                                               tGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCg                                                  <  1:3262075/100‑1 (MQ=255)
                                                ggCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCACCATTGTTGTTAA‑AACAACTGTCgg                                                 >  1:3024132/1‑100 (MQ=255)
                                                 gcagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGa                                                <  2:2369013/100‑1 (MQ=255)
                                                  cagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTT‑TTaacaacaacggcgcgta                                               <  2:894623/100‑9 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATTAACATTCCTGGCGATACGACCATT‑CCCGCTGACGCCATCTACGCGCCTTTGTTGTTAA‑AACAACtctc                                                   >  2:1044161/1‑94 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCg                                             >  2:238942/1‑100 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCg                                             >  1:2626806/1‑100 (MQ=255)
                                                     caACACTTCGCGCAGCTGCGTGATGAACACTCCTTGCGTTACAACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCg                                             <  1:634598/99‑1 (MQ=255)
                                                     caACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGt                                            <  2:749548/100‑1 (MQ=255)
                                                       acacTGCTCGAAGCTGCGTGATGAAAATTCCTGGCAGTGCGAAAATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTaa                                          <  2:362147/100‑1 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAAcc                                       >  1:3587886/1‑100 (MQ=255)
                                                           tGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCt                                      <  1:2941393/100‑1 (MQ=255)
                                                               cgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCgg                                  <  1:3188787/100‑1 (MQ=255)
                                                                gtaGCTGCGTGATGAACATTCTTGGCGTTATGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCgg                                  <  1:2589243/97‑1 (MQ=255)
                                                                gcAGCTGCGTGATGAACAGTCCTGGGGTTAGGACCA‑TGCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGa                                 <  2:3744359/100‑1 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGc                             >  2:3125838/1‑100 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGc                             <  2:3786909/100‑1 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                         <  1:1794041/100‑1 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                       <  1:488553/100‑1 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCCTAAACCTGCGGATTGCGTCGTg                       <  2:2268149/100‑1 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCAACGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGcgt                           >  2:3439114/1‑96 (MQ=255)
                                                                            tgaACATTCCTGGCGTTACGACCATT‑CCCGGTGACACTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAc                     >  2:3037781/1‑100 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                    <  1:1752597/100‑1 (MQ=255)
                                                                               aCATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGCTCGTAAACCTGCGGATTGCGTCGTGACcac                  >  2:1349333/1‑100 (MQ=255)
                                                                               aCATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACAGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                  >  1:887676/1‑100 (MQ=255)
                                                                                cATTCCTTGCGTTACGACCATT‑CCCGGTGACGCGATCTACGCGCCGTTGTTGGTAA‑AACCACTGTCGGATCGTAAACCTGCCGATTGCGTCGTGACcac                  <  1:1346715/99‑1 (MQ=255)
                                                                                cATTCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcaca                 >  2:3669614/1‑100 (MQ=255)
                                                                                 aTTCCTGGCGTTACGACCA‑T‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGa               <  1:1032828/100‑1 (MQ=255)
                                                                                  ttCCTGGCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                  >  2:1543188/1‑97 (MQ=255)
                                                                                       ggCGTTACGACCATT‑CCCGGTGACGCTATCTACTCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgg          <  1:3614904/100‑1 (MQ=255)
                                                                                       ggCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCCTGACCACAGATGTgg          >  2:3851643/1‑100 (MQ=255)
                                                                                        gCGTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTggg         <  1:2815826/100‑1 (MQ=255)
                                                                                          gTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  2:2775728/100‑1 (MQ=255)
                                                                                          gTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  2:764720/100‑1 (MQ=255)
                                                                                          gTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       >  1:377932/1‑100 (MQ=255)
                                                                                          gTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       >  2:605301/1‑100 (MQ=255)
                                                                                          gTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  2:3626539/100‑1 (MQ=255)
                                                                                          gTTACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  1:1909616/100‑1 (MQ=255)
                                                                                           ttACGACTATT‑CCCGGTGACTCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGGCCACAGATGTGGGgca      <  1:3672065/100‑1 (MQ=255)
                                                                                           ttACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgca      >  1:403111/1‑100 (MQ=255)
                                                                                           ttACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgca      <  1:412204/100‑1 (MQ=255)
                                                                                           ttACGACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgca      <  2:2929667/100‑1 (MQ=255)
                                                                                            tACGACCATT‑CCCAGCGACGGTATCTACGGCCCGTTGTTGTTAA‑AACAACTATCTGATCGTAAACCTGCGGATTCCGTCGTGTCCACAGATGTggggt       <  2:962508/98‑2 (MQ=255)
                                                                                               gACCATT‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAc  >  1:1429466/1‑100 (MQ=255)
                                                                                                    |                                                                                                
TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑T‑CCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAC  >  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: ‑