Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,408 A→T 17.5% intergenic (+35/+459) icd → / ← ymfD isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,4080AT17.5% 409.0 / 76.0 171intergenic (+35/+459)icd/ymfDisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (68/73);  new base T (15/15);  total (83/88)
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.59e-01

Only 100 of 178 total aligned reads displayed.

TGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGG  >  NC_000913/1196309‑1196506
                                                                                                   |                                                                                                  
tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                    >  2:1069505/1‑100 (MQ=255)
tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                    >  2:1867634/1‑100 (MQ=255)
tGATGGATGGCGCAAAACTGCTGACATGTTCAGAGTTTAGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                    >  1:2725075/1‑100 (MQ=255)
 gatggatgGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                    >  1:2589938/1‑99 (MQ=255)
  atggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAAc                                                                                                  <  2:1172151/100‑1 (MQ=255)
  atggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAAc                                                                                                  <  1:209481/100‑1 (MQ=255)
  atggagggCGCGAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACAGGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                    <  1:2836821/98‑1 (MQ=255)
   tggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACg                                                                                                 >  1:4393884/1‑100 (MQ=255)
    ggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGc                                                                                                <  2:1569914/100‑1 (MQ=255)
    ggagggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTTACGCGATCATCGAAAACATGTAATGCCGTATTTTGTTAAATTTATTAACGGTAGCGTAACg                                                                                                 <  1:1989123/99‑1 (MQ=25)
       tggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTcc                                                                                             >  1:3405392/1‑100 (MQ=255)
        ggCGCTAAACTGCGGAAATGTTCAGAGGTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTcc                                                                                             <  1:4164808/99‑1 (MQ=255)
         gcgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCg                                                                                           <  1:359443/100‑1 (MQ=255)
         gcgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCg                                                                                           <  2:4154047/100‑1 (MQ=255)
          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTTACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgt                                                                                          <  1:3280780/100‑1 (MQ=255)
          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgt                                                                                          <  2:167900/100‑1 (MQ=255)
          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgt                                                                                          >  2:2058770/1‑100 (MQ=255)
           gcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtt                                                                                         <  1:1052129/100‑1 (MQ=255)
           gcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtt                                                                                         <  1:1143681/100‑1 (MQ=255)
           gcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtt                                                                                         <  2:4005849/100‑1 (MQ=255)
           gcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAACTTATTAACGGGAGCGTTACGCTCCCgtt                                                                                         <  1:2045246/100‑1 (MQ=255)
            cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtta                                                                                        <  2:762591/100‑2 (MQ=255)
               aaCTGCTGAAATGTTCAGAGTTTGGTTACGCGATCATCGAAAACATGTACTGCCGTAGTTTGTTATATTTATTAACGGGAGCGTAACGCTCCCGTTGttt                                                                                     <  1:921803/100‑1 (MQ=255)
                 ctgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGttttt                                                                                   <  1:2783302/100‑1 (MQ=255)
                 ctgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGttttt                                                                                   >  1:1114031/1‑100 (MQ=255)
                  tgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGCAACGCTCCCGTTGtttttt                                                                                  <  2:1040333/100‑1 (MQ=255)
                    ctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGt                                                                                >  2:3245452/1‑100 (MQ=255)
                    ctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGt                                                                                >  1:2420807/1‑100 (MQ=255)
                           gTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAGCGGGAGCGTAACGCGCCCGTTGTTTTTTGTTAGgctg                                                                         >  1:2963285/1‑100 (MQ=255)
                           gTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctg                                                                         >  2:3887882/1‑100 (MQ=255)
                           gTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctg                                                                         >  2:3847559/1‑100 (MQ=255)
                           gTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctg                                                                         >  2:1412298/1‑100 (MQ=255)
                            ttCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctgc                                                                        >  2:688737/1‑100 (MQ=255)
                              cAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTa                                                                      >  2:949184/1‑100 (MQ=255)
                                   tttGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGt                                                                 <  2:198137/100‑1 (MQ=255)
                                       gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATc                                                             <  2:490297/100‑1 (MQ=255)
                                        tGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCa                                                            >  2:1296513/1‑100 (MQ=255)
                                           cgcgTTCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCaaaa                                                         >  1:2448849/1‑100 (MQ=255)
                                           cgcgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCaaa                                                          >  1:2850916/1‑99 (MQ=255)
                                             cgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTATGCTGCTAACGGTTATCAAAAtt                                                       >  2:4286126/1‑100 (MQ=255)
                                             cgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAtt                                                       >  1:2063961/1‑100 (MQ=255)
                                              gATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAttt                                                      >  1:102792/1‑100 (MQ=255)
                                                tcatcGAAAACATGTAATTCCGTAGGTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTa                                                    <  1:3205049/100‑1 (MQ=255)
                                                tcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTa                                                    <  1:4242080/100‑1 (MQ=255)
                                                tcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTa                                                    <  1:986594/100‑1 (MQ=255)
                                                 catcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCGCGTTGTTTTTTGTTGGGCTGCTAACGGTTATCAAAATTTTAt                                                   >  2:3301771/1‑100 (MQ=255)
                                                 catcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTAt                                                   <  2:2443501/100‑1 (MQ=255)
                                                 catcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTAt                                                   >  2:590887/1‑100 (MQ=255)
                                                 catcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTAt                                                   <  2:799864/100‑1 (MQ=255)
                                                 catcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTAt                                                   <  1:2489883/100‑1 (MQ=255)
                                                  atcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATc                                                  >  1:2129481/1‑100 (MQ=255)
                                                    cGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTTTTTTTTTGTTAACGGTTATCAAAATTTTATaaa                                                >  1:2673300/1‑100 (MQ=255)
                                                        aaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTCGGCTGCTAACGGTTATCAAAATTTTATCaaaaaa                                            >  1:2314798/1‑100 (MQ=255)
                                                        aaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaaaa                                            >  2:3914327/1‑100 (MQ=255)
                                                           aTGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGt                                         <  2:2129489/100‑1 (MQ=255)
                                                             gTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTa                                       >  2:2617877/1‑100 (MQ=255)
                                                                   ccGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCaaaa                                 >  2:3708726/1‑100 (MQ=255)
                                                                   ccGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCaaaa                                 >  2:2858463/1‑100 (MQ=255)
                                                                    cGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCaaagc                                <  2:2982772/100‑3 (MQ=255)
                                                                    cGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAc                                <  2:3753423/100‑1 (MQ=255)
                                                                    cGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTGTTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAc                                <  1:1077192/100‑1 (MQ=255)
                                                                          ttGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCg                          <  1:304152/100‑1 (MQ=255)
                                                                              tcaaTTTATTAACGGGAGCGTAACGCTCCACTTCTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAg                      <  1:3244545/98‑1 (MQ=255)
                                                                              tAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAg                      <  1:3268307/100‑1 (MQ=255)
                                                                               acaaaaaacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGt                     >  1:1655559/10‑100 (MQ=38)
                                                                                caaaaaacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGtt                    >  2:1855967/9‑100 (MQ=38)
                                                                                   ttattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTg                 >  1:1455651/1‑100 (MQ=255)
                                                                                    tattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCACAAACAGTTATCAAAACCCCTCGGTAGtttt                  >  1:1689704/1‑98 (MQ=255)
                                                                                    tattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgg                <  1:3719744/100‑1 (MQ=255)
                                                                                    nattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgg                >  2:2480502/2‑100 (MQ=255)
                                                                                    aaacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgg                >  1:1674947/5‑100 (MQ=255)
                                                                                     attaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg               <  2:1655882/100‑1 (MQ=255)
                                                                                     aacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg               >  2:688527/4‑100 (MQ=255)
                                                                                      ttaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTtggg              >  1:2609047/1‑100 (MQ=255)
                                                                                      ttaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgggg              <  2:3405410/100‑1 (MQ=255)
                                                                                      ttaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgggg              <  2:1437865/100‑1 (MQ=255)
                                                                                      ttaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg               <  2:2695947/99‑1 (MQ=255)
                                                                                       taccGGGAGCGTAACGCACCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAATATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgcgg              <  1:1125742/96‑1 (MQ=255)
                                                                                       taACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGt             >  2:3775674/1‑100 (MQ=255)
                                                                                       taACGGGAGCGTAACGCTCACGTGGTTTGTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgggg              <  1:506417/99‑1 (MQ=255)
                                                                                       caaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGt             >  2:1326032/2‑100 (MQ=255)
                                                                                       caaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGt             >  1:4013324/2‑100 (MQ=255)
                                                                                        aaCGGGAGCGTTACGCTCCCGATGTTTTTTGTTATGCTGCCAACGGTTAACACACTTTTATCAAAAAAAGTTATCAAAACCCCTCGATAGTATTAGGGTa            >  1:3816717/1‑100 (MQ=255)
                                                                                        aaCGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATTAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTa            >  1:4359955/1‑100 (MQ=255)
                                                                                         aCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg           >  1:3558568/1‑100 (MQ=255)
                                                                                         aCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg           >  1:2504546/1‑100 (MQ=255)
                                                                                          cGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTaag          >  1:1260873/1‑98 (MQ=255)
                                                                                          cGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAgg          >  2:3518901/1‑100 (MQ=255)
                                                                                          cGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAgg          <  2:180337/100‑1 (MQ=255)
                                                                                          cGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAgg          >  2:203780/1‑100 (MQ=255)
                                                                                           gggAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGc         <  2:3360447/100‑1 (MQ=255)
                                                                                           gggAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGc         <  2:4399445/100‑1 (MQ=255)
                                                                                             gAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGGTTTGGGGTAGGCTg       >  2:2731315/1‑100 (MQ=255)
                                                                                               gCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAACAACCTCGGAAGTTTTGGGGTAGGCTGGc     >  2:2563114/1‑100 (MQ=255)
                                                                                                cGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGcc    <  1:1857822/100‑1 (MQ=255)
                                                                                                cGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGcc    <  2:3821650/100‑1 (MQ=255)
                                                                                                cGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAAAAAAAAAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGcc    <  1:580728/100‑1 (MQ=255)
                                                                                                cGTAACGCTCCCGTTGATTTTTGTTAGGCTACTAACGGTTATCACAATTATATCAAAAAAAGTTATCAAAACCCCCCGGTAGTTTTGGGGTAGGCTg       >  1:723024/1‑97 (MQ=255)
                                                                                                 gTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCg   <  2:4226863/100‑1 (MQ=255)
                                                                                                  tAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg  <  1:3903226/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGG  >  NC_000913/1196309‑1196506

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