Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,408 A→T 20.7% 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,4080AT20.7% 168.0 / 41.4 82intergenic (+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 (31/34);  new base T (8/9);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

TGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCG  >  NC_000913/1196309‑1196505
                                                                                                   |                                                                                                 
tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTACTGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                   >  2:322061/1‑100 (MQ=255)
tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                   <  2:1286229/100‑1 (MQ=255)
 gatggatgGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTaa                                                                                                  <  1:262018/100‑1 (MQ=255)
         gcgcTAAACTGCTGAGATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCg                                                                                          <  1:1889414/100‑1 (MQ=255)
         gagcTAGACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGACAACATCTAATGCCGTACTCTGTTAAATTTATTAACGGAAGCGTAACGCTcc                                                                                            >  1:1145778/3‑98 (MQ=25)
          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgt                                                                                         <  2:2428080/100‑1 (MQ=255)
          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgt                                                                                         <  2:339136/100‑1 (MQ=255)
          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgt                                                                                         <  2:1859508/100‑1 (MQ=255)
          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGGGATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgt                                                                                         <  1:1286911/100‑1 (MQ=255)
           gcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGTGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtt                                                                                        <  1:1570560/100‑1 (MQ=255)
           gcTAAACTGCTGAAATGTTCAGAGTTTGGGGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTACATTTATTAACGGGAGCGTTACGCTCCCgtt                                                                                        <  1:1698837/100‑1 (MQ=25)
           gcTAAACTGCTGAAAGGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTGATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtt                                                                                        <  1:309090/100‑1 (MQ=255)
            cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtta                                                                                       <  2:2379496/100‑2 (MQ=255)
            cTAAACGGCTGAAATGTTCAGAGTTTGGTGACCCGAGCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCg                                                                                          <  1:2136486/97‑1 (MQ=255)
             tAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgttgt                                                                                      <  2:890687/100‑1 (MQ=255)
              aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgttgtt                                                                                     >  2:1919713/1‑100 (MQ=255)
               aaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaata                                                                                    >  2:2186766/1‑96 (MQ=255)
                aCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaataa                                                                                   <  1:1265540/100‑6 (MQ=255)
                aCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGtttt                                                                                   <  2:277993/100‑1 (MQ=255)
                 ctgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGttttt                                                                                  >  2:2196054/1‑100 (MQ=255)
                  tgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTGCCGTTGtttttt                                                                                 <  2:1705095/100‑1 (MQ=255)
                   gctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGtttttt                                                                                 >  1:1096496/1‑99 (MQ=255)
                   gctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTg                                                                                <  2:2362184/100‑1 (MQ=255)
                     tgAAATGTTCAGAGTTTGGTTACGTGATCATCGAAAACGTGTAATTCCCTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaataaattta                                                                              <  1:2545563/100‑11 (MQ=25)
                     tgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGtt                                                                              >  2:496919/1‑100 (MQ=255)
                     tgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGtt                                                                              >  2:874672/1‑100 (MQ=255)
                      gAAATGTTCAGTGTTTGGTGACACGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGTGTAACGCTCCCGTTGttttt                                                                                  <  1:2571852/95‑1 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTa                                                                             <  1:2527335/100‑1 (MQ=255)
                         aTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgc                                                                          <  2:1338346/100‑1 (MQ=255)
                           gTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctg                                                                        >  2:1988947/1‑100 (MQ=255)
                           gTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATCCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCTCGTTGTTTTTTGTTAGgctg                                                                        <  1:2145696/100‑1 (MQ=255)
                             tCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctgct                                                                      >  2:760549/1‑100 (MQ=255)
                              cAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTa                                                                     >  1:728738/1‑100 (MQ=255)
                               agagTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTaa                                                                    >  2:835153/1‑100 (MQ=255)
                               agagTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTaa                                                                    >  1:2199996/1‑100 (MQ=255)
                                  gTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACgg                                                                 <  2:472948/100‑1 (MQ=255)
                                    ttGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGtt                                                               <  1:1974211/100‑1 (MQ=255)
                                    ttGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGtt                                                               >  2:1070125/1‑100 (MQ=255)
                                       gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATc                                                            <  1:1761401/100‑1 (MQ=255)
                                             cgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAtt                                                      >  1:305995/1‑100 (MQ=255)
                                                 tatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTAt                                                  <  2:499905/99‑1 (MQ=255)
                                                 catcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCATAACGCTCTCGTTGTTATTTGTTAGACTGCTAACGGATATCAAATTTTGAt                                                  >  1:176069/1‑100 (MQ=255)
                                                  atcGAAAACATGTAATGCCGTAGTTTGTTAATTTTTTTAACGGGAGCGTAACGCTTCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATc                                                 <  1:1492723/100‑1 (MQ=255)
                                                  atcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTAta                                                 >  1:499903/1‑99 (MQ=255)
                                                  atcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATc                                                 <  1:2590895/100‑1 (MQ=255)
                                                   tcGAACACATGTAATGTCGTAGTTTGTTAATGTTTTTAACGGGAGCGTAACGCTCCAGTTGTTTTTTGTTAGGCTGATAACGGTTATCAAAATTGTAt                                                  <  1:321700/98‑1 (MQ=255)
                                                      aaaaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaa                                             >  2:616108/1‑100 (MQ=255)
                                                       aaaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaaa                                            <  1:168112/100‑1 (MQ=255)
                                                        aaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaaaa                                           <  1:899277/100‑1 (MQ=255)
                                                          cATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAg                                         >  1:458962/1‑100 (MQ=255)
                                                          cATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAg                                         <  1:1370972/100‑1 (MQ=255)
                                                            tGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGtt                                       >  1:2496230/1‑100 (MQ=255)
                                                             gTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTa                                      <  2:2185799/100‑1 (MQ=255)
                                                              tAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTAt                                     >  2:1667834/1‑100 (MQ=255)
                                                               aaTGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATc                                    >  2:1854427/1‑100 (MQ=255)
                                                               aaTGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATc                                    <  2:1935/100‑1 (MQ=255)
                                                                     gTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAcc                              >  1:2218662/1‑100 (MQ=255)
                                                                       aGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAcccc                            >  2:184789/1‑100 (MQ=255)
                                                                        gTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCt                           <  2:951821/100‑1 (MQ=255)
                                                                          ttGTTAAATTTATTAACGGGAGCGTAACGCTCACGTTGCTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCg                         >  1:1138741/1‑100 (MQ=255)
                                                                           tGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCATCg                         >  2:2574296/1‑99 (MQ=255)
                                                                           tGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAAATTATAAAAAAAAGTTATCAAAACCCCTCgg                        <  1:760546/100‑1 (MQ=255)
                                                                            gTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGt                       <  2:1936192/100‑1 (MQ=255)
                                                                               acaaaaaacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCATAAAAAGTTATCAAAACCCCTCGGTAg                     >  1:2329731/10‑99 (MQ=38)
                                                                                     attaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg              >  2:1217691/1‑100 (MQ=255)
                                                                                     attaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg              <  2:1216004/100‑1 (MQ=255)
                                                                                     attaACGGAAGGGTAACGCTCTCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg              <  2:1231060/100‑1 (MQ=255)
                                                                                     aatgaaggtaGCGCAACGCTCCAGTTGTTATTTGATAGGCTGCTACACGTTATCAAAAGTTTATGAAAAAAAGTTATCAAAACCCCGCGGGAGtttt                 <  1:24301/88‑1 (MQ=255)
                                                                                      ttaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgggg             >  2:1291960/1‑100 (MQ=255)
                                                                                      acaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgggg             >  2:921854/3‑100 (MQ=255)
                                                                                      acaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTg                >  1:2234621/3‑97 (MQ=255)
                                                                                       taACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGt            >  1:278946/1‑100 (MQ=255)
                                                                                        aaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTa           >  2:360474/1‑100 (MQ=255)
                                                                                        aaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTa           >  1:3445/1‑100 (MQ=255)
                                                                                         aCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg          >  2:861121/1‑100 (MQ=255)
                                                                                         aCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg          >  2:2105294/1‑100 (MQ=255)
                                                                                         aCGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg          >  2:2196150/1‑100 (MQ=255)
                                                                                           gggAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGc        <  2:458964/100‑1 (MQ=255)
                                                                                           gggAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGATATCAAAACCCCTCGGTAGTTTTGGGGTAGGc        <  1:235233/100‑1 (MQ=255)
                                                                                             gAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTg      >  2:248883/1‑100 (MQ=255)
                                                                                              aGCGTAACGCTCCCGTTGTTTTTTGTTTTTTTTCTAACGGTAATCAAAATTTTATCAAAAAACGTTATCAAAACCCCTCGGTAGTTTTCGGGTAGGCt       >  1:837518/1‑98 (MQ=255)
                                                                                               gCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGc    <  2:882432/100‑1 (MQ=255)
                                                                                                cGTAACGCTCCCGTTATTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg              <  1:1481330/89‑1 (MQ=255)
                                                                                                 gTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCg  <  2:1290445/100‑1 (MQ=255)
                                                                                                 gTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGGTATCAAAATCTTATCAAAAAAAGTTATCAAAACCCCCCGGTAGTTTTGGAGTAGGCTGGcc   >  1:1536764/1‑99 (MQ=255)
                                                                                                   |                                                                                                 
TGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCG  >  NC_000913/1196309‑1196505

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