Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,598 +T 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.T95.3% 442.3 / 12.5 127pseudogene (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 . (4/2);  new base T (60/61);  total (64/63)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAC  >  NC_000913/3951498‑3951691
                                                                                                    |                                                                                               
tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                                                  >  2:2653502/1‑100 (MQ=255)
tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                                                  >  1:1905638/1‑100 (MQ=255)
tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                                                  <  2:677021/100‑1 (MQ=255)
tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                                                  >  1:1180019/1‑100 (MQ=255)
  tAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑tt                                                                                               <  1:831990/100‑2 (MQ=255)
     aTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑ttccc                                                                                            >  2:164339/1‑96 (MQ=255)
       gCTCTGTTACCAGCATTACAGCAGCCGTTAGATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCg‑g                                                                                          <  1:210665/100‑1 (MQ=255)
        ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑Gt                                                                                         >  2:2624423/1‑100 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTg                                                                                        <  2:416320/100‑1 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTg                                                                                        <  1:1360451/100‑1 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTg                                                                                        >  2:1991465/1‑100 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTg                                                                                        <  2:2450387/100‑1 (MQ=255)
          ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGa                                                                                       >  2:2151596/1‑100 (MQ=255)
          ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGa                                                                                       <  2:2334305/100‑1 (MQ=255)
           tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCG‑GTGACg                                                                                     >  1:475840/1‑100 (MQ=255)
           ggTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGAAGCTCCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑ATGa                                                                                       >  2:1239451/2‑99 (MQ=255)
              tACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCt                                                                                   >  1:2162268/1‑100 (MQ=255)
               aCCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTa                                                                                  >  1:1054160/1‑100 (MQ=255)
               aCCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTa                                                                                  >  1:128161/1‑100 (MQ=255)
                ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTAt                                                                                 >  2:1497858/1‑100 (MQ=255)
                 cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATc                                                                                >  2:720386/1‑100 (MQ=255)
                   gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTa                                                                              >  1:299179/1‑100 (MQ=255)
                    cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTAc                                                                             <  1:428802/100‑1 (MQ=255)
                       tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTAcgcg                                                                          <  2:1697958/100‑1 (MQ=255)
                        aCAGCTGCCGTTTAATCAATTATTGTCAGTAACTGTGTTCGCAGCTGCGTTATGAACTTTCCTGGCGTTACGACCATTCCCG‑GTGACTCTATCTAcgctc                                                                         <  2:125390/100‑3 (MQ=255)
                        aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACgcgc                                                                         >  1:1552406/1‑100 (MQ=255)
                          agcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCg                                                                       >  2:1756049/1‑100 (MQ=255)
                          agcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCg                                                                       >  2:2926041/1‑100 (MQ=255)
                           gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCgt                                                                      >  1:175603/1‑100 (MQ=255)
                              gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCgttgt                                                                   >  2:2829668/1‑100 (MQ=255)
                              gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCgttgt                                                                   <  1:1087321/100‑1 (MQ=255)
                               ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCgttgt                                                                   >  1:1674932/1‑99 (MQ=255)
                                    aaaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTaa                                                             >  2:890728/1‑100 (MQ=255)
                                     aaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTaaa                                                            <  1:2720115/100‑1 (MQ=255)
                                     aaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTaaa                                                            <  1:972063/100‑1 (MQ=255)
                                       tCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAaac                                                          >  2:1373967/1‑100 (MQ=255)
                                       tCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                                                         <  2:475837/100‑1 (MQ=255)
                                        cAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                                                         >  1:1758563/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                                                        >  2:1225992/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                       >  2:114015/1‑100 (MQ=255)
                                          aTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                       <  1:2332591/100‑1 (MQ=255)
                                            gACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                                                     <  2:2246923/100‑1 (MQ=255)
                                             aCTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGt                                                    >  1:12459/1‑100 (MQ=255)
                                                 gcagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa                                                >  2:1613709/1‑100 (MQ=255)
                                                 gcagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa                                                >  1:590288/1‑100 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             <  2:1897999/100‑1 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             >  2:846514/1‑100 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             <  1:2653517/100‑1 (MQ=255)
                                                    gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                             <  2:1758552/100‑1 (MQ=255)
                                                     caACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGt                                            >  2:2220607/1‑100 (MQ=255)
                                                     caACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGt                                            >  2:1450935/1‑100 (MQ=255)
                                                      aaCACTGTGCGAAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTa                                           >  2:1565182/1‑100 (MQ=255)
                                                      aaCACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTa                                           <  1:2540914/100‑1 (MQ=255)
                                                       acacTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaa                                          >  2:1537631/1‑100 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc                                       >  2:1155366/1‑100 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc                                       <  1:1313240/100‑1 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc                                       >  1:1701944/1‑100 (MQ=255)
                                                               cgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCgg                                  <  2:590285/100‑1 (MQ=255)
                                                               cgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCgg                                  <  1:637962/100‑1 (MQ=255)
                                                               cgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCgg                                  <  2:1365546/100‑1 (MQ=255)
                                                                gcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTCTTAAAACAACTGTCGGATCGTAAACCTGCGGa                                 <  1:2829684/100‑1 (MQ=255)
                                                                    cTGCATGATGAACTTTCCTGGAGTTACAGCAATTACCG‑GTGACGCTAACTACGCGCCGTTGTTGTTAAAACAACTGTGTGATCGTAAACCTGCGGATTGc                             <  2:547881/100‑1 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                            <  1:2310177/100‑1 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                            >  1:887782/1‑100 (MQ=255)
                                                                      gCGTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgt                           <  2:1919049/100‑1 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                         <  1:846517/100‑1 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                         <  2:721897/100‑1 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                         >  1:2967706/1‑100 (MQ=255)
                                                                         tgatgaACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcgt                        >  2:48918/1‑100 (MQ=255)
                                                                         tgatgaACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcgt                        >  1:1904327/1‑100 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                       >  1:1527365/1‑100 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCA‑TCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcgt                        >  2:342360/1‑98 (MQ=255)
                                                                            tgaACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAc                     <  2:43775/100‑1 (MQ=255)
                                                                            tgaACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAc                     <  1:1327945/100‑1 (MQ=255)
                                                                            tgaACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAc                     <  1:2197919/100‑1 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                    <  1:2220619/100‑1 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                    <  2:408696/100‑1 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATTCACG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                    >  1:1817708/1‑100 (MQ=255)
                                                                              aaCATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACca                   >  1:1789020/1‑100 (MQ=255)
                                                                              aaCATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACca                   >  2:1952943/1‑100 (MQ=255)
                                                                                cATTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcaca                 >  1:524745/1‑100 (MQ=255)
                                                                                 aTTCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAg                <  2:717356/100‑1 (MQ=255)
                                                                                  ttCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGa               <  1:2124888/100‑1 (MQ=255)
                                                                                   tCCTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAt              >  2:1447191/1‑100 (MQ=255)
                                                                                    ccTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAtg             >  1:1549036/1‑100 (MQ=255)
                                                                                    ccTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAtg             <  1:1450943/100‑1 (MQ=255)
                                                                                     cTGGCGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAtgt            >  1:592137/1‑100 (MQ=255)
                                                                                         ngTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgggg        >  1:1452889/2‑100 (MQ=255)
                                                                                         cGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgggg        <  2:43935/100‑1 (MQ=255)
                                                                                         cGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgggg        <  2:43596/100‑1 (MQ=255)
                                                                                         cGTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgggg        <  2:1701933/100‑1 (MQ=255)
                                                                                          gTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  2:2540852/100‑1 (MQ=255)
                                                                                          gTTACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc       <  2:1407733/100‑1 (MQ=255)
                                                                                           ttACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgca      <  1:93519/100‑1 (MQ=255)
                                                                                            tACGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcag     >  2:1612930/1‑100 (MQ=255)
                                                                                             aCGACCA‑TTCCGTATGACGCTATCTACGCGCCGTTGTTGTGAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagc    <  2:2864745/100‑1 (MQ=255)
                                                                                              cGACCATTCCCG‑GTGACGCTATTTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagca   <  2:2017459/100‑1 (MQ=255)
                                                                                              cGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagca   >  1:2580666/1‑100 (MQ=255)
                                                                                              cGACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagca   <  1:1022083/100‑1 (MQ=255)
                                                                                               gACCATTCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAc  >  1:2545580/1‑100 (MQ=255)
                                                                                                    |                                                                                               
TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCG‑GTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAC  >  NC_000913/3951498‑3951691

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: ‑