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.

TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGAT  >  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)
   aaaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑ttc                                                                                                  >  2:1771334/1‑98 (MQ=255)
   aaaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑ttc                                                                                                  >  2:1578022/1‑98 (MQ=255)
    aaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑Tccc                                                                                                 >  2:1326595/1‑100 (MQ=255)
      tGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCgg                                                                                               <  1:1717846/100‑1 (MQ=255)
       gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGt                                                                                              <  2:551245/100‑1 (MQ=255)
         tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCTTTACGACCA‑TCCCGGTGa                                                                                            >  2:780326/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACg                                                                                          >  1:2024302/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGc                                                                                         >  1:1465942/1‑100 (MQ=255)
            gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGc                                                                                         <  1:576123/100‑1 (MQ=255)
              tACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCTGCTGCGTGAAGAACAGTCCTGGCGGTACGACCATTCCCGGTGACGCt                                                                                        >  2:1586469/1‑100 (MQ=255)
                ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATc                                                                                     >  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)
                  aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCt                                                                                    >  2:490338/1‑99 (MQ=255)
                   gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTAc                                                                                  >  1:1846127/1‑100 (MQ=255)
                    cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTAcg                                                                                 <  2:2305735/100‑1 (MQ=255)
                    cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTAcg                                                                                 >  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)
                        aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGcc                                                                             <  1:1082932/100‑1 (MQ=255)
                         cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGcc                                                                             <  2:1688872/100‑1 (MQ=255)
                         cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCg                                                                            <  2:1546571/100‑1 (MQ=255)
                            cagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCgttg                                                                         <  1:822821/100‑1 (MQ=255)
                            cagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCgttg                                                                         >  2:1270083/1‑100 (MQ=255)
                            cagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCgttg                                                                         >  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)
                                  ttAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                                                                   <  1:1515755/100‑1 (MQ=255)
                                    aaaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                                                 <  1:2201128/100‑1 (MQ=255)
                                      aTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                                                               >  1:608234/1‑100 (MQ=255)
                                      aTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                                                               <  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)
                                        cAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                                                             >  2:161541/1‑100 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                            <  1:1270088/100‑1 (MQ=255)
                                         aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                            <  1:1270181/100‑1 (MQ=255)
                                          aTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                            >  2:2078339/1‑100 (MQ=255)
                                          aTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                                                            <  1:1955316/100‑1 (MQ=255)
                                          aTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACt                                                           >  1:998982/1‑100 (MQ=255)
                                           tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                                                          >  2:289587/1‑100 (MQ=255)
                                                 gcagcaACACTGCGCGCAGCTGCGTGATGAACATTGCTGTCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                                                       >  2:1965796/1‑97 (MQ=255)
                                                 gcagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgaa                                                     >  2:1429936/1‑98 (MQ=255)
                                                  cagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc                                                   <  2:1987571/100‑1 (MQ=255)
                                                       acacTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaa                                               <  1:1169187/100‑1 (MQ=255)
                                                       acacTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaaa                                              <  2:719230/100‑1 (MQ=255)
                                                         acTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAc                                             <  2:1223632/100‑1 (MQ=255)
                                                          cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCt                                           >  2:838440/1‑100 (MQ=255)
                                                            gcgcgcAGCTGCGTGATGAACATTCCTGTCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTg                                          <  2:204766/100‑1 (MQ=255)
                                                            gcgcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGc                                         <  2:1051543/100‑1 (MQ=255)
                                                                gcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGa                                      <  2:2029625/100‑1 (MQ=255)
                                                                gcAGCTGCGTGATGAACATGCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGAt                                     <  2:1573443/100‑1 (MQ=255)
                                                                  aGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTg                                   >  1:1075613/1‑100 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGc                                  <  1:321289/100‑1 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                                 <  2:1225689/100‑1 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                                 <  2:174744/100‑1 (MQ=255)
                                                                    cTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGCGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                                 <  1:911084/100‑1 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTTTTGTTAAAACAACTGTCGGATCGTAATCCTGCGGATTGcg                                 >  1:34536/1‑100 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcg                                 >  1:2044396/1‑100 (MQ=255)
                                                                     tGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgt                                <  2:2046412/100‑1 (MQ=255)
                                                                      gCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtc                               <  1:1931074/100‑1 (MQ=255)
                                                                      gCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtc                               >  2:1627178/1‑100 (MQ=255)
                                                                      gCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtc                               >  1:2070031/1‑100 (MQ=255)
                                                                       cGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                              <  2:1075608/100‑1 (MQ=255)
                                                                       cGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                              >  2:1676314/1‑100 (MQ=255)
                                                                       cGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCCGCGGATTGcgt                                <  2:610360/98‑1 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgtcg                              >  1:884147/1‑100 (MQ=255)
                                                                        gTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCTGATCGTAAACCTGCGGATTGcgtcg                              <  2:282929/99‑1 (MQ=255)
                                                                         tgatgaACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                            >  2:200675/1‑100 (MQ=255)
                                                                          gatgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                            <  2:1222337/100‑1 (MQ=255)
                                                                           atgaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAATCAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGa                           <  2:1258635/100‑1 (MQ=255)
                                                                           atgaACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAc                          <  1:2203089/100‑1 (MQ=255)
                                                                            tgaACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                         <  1:780329/100‑1 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                         <  1:2129659/100‑1 (MQ=255)
                                                                             gaACATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                         >  2:1647515/1‑100 (MQ=255)
                                                                              aaCATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                       <  2:1514559/100‑1 (MQ=255)
                                                                              aaCATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                       <  2:1949113/100‑1 (MQ=255)
                                                                               aCATTCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                       <  2:1630763/100‑1 (MQ=255)
                                                                                 aTTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGa                    <  1:1454793/100‑1 (MQ=255)
                                                                                  ttCCTGGCGTTACGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGa                    >  1:1644059/1‑100 (MQ=255)
                                                                                     cTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGAGACCACAGAtgtg                <  1:1599080/100‑1 (MQ=255)
                                                                                      tGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgg               <  1:1960889/100‑1 (MQ=255)
                                                                                         cGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAtgtg                >  2:1874156/1‑96 (MQ=255)
                                                                                         cGTAACGACCA‑TCCGGGTGACGATCTCTACGCGCCGTTGATGTTAAGACAACTGTAGGATCGTAAACCTGCTGATTGCGTCGTGACCACAGATGTggg              >  2:1046677/1‑98 (MQ=255)
                                                                                           ttACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcag          >  1:86138/1‑100 (MQ=255)
                                                                                             aCGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagca        <  2:2070014/100‑1 (MQ=255)
                                                                                             aCGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagca        >  2:854223/1‑100 (MQ=255)
                                                                                              cGACCATTCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagc         >  1:1428816/1‑99 (MQ=255)
                                                                                                aCCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCa     <  2:1218480/100‑1 (MQ=255)
                                                                                                aCCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACACATGTGGGGCAGCACCa     >  2:1376394/1‑100 (MQ=255)
                                                                                                 ccA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAg    >  2:1619678/1‑100 (MQ=255)
                                                                                                 ccA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAg    >  1:1205205/1‑100 (MQ=255)
                                                                                                  cA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGa   >  2:1878889/1‑100 (MQ=255)
                                                                                                   a‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGAT  >  2:641770‑M1/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
TTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCA‑TCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGAT  >  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: ‑