Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,851,510 T→C intergenic (‑895/‑33) 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,851,5100TC99.7% 977.3 / ‑4.9 298intergenic (‑895/‑33)ivbL/emrDilvB operon leader peptide/2‑module integral membrane pump
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (121/176);  total (121/177)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

Only 100 of 298 total aligned reads displayed.

GCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTG  >  NC_000913/3851420‑3851599
                                                                                          |                                                                                         
gCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGa                                                                                         <  2:5598724/93‑1 (MQ=255)
 ccGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGAc                                                                                        <  2:6565360/93‑1 (MQ=255)
 ccGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGAc                                                                                        <  1:5516056/93‑1 (MQ=255)
  cGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACt                                                                                       <  2:2455694/93‑1 (MQ=255)
  cGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACt                                                                                       <  1:5616142/93‑1 (MQ=255)
  cGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACt                                                                                       <  2:3884607/93‑1 (MQ=255)
    cTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTa                                                                                     <  1:1067680/93‑1 (MQ=255)
     tAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTAt                                                                                    >  1:1334513/1‑93 (MQ=255)
     tAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTAt                                                                                    <  2:5980913/93‑1 (MQ=255)
     tAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTCATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTAt                                                                                    <  2:6018058/93‑1 (MQ=255)
      aaaCAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTAtt                                                                                   <  1:5195148/93‑1 (MQ=255)
       aaCAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTc                                                                                  <  1:1576769/93‑1 (MQ=255)
       aaCAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTc                                                                                  >  2:3431351/1‑93 (MQ=255)
       aaCAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTc                                                                                  <  1:4507441/93‑1 (MQ=255)
        aCAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCt                                                                                 >  1:1618579/1‑93 (MQ=255)
         cAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTa                                                                                <  2:1294863/93‑1 (MQ=255)
          aaTTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTaa                                                                               <  1:5007590/93‑1 (MQ=255)
          aaTTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTaa                                                                               >  2:1143118/1‑93 (MQ=255)
          aaTTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTaa                                                                               <  1:6232033/93‑1 (MQ=255)
            ttCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAatt                                                                             <  1:6712284/93‑1 (MQ=255)
             tCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAatta                                                                            <  1:335970/93‑1 (MQ=255)
             tCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAatta                                                                            >  2:6634612/1‑93 (MQ=255)
              ccccATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAattat                                                                           <  1:3546275/93‑1 (MQ=255)
              ccccATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAattat                                                                           >  2:4796570/1‑93 (MQ=255)
              ccccATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAattat                                                                           <  1:5291061/93‑1 (MQ=255)
               cccATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAattatt                                                                          >  2:4744899/1‑93 (MQ=255)
                 cATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTAtttt                                                                        >  2:2784011/1‑93 (MQ=255)
                  aTTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTAttttt                                                                       >  1:5056857/1‑93 (MQ=255)
                  aTTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTAttttt                                                                       <  2:2599356/93‑1 (MQ=255)
                   ttCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGTCTTATTCTAATTATTTTta                                                                      <  1:2928198/93‑1 (MQ=255)
                   ttCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTta                                                                      >  1:5973346/1‑93 (MQ=255)
                   ttCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTta                                                                      <  1:3290076/93‑1 (MQ=255)
                   ttCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTta                                                                      <  2:5514198/93‑1 (MQ=255)
                   ttCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTta                                                                      >  2:5913512/1‑93 (MQ=255)
                    tCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTtat                                                                     <  2:2162690/93‑1 (MQ=255)
                     ccTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTtata                                                                    <  1:6476494/93‑1 (MQ=255)
                     ccTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTtata                                                                    <  1:3642107/93‑1 (MQ=255)
                      cTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATaa                                                                   <  2:6256650/93‑1 (MQ=255)
                       tGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATaaa                                                                  >  1:4180113/1‑93 (MQ=255)
                       tGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATaaa                                                                  >  2:3242663/1‑93 (MQ=255)
                        ggCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATaaaa                                                                 <  1:1785326/93‑1 (MQ=255)
                        ggCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATaaaa                                                                 >  2:4372424/1‑93 (MQ=255)
                        ggCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATaaaa                                                                 >  1:2716840/1‑93 (MQ=255)
                         gCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAg                                                                <  1:1948084/93‑1 (MQ=255)
                         gCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAg                                                                >  2:6171524/1‑93 (MQ=255)
                          cGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGc                                                               >  2:529057/1‑93 (MQ=255)
                          cGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGc                                                               <  1:3139798/93‑1 (MQ=255)
                           gTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCa                                                              <  1:443923/93‑1 (MQ=255)
                            tataTCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCAt                                                             >  2:1600185/1‑93 (MQ=255)
                              tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATcc                                                           >  2:3522964/1‑93 (MQ=255)
                              tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATcc                                                           <  1:6246190/93‑1 (MQ=255)
                              tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATcc                                                           >  2:3779710/1‑93 (MQ=255)
                                tCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGt                                                         >  1:1004985/1‑93 (MQ=255)
                                     cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATaa                                                    <  2:1790767/93‑1 (MQ=255)
                                       aaCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATg                                                  <  2:798937/93‑1 (MQ=255)
                                       aaCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATg                                                  <  2:5851960/93‑1 (MQ=255)
                                       aaCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATg                                                  >  1:2762076/1‑93 (MQ=255)
                                        aCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGa                                                 <  2:642112/93‑1 (MQ=255)
                                          aTTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGaaa                                               >  1:262228/1‑93 (MQ=255)
                                           ttCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGaaaa                                              <  2:4218626/93‑1 (MQ=255)
                                            tcatcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAg                                             <  1:887158/93‑1 (MQ=255)
                                             catcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAgg                                            >  2:3190048/1‑93 (MQ=255)
                                               tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCa                                          >  2:1977941/1‑93 (MQ=255)
                                               tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCa                                          >  2:4937409/1‑93 (MQ=255)
                                                 aaTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCaaa                                        <  1:2211993/93‑1 (MQ=255)
                                                 aaTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCaaa                                        >  1:3051723/1‑93 (MQ=255)
                                                  aTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAg                                       <  2:3132819/93‑1 (MQ=255)
                                                   tgtgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGa                                      >  1:1215201/1‑93 (MQ=255)
                                                     tgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGaaa                                    <  1:5120903/93‑1 (MQ=255)
                                                     tgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGaaa                                    <  1:1398216/93‑1 (MQ=255)
                                                      gATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAAc                                   >  2:4029393/1‑93 (MQ=255)
                                                         aGATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTc                                >  2:4513240/1‑93 (MQ=255)
                                                          gATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCa                               >  1:1278867/1‑93 (MQ=255)
                                                          gATTCCTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCa                               >  1:2064539/1‑93 (MQ=255)
                                                              ccTCTCCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAttt                           >  2:3331190/1‑93 (MQ=255)
                                                                tctcCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGt                         >  2:1034537/1‑93 (MQ=255)
                                                                tctcCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGt                         >  1:3926991/1‑93 (MQ=255)
                                                                 ctcCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                        <  1:3786276/93‑1 (MQ=255)
                                                                 ctcCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                        <  2:2216240/93‑1 (MQ=255)
                                                                  tcCCGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTa                       >  1:4215537/1‑93 (MQ=255)
                                                                   cccGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAt                      >  1:2590242/1‑93 (MQ=255)
                                                                    ccGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAtt                     <  1:2285029/93‑1 (MQ=255)
                                                                    ccGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAtt                     >  2:5802733/1‑93 (MQ=255)
                                                                    ccGCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAtt                     <  1:6062685/93‑1 (MQ=255)
                                                                    ccGCATTTATGGGAATGCGTAGCGAATTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAtt                     >  1:487114/1‑93 (MQ=255)
                                                                      gCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                   <  1:5755441/93‑1 (MQ=255)
                                                                      gCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                   >  1:6645910/1‑93 (MQ=255)
                                                                      gCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                   <  1:6381021/93‑1 (MQ=255)
                                                                      gCATTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                   >  2:2832473/1‑93 (MQ=255)
                                                                        aTTTATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATg                 >  1:4587576/1‑93 (MQ=255)
                                                                          ttATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt               >  1:1676627/1‑93 (MQ=255)
                                                                          ttATGGGAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt               >  2:391949/1‑93 (MQ=255)
                                                                              gggAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTa           >  2:3046868/1‑93 (MQ=255)
                                                                               ggAATGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAACTTGTTATTGATGTTGGTAt          <  1:6570227/93‑1 (MQ=255)
                                                                                   tGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACt      >  1:3672666/1‑93 (MQ=255)
                                                                                   tGCGTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACt      >  2:3706400/1‑93 (MQ=255)
                                                                                      gTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGt   <  2:5150219/93‑1 (MQ=255)
                                                                                      gTAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGt   >  2:4584299/1‑93 (MQ=255)
                                                                                       tAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTg  <  1:3463178/93‑1 (MQ=255)
                                                                                       tAGCGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTg  >  2:6725147/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTG  >  NC_000913/3851420‑3851599

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