Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,951,598:1 +T 35.8% 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.T35.8% 87.7 / 167.8 95pseudogene (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 . (28/33);  new base T (14/20);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 101 total aligned reads displayed.

TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGAT  >  NC_000913/3951498‑3951696
                                                                                                    |                                                                                                   
tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACATTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                                                      <  2:1286219/100‑1 (MQ=255)
tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                                                      >  1:1949127/1‑100 (MQ=255)
tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                                                      >  1:1258640/1‑100 (MQ=255)
   aaaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAttc                                                                                                  >  2:1771334/1‑98 (MQ=255)
   aaaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAttc                                                                                                  >  2:1578022/1‑98 (MQ=255)
    aaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATccc                                                                                                 >  2:1326595/1‑100 (MQ=255)
      tGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                                                               <  1:1717846/100‑1 (MQ=255)
       gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                                                              <  2:551245/100‑1 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCTTTACGACCATCCCGGTGa                                                                                            >  2:780326/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACg                                                                                          >  1:2024302/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGc                                                                                         >  1:1465942/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGc                                                                                         <  1:576123/100‑1 (MQ=255)
              tACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCTGCTGCGTGAAGAACAGTCCTGGCGGTACGACCATTCCCGGTGACGCt                                                                                        >  2:1586469/1‑100 (MQ=255)
                ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                                                     >  1:170733/1‑100 (MQ=255)
                  aGCATTACATCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCt                                                                                    <  2:972363/100‑1 (MQ=255)
                  aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCt                                                                                    >  2:2129641/1‑100 (MQ=255)
                  aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCt                                                                                    >  2:1955302/1‑100 (MQ=255)
                  aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCt                                                                                    >  2:490338/1‑99 (MQ=255)
                   gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAc                                                                                  >  1:1846127/1‑100 (MQ=255)
                    cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                                                                 <  2:2305735/100‑1 (MQ=255)
                    cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                                                                 >  2:255135/1‑100 (MQ=255)
                     aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTAcg                                                                                 <  1:1251188/100‑1 (MQ=255)
                       tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTAcgcg                                                                               <  1:1478829/100‑1 (MQ=255)
                       tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTAcgcg                                                                               >  1:274507/1‑100 (MQ=255)
                        aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGcc                                                                             <  1:1082932/100‑1 (MQ=255)
                         cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGcc                                                                             <  2:1688872/100‑1 (MQ=255)
                         cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCg                                                                            <  2:1546571/100‑1 (MQ=255)
                            cagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttg                                                                         <  1:822821/100‑1 (MQ=255)
                            cagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttg                                                                         >  2:1270083/1‑100 (MQ=255)
                            cagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttg                                                                         >  2:1270176/1‑100 (MQ=255)
                              gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgt                                                                        <  1:1026540/100‑1 (MQ=255)
                              gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttg                                                                         <  2:1745532/99‑1 (MQ=255)
                               ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCgttgtt                                                                       >  1:980690/1‑100 (MQ=255)
                                  ttAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                                                                   <  1:1515755/100‑1 (MQ=255)
                                    aaaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                                                 <  1:2201128/100‑1 (MQ=255)
                                      aTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                                                               >  1:608234/1‑100 (MQ=255)
                                      aTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                                                               <  1:1992590/100‑1 (MQ=255)
                                       taaaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                                                               <  2:979776/98‑1 (MQ=255)
                                       tCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                                                               <  1:344792/100‑1 (MQ=255)
                                        cAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                                                              <  1:2336291/100‑1 (MQ=255)
                                        cAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                                                             >  2:161541/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                            <  1:1270088/100‑1 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                            <  1:1270181/100‑1 (MQ=255)
                                          aTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                            >  2:2078339/1‑100 (MQ=255)
                                          aTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                            <  1:1955316/100‑1 (MQ=255)
                                          aTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACt                                                           >  1:998982/1‑100 (MQ=255)
                                           tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                                                          >  2:289587/1‑100 (MQ=255)
                                                 gcagcaACACTGCGCGCAGCTGCGTGATGAACATTGCTGTCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                                                       >  2:1965796/1‑97 (MQ=255)
                                                 gcagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgaa                                                     >  2:1429936/1‑98 (MQ=255)
                                                  cagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc                                                   <  2:1987571/100‑1 (MQ=255)
                                                       acacTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaa                                               <  1:1169187/100‑1 (MQ=255)
                                                       acacTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaaa                                              <  2:719230/100‑1 (MQ=255)
                                                         acTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAc                                             <  2:1223632/100‑1 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCt                                           >  2:838440/1‑100 (MQ=255)
                                                            gcgcgcAGCTGCGTGATGAACATTCCTGTCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTg                                          <  2:204766/100‑1 (MQ=255)
                                                            gcgcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGc                                         <  2:1051543/100‑1 (MQ=255)
                                                                gcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGa                                      <  2:2029625/100‑1 (MQ=255)
                                                                gcAGCTGCGTGATGAACATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGAt                                     <  2:1573443/100‑1 (MQ=255)
                                                                  aGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTg                                   >  1:1075613/1‑100 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGc                                  <  1:321289/100‑1 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                                 <  2:1225689/100‑1 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                                 <  2:174744/100‑1 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGCGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                                 <  1:911084/100‑1 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTTTTGTTAAAACAACTGTCGGATCGTAATCCTGCGGATTGcg                                 >  1:34536/1‑100 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                                 >  1:2044396/1‑100 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgt                                <  2:2046412/100‑1 (MQ=255)
                                                                      gCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtc                               <  1:1931074/100‑1 (MQ=255)
                                                                      gCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtc                               >  2:1627178/1‑100 (MQ=255)
                                                                      gCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtc                               >  1:2070031/1‑100 (MQ=255)
                                                                       cGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                              <  2:1075608/100‑1 (MQ=255)
                                                                       cGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                              >  2:1676314/1‑100 (MQ=255)
                                                                       cGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCCGCGGATTGcgt                                <  2:610360/98‑1 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                              >  1:884147/1‑100 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCTGATCGTAAACCTGCGGATTGcgtcg                              <  2:282929/99‑1 (MQ=255)
                                                                         tgatgaACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                            >  2:200675/1‑100 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                            <  2:1222337/100‑1 (MQ=255)
                                                                           atgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAATCAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGa                           <  2:1258635/100‑1 (MQ=255)
                                                                           atgaACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAc                          <  1:2203089/100‑1 (MQ=255)
                                                                            tgaACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                         <  1:780329/100‑1 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                         <  1:2129659/100‑1 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                         >  2:1647515/1‑100 (MQ=255)
                                                                              aaCATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                       <  2:1514559/100‑1 (MQ=255)
                                                                              aaCATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                       <  2:1949113/100‑1 (MQ=255)
                                                                               aCATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                       <  2:1630763/100‑1 (MQ=255)
                                                                                 aTTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGa                    <  1:1454793/100‑1 (MQ=255)
                                                                                  ttCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGa                    >  1:1644059/1‑100 (MQ=255)
                                                                                     cTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGAGACCACAGAtgtg                <  1:1599080/100‑1 (MQ=255)
                                                                                      tGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgg               <  1:1960889/100‑1 (MQ=255)
                                                                                         cGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAtgtg                >  2:1874156/1‑96 (MQ=255)
                                                                                         cGTAACGACCATCCGGGTGACGATCTCTACGCGCCGTTGATGTTAAGACAACTGTAGGATCGTAAACCTGCTGATTGCGTCGTGACCACAGATGTggg              >  2:1046677/1‑98 (MQ=255)
                                                                                           ttACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcag          >  1:86138/1‑100 (MQ=255)
                                                                                             aCGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagca        <  2:2070014/100‑1 (MQ=255)
                                                                                             aCGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagca        >  2:854223/1‑100 (MQ=255)
                                                                                              cGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagc         >  1:1428816/1‑99 (MQ=255)
                                                                                                aCCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCa     <  2:1218480/100‑1 (MQ=255)
                                                                                                aCCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACACATGTGGGGCAGCACCa     >  2:1376394/1‑100 (MQ=255)
                                                                                                 ccATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAg    >  2:1619678/1‑100 (MQ=255)
                                                                                                 ccATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAg    >  1:1205205/1‑100 (MQ=255)
                                                                                                  cATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGa   >  2:1878889/1‑100 (MQ=255)
                                                                                                   aTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGAT  >  2:641770‑M1/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGAT  >  NC_000913/3951498‑3951696

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