Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,951,559	+A	pseudogene (18/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_000913	3,951,559	1	.	A	100.0%	305.1 / NA	85	pseudogene (18/663 nt)	ilvG	pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (46/39);  total (46/39)

AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAA  >  NC_000913/3951460‑3951652                                                                                                     |                                                                                              aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgc                                                                                                <  1:1836179/100‑1 (MQ=255) aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgc                                                                                                <  2:443452/100‑1 (MQ=255)      gCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgcagcgc                                                                                           <  2:148922/100‑6 (MQ=255)       cTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCa                                                                                          <  1:758328/100‑1 (MQ=255)          cGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCt                                                                                       <  2:1492174/100‑1 (MQ=255)             cAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCg                                                                                    <  2:2038116/100‑1 (MQ=255)              aGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGt                                                                                   <  2:156978/100‑1 (MQ=255)               ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                  >  1:628391/1‑100 (MQ=255)               ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                  >  2:1626275/1‑100 (MQ=255)               ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                  >  2:538873/1‑100 (MQ=255)                gCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtga                                                                                 <  2:2427158/100‑1 (MQ=255)                  acaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtgatg                                                                               <  2:487138/100‑1 (MQ=255)                     tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAAc                                                                            >  2:2113527/1‑100 (MQ=255)                       tgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAt                                                                          >  1:1756054/1‑100 (MQ=255)                        ggCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAtt                                                                         >  1:2229569/1‑100 (MQ=255)                          cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTcc                                                                       <  1:2113544/100‑1 (MQ=255)                          cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTcc                                                                       <  1:1280700/100‑1 (MQ=255)                          cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTcc                                                                       <  1:1504558/100‑1 (MQ=255)                           aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCt                                                                      >  1:209529/1‑100 (MQ=255)                            ttACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTg                                                                     <  2:60026/100‑1 (MQ=255)                            ttACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTg                                                                     <  2:718159/100‑1 (MQ=255)                               cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                  <  1:2097528/100‑1 (MQ=255)                               cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                  >  1:191791/1‑100 (MQ=255)                               cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                  >  2:2235205/1‑100 (MQ=255)                               cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                  >  2:2291945/1‑100 (MQ=255)                                aaGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGt                                                                 >  1:763055/1‑100 (MQ=255)                                 aGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGtt                                                                >  1:683753/1‑100 (MQ=255)                                 aGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGtt                                                                <  2:1007906/100‑1 (MQ=255)                                 aGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGtt                                                                >  1:840141/1‑100 (MQ=255)                                   gTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                              <  1:792894/100‑1 (MQ=255)                                       tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAcc                                                          >  1:1951891/1‑100 (MQ=255)                                         tAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAt                                                        >  2:958304/1‑100 (MQ=255)                                         tAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAt                                                        >  1:1545686/1‑100 (MQ=255)                                             tGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCg                                                    >  1:2367742/1‑100 (MQ=255)                                              gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                   <  1:2379345/100‑1 (MQ=255)                                              gCTCTGTTACCAGCATTACAGCAGCCGTTAAATAAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                   <  2:1053356/100‑1 (MQ=255)                                               ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGCACCTTCCTGGCGTTACGACCATCCCGGt                                                  >  2:2419405/1‑100 (MQ=255)                                               ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                  >  2:847818/1‑100 (MQ=255)                                                tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTg                                                 <  2:683748/100‑1 (MQ=255)                                                 ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGa                                                <  1:1170558/100‑1 (MQ=255)                                                  tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                               <  2:369691/100‑1 (MQ=255)                                                  tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                               >  2:967226/1‑100 (MQ=255)                                                   gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACg                                              <  2:238679/100‑1 (MQ=255)                                                   gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACg                                              >  2:2011771/1‑100 (MQ=255)                                                      aCCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTa                                           >  1:2660185/1‑100 (MQ=255)                                                       ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTAt                                          >  2:44428/1‑100 (MQ=255)                                                        cAGCATTACAGCAGCCGTTAAATCAATGACTTGCAGCAACACCGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                         <  1:592987/100‑1 (MQ=255)                                                        cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                         >  2:1782538/1‑100 (MQ=255)                                                           cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAc                                      >  2:2592979/1‑100 (MQ=255)                                                             ttACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcgc                                    >  1:101138/1‑100 (MQ=255)                                                             ttACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcgc                                    >  1:1018360/1‑100 (MQ=255)                                                                 agcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCg                                >  2:1033325/1‑100 (MQ=255)                                                                  gcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                               >  1:2239882/1‑100 (MQ=255)                                                                  gcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                               >  2:2109781/1‑100 (MQ=255)                                                                   cagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgtt                              >  1:449492/1‑100 (MQ=255)                                                                     gcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                            >  2:1748685/1‑100 (MQ=255)                                                                     gcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                            >  2:120601/1‑100 (MQ=255)                                                                      ccGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                            <  2:680116/99‑1 (MQ=255)                                                                       cGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                          <  2:2560717/100‑1 (MQ=255)                                                                          tAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                       >  2:1182281/1‑100 (MQ=255)                                                                            aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                     >  2:2323588/1‑100 (MQ=255)                                                                            aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                     <  2:780990/100‑1 (MQ=255)                                                                            aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                     <  2:415543/100‑1 (MQ=255)                                                                             aTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaaa                    >  1:127029/1‑100 (MQ=255)                                                                              tCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                   <  2:2261180/100‑1 (MQ=255)                                                                              tCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                   >  2:959531/1‑100 (MQ=255)                                                                              tCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                   <  1:253708/100‑1 (MQ=255)                                                                              tCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                   <  1:538877/100‑1 (MQ=255)                                                                               cAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                  <  2:805694/100‑1 (MQ=255)                                                                                aaTGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                 >  2:1514626/1‑100 (MQ=255)                                                                                aaTGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                 >  2:2609242/1‑100 (MQ=255)                                                                                 aTGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                <  1:847824/100‑1 (MQ=255)                                                                                  tGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACt               <  2:2239864/100‑1 (MQ=255)                                                                                   gACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg              <  2:485756/100‑1 (MQ=255)                                                                                    aCTGGCAGCAACACTGCAGCGCAGCTGCGTGATGATCATTCCTGGCGTTACGAACATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACt               <  2:1339633/98‑1 (MQ=255)                                                                                      tGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg           <  2:2362844/100‑1 (MQ=255)                                                                                       ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg          >  1:2304826/1‑100 (MQ=255)                                                                                       ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg          <  2:628386/100‑1 (MQ=255)                                                                                       ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg          >  1:1449907/1‑100 (MQ=255)                                                                                        gcagcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa         >  1:2368737/1‑100 (MQ=255)                                                                                         cagcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGAt        >  2:1762409/1‑100 (MQ=255)                                                                                          agcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc       <  2:315279/100‑1 (MQ=255)                                                                                           gcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg      <  1:1848439/100‑1 (MQ=255)                                                                                           gcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg      <  2:2660165/100‑1 (MQ=255)                                                                                            caACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGt     <  1:2419424/100‑1 (MQ=255)                                                                                            caACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGt     >  1:2472829/1‑100 (MQ=255)                                                                                             aaCACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTa    >  1:245016/1‑100 (MQ=255)                                                                                               cacTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaaa  >  2:284505/1‑100 (MQ=255)                                                                                                     |                                                                                              AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAA  >  NC_000913/3951460‑3951652

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