Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,951,559:1 +A 59.6% 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.A59.6% 389.6 / 277.5 171pseudogene (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 . (41/28);  new base A (50/52);  total (91/80)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 180 total aligned reads displayed.

AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCG  >  NC_000913/3951460‑3951658
                                                                                                    |                                                                                                    
 aCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgca                                                                                                      <  2:1280690/100‑2 (MQ=255)
 aCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgca                                                                                                      <  2:1696057/100‑2 (MQ=255)
   aaGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgcgcg                                                                                                   <  2:979037/100‑1 (MQ=255)
    aGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgcagcg                                                                                                   >  2:1936836/1‑96 (MQ=255)
     gCTGTGTCAGGCACATGTGGCATTACAAAGTGATTTAAAAGCTCTGTCACCAGCAGTACAGCAGCCCTTAAATCAATGACTGGTAGCAACACTGCGCGCa                                                                                                 <  1:2662436/100‑1 (MQ=255)
      cTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAg                                                                                                <  1:2571031/100‑1 (MQ=255)
       tGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAg                                                                                                >  2:651400/1‑100 (MQ=255)
          gTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCGTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGAAGCTg                                                                                             <  1:2240678/100‑1 (MQ=255)
          gTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTg                                                                                             <  2:2212907/100‑1 (MQ=255)
             aGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGt                                                                                          <  1:1595249/100‑1 (MQ=255)
             aGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtg                                                                                         <  1:3078048/100‑1 (MQ=255)
             aGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGt                                                                                          >  1:452959/1‑99 (MQ=255)
              ngCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                         >  1:1433056/2‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                         >  2:1653002/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                         >  1:2391707/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtga                                                                                        >  1:2901914/1‑100 (MQ=255)
                cacaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtgatg                                                                                      <  2:1408584/100‑1 (MQ=255)
                   aTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAAc                                                                                   >  1:3287568/1‑100 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAAc                                                                                   >  2:1392476/1‑100 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAAc                                                                                   <  2:2468955/100‑1 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACa                                                                                  >  1:1199446/1‑100 (MQ=255)
                     gtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACa                                                                                  >  2:2848709/1‑100 (MQ=255)
                       ggCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAtt                                                                                >  2:2988899/1‑100 (MQ=255)
                        gCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTcc                                                                              >  2:576047/1‑100 (MQ=255)
                        gCATTACAAGGTGATTTAAATGCTCGGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTc                                                                               <  1:3140114/100‑1 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCTGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTcc                                                                              <  2:1155374/100‑1 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCt                                                                             >  2:1455644/1‑100 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTg                                                                            <  2:2875243/100‑1 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTg                                                                            >  2:1486830/1‑100 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTg                                                                            >  2:1867485/1‑100 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGc                                                                          >  1:1033215/1‑100 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                         <  1:640791/100‑1 (MQ=255)
                              cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                         <  2:3376819/100‑1 (MQ=255)
                               aaGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGt                                                                        >  1:1160390/1‑100 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                     <  2:2503031/100‑1 (MQ=255)
                                    gATTTAAGTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGa                                                                   <  1:1264647/100‑1 (MQ=255)
                                    gATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAc                                                                  >  2:605418/1‑100 (MQ=255)
                                     aTTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAc                                                                  <  1:3191129/100‑1 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAc                                                                  >  2:2967364/1‑99 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                >  1:3357872/1‑100 (MQ=255)
                                         aaaTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATc                                                              >  2:3186189/1‑100 (MQ=255)
                                             gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                          <  2:290410/100‑1 (MQ=255)
                                             gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                          <  2:955871/100‑1 (MQ=255)
                                                ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGa                                                       >  2:727979/1‑100 (MQ=255)
                                                ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGa                                                       <  2:1060429/100‑1 (MQ=255)
                                                 tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                      <  2:2189279/100‑1 (MQ=255)
                                                 tGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                      >  2:1844938/1‑100 (MQ=255)
                                                 tGTTACCAGAAGTAACCCAGACGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAAGATTCCTGGCGTTACGACCATCCCGGTGAc                                                      <  1:2190209/100‑1 (MQ=255)
                                                  gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACg                                                     >  1:818750/1‑100 (MQ=255)
                                                  gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACg                                                     >  1:160122/1‑100 (MQ=255)
                                                  gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGc                                                    >  1:600381/1‑100 (MQ=255)
                                                   ttACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAATACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGc                                                    <  1:2218385/100‑1 (MQ=255)
                                                   ttACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGc                                                    <  1:498831/100‑1 (MQ=255)
                                                   ttACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGc                                                    <  2:3200064/100‑1 (MQ=255)
                                                   ttACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCt                                                   <  2:689998/100‑1 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                >  2:2164755/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                >  2:3067258/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                >  2:1101266/1‑100 (MQ=255)
                                                       cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                <  1:1903228/100‑1 (MQ=255)
                                                        aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCt                                               >  1:2777684/1‑100 (MQ=255)
                                                           aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                            <  2:402385/100‑1 (MQ=255)
                                                           aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                            <  1:3431922/100‑1 (MQ=255)
                                                           aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                            >  2:1043858/1‑100 (MQ=255)
                                                           aTTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                            >  2:1043749/1‑100 (MQ=255)
                                                             tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACgcgc                                         >  1:1023178/1‑100 (MQ=255)
                                                               cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCg                                       >  1:2718130/1‑100 (MQ=255)
                                                                 gcagcCGTTAAATGAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                      >  2:2351444/1‑100 (MQ=255)
                                                                  cagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttg                                    >  1:1260581/1‑100 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                   >  2:2580135/1‑100 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                   >  1:3004369/1‑100 (MQ=255)
                                                                     ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                                 <  2:465279/100‑1 (MQ=255)
                                                                      cGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                                 <  1:1136047/100‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                               >  2:2096949/1‑100 (MQ=255)
                                                                        ttAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                               <  2:3436611/100‑1 (MQ=255)
                                                                        ttAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                              <  1:1725353/100‑1 (MQ=255)
                                                                        ntAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                               >  1:738655/2‑100 (MQ=255)
                                                                          aaaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                             <  2:2557349/100‑1 (MQ=255)
                                                                          aaaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                             <  1:1348867/100‑1 (MQ=255)
                                                                          aaaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                           <  1:1558213/100‑1 (MQ=255)
                                                                          aaaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                           <  1:2858627/100‑1 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                        >  1:1783198/1‑100 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                        >  1:3033069/1‑100 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                        >  1:274732/1‑100 (MQ=255)
                                                                               naTGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                       >  1:835068/2‑100 (MQ=255)
                                                                                  gACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGt                   >  1:2266561/1‑100 (MQ=255)
                                                                                   aCTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTTTaacaacaacggc                   <  2:2628090/100‑4 (MQ=255)
                                                                                    cTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTc                  >  2:2246934/1‑100 (MQ=255)
                                                                                    cTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTTTaacaacaacggcg                  <  2:1558781/100‑5 (MQ=255)
                                                                                    cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 <  1:576043/100‑1 (MQ=255)
                                                                                    cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 >  1:531753/1‑100 (MQ=255)
                                                                                      ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg                >  1:2154742/1‑100 (MQ=255)
                                                                                      ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACAATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg                >  1:704919/1‑100 (MQ=255)
                                                                                      ggCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa               >  2:88714/1‑100 (MQ=255)
                                                                                      ggCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa               >  2:2686958/1‑100 (MQ=255)
                                                                                       gcagcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa               >  2:1434828/1‑100 (MQ=255)
                                                                                        cagcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCgg                >  1:2407355/1‑98 (MQ=255)
                                                                                         agcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc             <  1:3160516/100‑1 (MQ=255)
                                                                                          gcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg            >  1:919242/1‑100 (MQ=255)
                                                                                            aaCACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTa          >  1:3052496/1‑100 (MQ=255)
                                                                                                   cgcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg  <  2:3287592/100‑1 (MQ=255)
                                                                                                    |                                                                                                    
AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCG  >  NC_000913/3951460‑3951658

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: