Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,850,750 A→G intergenic (‑135/‑793) ivbL ← / → emrD ilvB operon leader peptide/2‑module integral membrane pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,850,7500AG99.6% 693.8 / NA 224intergenic (‑135/‑793)ivbL/emrDilvB operon leader peptide/2‑module integral membrane pump
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (112/110);  total (112/112)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 239 total aligned reads displayed.

CAGCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTCC  >  NC_000913/3850660‑3850839
                                                                                          |                                                                                         
caGCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGcc                                                                                         <  1:6654095/93‑1 (MQ=255)
 aGCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGccc                                                                                        >  1:531906/1‑93 (MQ=255)
  gCTTTACAGGGGAGACAATGGAAAATTTTTGAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGcccc                                                                                       <  1:4621638/93‑1 (MQ=255)
  gCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGcccc                                                                                       >  1:3264981/1‑93 (MQ=255)
  gCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGcccc                                                                                       >  2:1799700/1‑93 (MQ=255)
   cTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccg                                                                                      >  1:105867/1‑93 (MQ=255)
    tttACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccgc                                                                                     >  2:6623930/1‑93 (MQ=255)
    tttACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccgc                                                                                     >  1:4838125/1‑93 (MQ=255)
    tttACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccgc                                                                                     <  2:211279/93‑1 (MQ=255)
    tttACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccgc                                                                                     <  2:1832010/93‑1 (MQ=255)
     ttACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccgcc                                                                                    <  2:3597895/93‑1 (MQ=255)
      tACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccgccg                                                                                   <  2:6535468/93‑1 (MQ=255)
       aCAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCgg                                                                                  <  2:1672992/93‑1 (MQ=255)
        cAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGa                                                                                 >  2:3370449/1‑93 (MQ=255)
        cAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGa                                                                                 >  2:4032233/1‑93 (MQ=255)
          ggggAGACAATGGAAAATTTTTCGGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACCAAAAGCCCCGCCGGAGc                                                                               >  1:1456120/1‑93 (MQ=255)
          ggggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAgg                                                                               >  1:2321748/1‑92 (MQ=255)
          ggggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGc                                                                               >  1:5400928/1‑93 (MQ=255)
          ggggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGc                                                                               >  1:5545241/1‑93 (MQ=255)
          ggggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGc                                                                               >  1:968869/1‑93 (MQ=255)
          ggggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGATGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGc                                                                               >  2:4479418/1‑93 (MQ=255)
           gggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCg                                                                              <  1:2058137/93‑1 (MQ=255)
           gggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCg                                                                              <  1:544806/93‑1 (MQ=255)
           gggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCg                                                                              <  2:798275/93‑1 (MQ=255)
            ggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGa                                                                             <  2:3906163/93‑1 (MQ=255)
             gagaCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAg                                                                            <  1:4738253/93‑1 (MQ=255)
              agaCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAgg                                                                           >  2:4091767/1‑93 (MQ=255)
               gaCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGt                                                                          <  1:2283621/93‑1 (MQ=255)
                aCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGtt                                                                         >  1:1275297/1‑93 (MQ=255)
                aCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGtt                                                                         <  2:4249569/93‑1 (MQ=255)
                  aaTGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTc                                                                       >  1:4943501/1‑93 (MQ=255)
                  aaTGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTc                                                                       >  2:3751773/1‑93 (MQ=255)
                  aaTGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTc                                                                       <  2:4680926/93‑1 (MQ=255)
                   aTGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCg                                                                      >  1:4891816/1‑93 (MQ=255)
                    tGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGt                                                                     <  2:2852704/93‑1 (MQ=255)
                    tGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGCTTCGt                                                                     >  2:3950619/1‑93 (MQ=255)
                       aaaaTTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAg                                                                  <  1:2741607/93‑1 (MQ=255)
                        aaaTTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGt                                                                 <  1:1449919/93‑1 (MQ=255)
                         aaTTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTc                                                                >  2:2223237/1‑93 (MQ=255)
                          aTTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCg                                                               >  2:674738/1‑93 (MQ=255)
                            ttttCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGcc                                                             <  2:2361673/93‑1 (MQ=255)
                            ttttCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGcc                                                             <  1:6423328/93‑1 (MQ=255)
                             tttCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCt                                                            <  1:1937149/93‑1 (MQ=255)
                              ttCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTg                                                           <  2:6590416/93‑1 (MQ=255)
                               tCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGc                                                          <  2:1207933/93‑1 (MQ=255)
                                 aGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCgg                                                        >  2:4613070/1‑93 (MQ=255)
                                 aGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCgg                                                        <  2:2441589/93‑1 (MQ=255)
                                  gCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGc                                                       <  1:1536252/93‑1 (MQ=255)
                                   cAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCt                                                      <  2:2788128/93‑1 (MQ=255)
                                     aGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTgg                                                    >  2:4979821/1‑93 (MQ=255)
                                      gggAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGt                                                   <  2:5704540/93‑1 (MQ=255)
                                       tgAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTa                                                  >  2:6316968/2‑93 (MQ=255)
                                        gAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTaa                                                 >  1:1114309/1‑93 (MQ=255)
                                         aaaaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAc                                                <  2:1539037/93‑1 (MQ=255)
                                          aaaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAGCCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAcc                                               <  1:38211/93‑1 (MQ=255)
                                             ttGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCa                                            >  1:3568470/1‑93 (MQ=255)
                                              tGAGGGGTTGATCACGTTTTGTACAGAAGTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGGTTCGTCAGTCGCCTGCGCCTGATAAACACaa                                           >  1:6462499/1‑93 (MQ=25)
                                               gAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCaaa                                          <  1:4172669/93‑1 (MQ=255)
                                                aGGGGTTGATCACGTTTTGTTCTGAATTGCAGATAACAAAAAGCCCCTCCGGAGCGAGGTTTCGTCAGTCGCTTGCGGCTGGGAACCGCTAAg                                         >  2:6729286/1‑93 (MQ=37)
                                                aGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAg                                         >  1:4491441/1‑93 (MQ=255)
                                                 ggggTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGc                                        >  2:723478/1‑93 (MQ=255)
                                                 ggggTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGc                                        >  1:1655253/1‑93 (MQ=255)
                                                   ggTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGcac                                      <  1:2593355/93‑1 (MQ=255)
                                                      tGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACt                                   <  1:5601051/93‑1 (MQ=255)
                                                       gATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTg                                  >  1:3182692/1‑93 (MQ=255)
                                                       gATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTg                                  >  1:4420441/1‑93 (MQ=255)
                                                         tCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGta                                >  1:607559/1‑93 (MQ=255)
                                                            cGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGtatta                             >  2:2327774/1‑93 (MQ=255)
                                                            cGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGtatta                             >  1:3326061/1‑93 (MQ=255)
                                                             gTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGtattat                            >  1:2451779/1‑93 (MQ=255)
                                                              ttttGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATg                           <  1:720829/93‑1 (MQ=255)
                                                              ttttGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATg                           >  1:4329901/1‑93 (MQ=255)
                                                               tttGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGt                          <  2:2888560/93‑1 (MQ=255)
                                                               tttGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGt                          >  1:6328537/1‑93 (MQ=255)
                                                                  gTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCaa                       <  1:3445103/93‑1 (MQ=255)
                                                                  gTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCaa                       >  1:550969/1‑93 (MQ=255)
                                                                   tACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAac                      <  2:1711447/93‑1 (MQ=255)
                                                                    aCTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAaca                     >  2:256840/1‑93 (MQ=255)
                                                                       gAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTg                  >  1:5168060/1‑93 (MQ=255)
                                                                        aaTTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGa                 >  2:907229/1‑93 (MQ=255)
                                                                         aTTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGaa                >  2:513407/1‑93 (MQ=255)
                                                                         aTTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGaa                <  1:3093403/93‑1 (MQ=255)
                                                                          ttGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGaaa               <  2:1913718/93‑1 (MQ=255)
                                                                          ttGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGaaa               >  2:6186124/1‑93 (MQ=255)
                                                                           tGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAg              <  2:645937/93‑1 (MQ=255)
                                                                            gCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGt             <  1:2235527/93‑1 (MQ=255)
                                                                            gCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGt             >  1:6225952/1‑93 (MQ=255)
                                                                             cAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGta            <  2:1267578/93‑1 (MQ=255)
                                                                             cAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGta            <  1:2968450/93‑1 (MQ=255)
                                                                                aTAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATAc         <  2:3081503/93‑1 (MQ=255)
                                                                                aTAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATAc         <  2:5667883/93‑1 (MQ=255)
                                                                                aTAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATAc         <  2:4518164/93‑1 (MQ=255)
                                                                                  aaCAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACgt       <  2:5626688/93‑1 (MQ=255)
                                                                                  aaCAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACgt       >  1:1585741/1‑93 (MQ=255)
                                                                                    cAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACgtgt     <  1:4099996/93‑1 (MQ=255)
                                                                                    cAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACgtgt     >  1:3209161/1‑93 (MQ=255)
                                                                                      aaaagCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTc   <  1:671386/93‑1 (MQ=255)
                                                                                      aaaagCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTc   >  1:5820527/1‑93 (MQ=255)
                                                                                      aaaagCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTc   <  1:3794758/93‑1 (MQ=255)
                                                                                       aaagCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTcc  <  2:5879980/93‑1 (MQ=255)
                                                                                          |                                                                                         
CAGCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTCC  >  NC_000913/3850660‑3850839

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