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.

AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCG  >  NC_000913/3951460‑3951658
                                                                                                    |                                                                                                    
 aCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgca                                                                                                      <  2:1280690/100‑2 (MQ=255)
 aCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgca                                                                                                      <  2:1696057/100‑2 (MQ=255)
   aaGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgc‑gcg                                                                                                   <  2:979037/100‑1 (MQ=255)
    aGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgcagcg                                                                                                   >  2:1936836/1‑96 (MQ=255)
     gCTGTGTCAGGCACATGTGGCATTACAAAGTGATTTAAAAGCTCTGTCACCAGCAGTACAGCAGCCCTTAAATCAATGACTGGTAGCAACACTGC‑GCGCa                                                                                                 <  1:2662436/100‑1 (MQ=255)
      cTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAg                                                                                                <  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)
             aGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGtg                                                                                         <  1:3078048/100‑1 (MQ=255)
             aGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGt                                                                                          >  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)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGtga                                                                                        >  1:2901914/1‑100 (MQ=255)
                cacaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGtgatg                                                                                      <  2:1408584/100‑1 (MQ=255)
                   aTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAAc                                                                                   >  1:3287568/1‑100 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAAc                                                                                   >  2:1392476/1‑100 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAAc                                                                                   <  2:2468955/100‑1 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACa                                                                                  >  1:1199446/1‑100 (MQ=255)
                     gtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACa                                                                                  >  2:2848709/1‑100 (MQ=255)
                       ggCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAtt                                                                                >  2:2988899/1‑100 (MQ=255)
                        gCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTcc                                                                              >  2:576047/1‑100 (MQ=255)
                        gCATTACAAGGTGATTTAAATGCTCGGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTc                                                                               <  1:3140114/100‑1 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCTGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTcc                                                                              <  2:1155374/100‑1 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCt                                                                             >  2:1455644/1‑100 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTg                                                                            <  2:2875243/100‑1 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTg                                                                            >  2:1486830/1‑100 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTg                                                                            >  2:1867485/1‑100 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGc                                                                          >  1:1033215/1‑100 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCg                                                                         <  1:640791/100‑1 (MQ=255)
                              cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                         <  2:3376819/100‑1 (MQ=255)
                               aaGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGt                                                                        >  1:1160390/1‑100 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                     <  2:2503031/100‑1 (MQ=255)
                                    gATTTAAGTGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGa                                                                   <  1:1264647/100‑1 (MQ=255)
                                    gATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAc                                                                  >  2:605418/1‑100 (MQ=255)
                                     aTTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAc                                                                  <  1:3191129/100‑1 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAc                                                                  >  2:2967364/1‑99 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                                >  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)
                                                  gTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGc                                                    >  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)
                                                   ttACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCt                                                   <  2:689998/100‑1 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                >  2:2164755/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                >  2:3067258/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                >  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)
                                                             tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACgcgc                                         >  1:1023178/1‑100 (MQ=255)
                                                               cagcagCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCg                                       >  1:2718130/1‑100 (MQ=255)
                                                                 gcagcCGTTAAATGAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                      >  2:2351444/1‑100 (MQ=255)
                                                                  cagcCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttg                                    >  1:1260581/1‑100 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                   >  2:2580135/1‑100 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                   >  1:3004369/1‑100 (MQ=255)
                                                                     ccGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                                 <  2:465279/100‑1 (MQ=255)
                                                                      cGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                                 <  1:1136047/100‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                               >  2:2096949/1‑100 (MQ=255)
                                                                        ttAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                               <  2:3436611/100‑1 (MQ=255)
                                                                        ttAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                              <  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)
                                                                          aaaTCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa‑a                           <  1:1558213/100‑1 (MQ=255)
                                                                          aaaTCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa‑a                           <  1:2858627/100‑1 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aaca                        >  1:1783198/1‑100 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aaca                        >  1:3033069/1‑100 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aaca                        >  1:274732/1‑100 (MQ=255)
                                                                               naTGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑aacaa                       >  1:835068/2‑100 (MQ=255)
                                                                                  gACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGt                   >  1:2266561/1‑100 (MQ=255)
                                                                                   aCTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTT‑TTaacaacaacggc                   <  2:2628090/100‑4 (MQ=255)
                                                                                    cTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTc                  >  2:2246934/1‑100 (MQ=255)
                                                                                    cTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTT‑TTaacaacaacggcg                  <  2:1558781/100‑5 (MQ=255)
                                                                                    cTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCg                 <  1:576043/100‑1 (MQ=255)
                                                                                    cTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCg                 >  1:531753/1‑100 (MQ=255)
                                                                                      ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCgg                >  1:2154742/1‑100 (MQ=255)
                                                                                      ggCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACAATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCgg                >  1:704919/1‑100 (MQ=255)
                                                                                      ggCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGa               >  2:88714/1‑100 (MQ=255)
                                                                                      ggCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGa               >  2:2686958/1‑100 (MQ=255)
                                                                                       gcagcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGa               >  2:1434828/1‑100 (MQ=255)
                                                                                        cagcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCgg                >  1:2407355/1‑98 (MQ=255)
                                                                                         agcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATc             <  1:3160516/100‑1 (MQ=255)
                                                                                          gcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCg            >  1:919242/1‑100 (MQ=255)
                                                                                            aaCACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTa          >  1:3052496/1‑100 (MQ=255)
                                                                                                   c‑gcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCg  <  2:3287592/100‑1 (MQ=255)
                                                                                                    |                                                                                                    
AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAA‑AACAACTGTCGGATCGTAAACCTGCG  >  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: ‑