Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,951,559:1 +A 100% 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_0009133,951,5591.A86.7% 332.7 / 46.3 113pseudogene (18/663 nt)ilvGpseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base . (5/10);  new base A (49/49);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGC  >  NC_000913/3951460‑3951657
                                                                                                      |                                                                                                   
aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgc                                                                                                      >  2:2218963/1‑100 (MQ=255)
aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgc                                                                                                      <  1:2096921/100‑1 (MQ=255)
aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgc                                                                                                      >  2:2456137/1‑100 (MQ=255)
 aCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgca                                                                                                     >  2:568491/1‑99 (MQ=255)
    aGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgcagcg                                                                                                  <  2:44251/100‑5 (MQ=255)
    aGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCGGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgcagcg                                                                                                  >  2:131635/1‑96 (MQ=255)
     gCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgcagcgc                                                                                                 <  2:49174/100‑6 (MQ=255)
     gCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGACGATAAATCAATGACTGGCAGCAACACTgcagcgc                                                                                                 >  1:1281028/1‑95 (MQ=255)
           tCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGc                                                                                           >  2:1730986/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGTTGCGtg                                                                                        <  1:451149/100‑1 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                        <  2:1817212/100‑1 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                        >  2:1124995/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                        >  2:2394052/1‑100 (MQ=255)
                cacaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtgat                                                                                      <  2:665813/100‑1 (MQ=255)
                 acttgttGCATTACAATGTGATTTAAATGCTCTGTTACCAGCAGTACAACAGCCTTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtgat                                                                                      <  1:422258/96‑1 (MQ=255)
                  caTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGAAACACTGCAGCGCAGCTGCGtgatga                                                                                    <  1:970198/100‑1 (MQ=255)
                     gtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACa                                                                                 >  1:1712939/1‑100 (MQ=255)
                     gtgGCATTACAAGGTCATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACa                                                                                 <  2:2876367/100‑1 (MQ=255)
                      tgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCATCTGCGTGATGAACAt                                                                                >  2:2870053/1‑100 (MQ=255)
                       ggCATTACAAGGAGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAtt                                                                               >  2:1912962/1‑100 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTcc                                                                             >  1:121036/1‑100 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTcc                                                                             <  1:133918/100‑1 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCt                                                                            >  2:1464792/1‑100 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCCATGACTGGCAGCAACACTGCAGCGCAGCTGCGTAATGAACATTCCTgg                                                                          >  1:1293230/1‑100 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTgg                                                                          <  2:1202652/100‑1 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGc                                                                         >  2:906222/1‑100 (MQ=255)
                              cgaGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                        <  1:1602830/98‑1 (MQ=255)
                              cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                        >  2:498466/1‑100 (MQ=255)
                               aaGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGtt                                                                      >  2:1738439/1‑100 (MQ=255)
                               aaGGTGATTTAAATGCTCTGTTACCAGCACCACAGCAGCCGTTAAATCAATGACTGGCAGCAACACCGCAGCGCAGCTGCGTGATGAACATTCCTGGCGt                                                                       <  1:1756450/100‑1 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTa                                                                     <  1:2864690/100‑1 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTa                                                                     <  1:607713/100‑1 (MQ=255)
                                  gTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                    <  1:1227089/100‑1 (MQ=255)
                                     aTTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAc                                                                 <  1:782892/100‑1 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCTCCGCTGCATGATGAACATTCCTGGCGTTACGAcc                                                                >  1:2753928/1‑100 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAcc                                                                >  1:1503115/1‑100 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAcc                                                                >  1:2577062/1‑100 (MQ=255)
                                        tAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAt                                                              >  1:1728280/1‑100 (MQ=255)
                                        tAAATGCTCTGTTACCAGCATTACAGAAGCCGTTAAATCAATCACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAt                                                              <  2:2775083/100‑1 (MQ=255)
                                           aTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCTGCTGCGTGATGAACATTCCTGGCGTTACGACCATccc                                                           <  1:1556300/100‑1 (MQ=255)
                                             gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                         <  2:2448306/100‑1 (MQ=255)
                                             gCTCTGTTACCACCATTACAGCAGCCGTTCAATCAATCCCTCGCACCACCACTCCAGCCCAGCAGCGTCACGAACCTTCCTGGCGTCACGACCGTCCCg                                                          >  1:66875/1‑99 (MQ=255)
                                              ctctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                        >  2:1232813/1‑100 (MQ=255)
                                                 tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                     >  1:678311/1‑100 (MQ=255)
                                                  gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGCCCATCCCGGTGACGc                                                   >  2:2231818/1‑100 (MQ=255)
                                                  gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGc                                                   <  2:426516/100‑1 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTAt                                                >  2:2203692/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                               <  2:1446745/100‑1 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                               >  2:1099545/1‑100 (MQ=255)
                                                       cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCt                                              >  2:2428478/1‑100 (MQ=255)
                                                        aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCt                                              >  1:1483625/1‑100 (MQ=255)
                                                        aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCt                                              >  2:1449282/1‑100 (MQ=255)
                                                          ccatacacaagcCGTAATATCACTCACTGGCAGCAACACCGCCGCGCATATGCGTGATGAACATTCCACGCGTAACGCCCATCCCGCTGACGCTATc                                               <  1:1308044/88‑1 (MQ=255)
                                                           aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                           <  2:566024/100‑1 (MQ=255)
                                                            ttACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcgc                                          >  1:2542800/1‑100 (MQ=255)
                                                             tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcgcg                                         <  2:2592638/100‑1 (MQ=255)
                                                               cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGcc                                       >  2:77090/1‑100 (MQ=255)
                                                                 gcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                     <  1:1029789/100‑1 (MQ=255)
                                                                   agcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttg                                   <  2:1863621/100‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                  <  1:1349405/100‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                  <  1:1349350/100‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTCCAGCGCAGCTGCGTGATGAACATCCCTGGCGTTACGACCATCCCGGTGCCGCTATCTACGCGCCgttg                                   >  1:2562345/1‑99 (MQ=255)
                                                                     ccGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                 >  2:491419/1‑100 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGGAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                 <  2:1094469/98‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               <  2:815789/100‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               <  2:2003859/100‑1 (MQ=255)
                                                                        ttAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                              <  1:1730978/100‑1 (MQ=255)
                                                                         tAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                             >  2:2544607/1‑100 (MQ=255)
                                                                          aaaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                            <  2:2629428/100‑1 (MQ=255)
                                                                          aaaTCAATGACTAGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                             >  2:841132/1‑99 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                          >  2:1521799/1‑100 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                          <  2:29262/100‑1 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                          <  2:206297/100‑1 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                          >  2:1276516/1‑100 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                             <  1:424260/98‑1 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                       >  1:602327/1‑100 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                        <  1:1739581/99‑1 (MQ=255)
                                                                               aaTGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCGCGGTGACGCTATCTGCGCGCCGTTGTTGTTAAaacaa                      >  1:1963253/1‑100 (MQ=255)
                                                                                aTGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                     <  1:691761/100‑1 (MQ=255)
                                                                                  gACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTTTaacaacaacgg                   <  1:988805/100‑3 (MQ=255)
                                                                                  gACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGt                  <  1:2231804/100‑1 (MQ=255)
                                                                                   aCTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                   <  1:2655793/99‑1 (MQ=255)
                                                                                   aCTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTc                 >  1:2219056/1‑100 (MQ=255)
                                                                                    cTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTTTaacaacaacggcg                 <  2:1404191/100‑5 (MQ=255)
                                                                                    cTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTTTaacaacaacggcg                 <  1:1512116/100‑5 (MQ=255)
                                                                                      ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTTTaacaacaacggc               <  1:896793/100‑7 (MQ=255)
                                                                                      ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg               >  1:2524719/1‑100 (MQ=255)
                                                                                      ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg               >  2:2356298/1‑100 (MQ=255)
                                                                                       gcagcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa              >  1:1601644/1‑100 (MQ=255)
                                                                                        cagcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGAt             <  2:1642286/100‑1 (MQ=255)
                                                                                         agcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc            >  2:173193/1‑100 (MQ=255)
                                                                                          gcagcaCTGCAGCGCAGATACGCGATGCACCCTCCTAAGGTTACGACAATCCCGGTGACGCTATCTACGCGACGTTGTTGTTAAAACAACTGTGGGATCg           <  1:966797/100‑1 (MQ=255)
                                                                                           caACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGt          >  1:2802037/1‑100 (MQ=255)
                                                                                           caACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGt          <  2:1269992/100‑1 (MQ=255)
                                                                                           caACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGt          <  2:891249/100‑1 (MQ=255)
                                                                                            aaCACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg           >  2:2617478/1‑98 (MQ=255)
                                                                                            aaCACTGCAGCGCAGCTGCGTGATCAACATTCCTGGCTTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACACCTGTCGGATCGTa         <  1:1799145/100‑1 (MQ=255)
                                                                                             acacTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaaa       <  2:2540573/100‑1 (MQ=255)
                                                                                               acTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc     <  1:1664745/100‑1 (MQ=255)
                                                                                                   gcgcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGc  >  1:1564373/1‑100 (MQ=255)
                                                                                                      |                                                                                                   
AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGC  >  NC_000913/3951460‑3951657

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