Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,408 A→T 13.1% 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,4080AT13.1% 205.9 / 20.4 84intergenic (+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 (35/38);  new base T (7/4);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

TATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTT  >  NC_000913/1196259‑1196557
                                                                                                                                                     |                                                                                                                                                     
tATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                                                                       >  1:1699146/1‑150 (MQ=255)
   tctaGGCGCAATCACCGCGAAACCCGTAACATATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACg                                                                                                                                                    <  2:103606/147‑1 (MQ=21)
   ggAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACg                                                                                                                                                    <  2:202887/150‑1 (MQ=255)
    gAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGAGCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGc                                                                                                                                                   >  2:1109721/1‑150 (MQ=255)
      aGGCGCAATCAACGCGAACACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTACACTGCTGAAAGGGTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTGAAAGTTATTAACGGGAGCGTAACgaga                                                                                                                                                 >  1:103606/1‑147 (MQ=21)
         cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCg                                                                                                                                              <  1:671364/150‑1 (MQ=255)
           cAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTATTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCttt                                                                                                                                            >  1:146120/1‑150 (MQ=255)
           cAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtt                                                                                                                                            <  1:1714704/150‑1 (MQ=255)
           cAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgtt                                                                                                                                            <  2:1295731/150‑1 (MQ=255)
             aTCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaa                                                                                                                                          <  2:1063187/150‑3 (MQ=255)
              tCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaat                                                                                                                                         <  2:1039514/150‑4 (MQ=255)
              tCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATAGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgttgtt                                                                                                                                         <  1:98686/150‑1 (MQ=255)
               cAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaata                                                                                                                                        <  1:396995/150‑5 (MQ=255)
                aaCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGtttt                                                                                                                                       >  2:1235153/1‑150 (MQ=255)
                 aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCTCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGttttt                                                                                                                                      <  1:326759/150‑1 (MQ=255)
                   gcaaaaaCCGTACACTATGCCTTCGAGCGGCTGATGGCTGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTg                                                                                                                                    <  2:180242/147‑1 (MQ=37)
                   ccgAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTg                                                                                                                                    <  1:1445069/149‑1 (MQ=255)
                    cgAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGt                                                                                                                                   >  1:1028608/1‑150 (MQ=255)
                         aCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgc                                                                                                                              <  2:1699146/150‑1 (MQ=255)
                                 cTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACg                                                                                                                      <  1:424435/150‑1 (MQ=255)
                                  tATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACgg                                                                                                                     >  2:8008/1‑150 (MQ=255)
                                   aTGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGt                                                                                                                    <  1:1169556/150‑1 (MQ=255)
                                   aTGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGt                                                                                                                    >  1:375676/1‑150 (MQ=255)
                                              cGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAttt                                                                                                         >  2:1418386/1‑150 (MQ=255)
                                                tCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGATAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTa                                                                                                       <  2:485831/150‑1 (MQ=255)
                                                  tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATc                                                                                                     <  1:1501254/150‑1 (MQ=255)
                                                  tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATc                                                                                                     >  2:957158/1‑150 (MQ=255)
                                                   gatggatgGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTGGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCa                                                                                                    >  2:479720/1‑150 (MQ=255)
                                                      ggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaa                                                                                                 <  2:964058/150‑1 (MQ=255)
                                                           gcgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGt                                                                                            <  2:1298867/150‑1 (MQ=255)
                                                               tAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATc                                                                                        <  2:547401/150‑1 (MQ=255)
                                                               tAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATc                                                                                        <  2:746980/150‑1 (MQ=255)
                                                                 aaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCaa                                                                                      <  2:492018/150‑1 (MQ=255)
                                                                     tctgAAATTTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTATGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAcc                                                                                  <  1:1264973/149‑1 (MQ=255)
                                                                      ctgAAATGTTCAGAGTTTGGTGACGCAATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAccc                                                                                 >  2:368187/1‑150 (MQ=255)
                                                                        gAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCt                                                                               >  1:852584/1‑150 (MQ=255)
                                                                        gAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCt                                                                               >  1:661307/1‑150 (MQ=255)
                                                                         aaaTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTc                                                                              <  1:735914/150‑1 (MQ=255)
                                                                         aaaTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTc                                                                              <  1:553040/150‑1 (MQ=255)
                                                                              ttCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAg                                                                         >  2:1409957/1‑150 (MQ=255)
                                                                               tCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAgc                                                                        >  2:1513222/1‑149 (MQ=255)
                                                                                  tagTTTGGTGACGCGATCATCGAAAAAATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCTTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGtttt                                                                     <  1:1314813/149‑1 (MQ=255)
                                                                                    gTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTATTTTTgg                                                                   >  1:1069553/1‑150 (MQ=255)
                                                                                     tttGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg                                                                  <  1:1695631/150‑1 (MQ=255)
                                                                                     tttGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGAGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTTGTTTTggg                                                                  <  1:1019822/150‑1 (MQ=255)
                                                                                         gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg                                                              >  2:1313188/1‑150 (MQ=255)
                                                                                          tGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAgg                                                             <  1:1140117/150‑1 (MQ=255)
                                                                                           gACGCGATCATCGAAAACATGTAATGCAGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGATTTTTGTTAGGATGATAACGGTTATCAAAATTTTATCAAAAAAAGATATCAAAACCCCTCGGTATTTTTGTGGTAgtc                                                            >  1:611818/1‑148 (MQ=255)
                                                                                            aCGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCt                                                           <  1:608866/150‑1 (MQ=255)
                                                                                             cgcgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAACCCCCTCGGTAGTTTTGGGGTAGGCTg                                                          >  2:329333/1‑150 (MQ=255)
                                                                                             cgcgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTg                                                          >  1:754894/1‑150 (MQ=255)
                                                                                             cgcgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTg                                                          >  2:583730/1‑150 (MQ=255)
                                                                                               cgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGc                                                        >  1:294053/1‑150 (MQ=255)
                                                                                                 atcatcGAAAACATGTAATGCCGTTGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCg                                                      <  1:479720/150‑1 (MQ=255)
                                                                                                 atcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCg                                                      >  2:433408/1‑150 (MQ=255)
                                                                                                  tcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg                                                     <  1:8008/150‑1 (MQ=255)
                                                                                                     tcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCa                                                  <  1:1235153/150‑1 (MQ=255)
                                                                                                      cGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAg                                                 >  1:845671/1‑150 (MQ=255)
                                                                                                         aaaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGCGCGTACCGCTCCCGTGGTTATTTGTTAGGCTGCTAACGGTTATCAAAATTTTAACAAAAAAACTAATCAAAACCCATCGGTAGTTTTGGGGTTGGCTGGCCGGTCAggtg                                              >  1:1418588/1‑150 (MQ=255)
                                                                                                          accATCTAATGCCGTAGTTTGTTAATTTTATTACCGGGAGCGTAACGCTCCAGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAACTTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtgg                                             <  2:1608469/148‑1 (MQ=255)
                                                                                                           aCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtggt                                            <  2:294053/150‑1 (MQ=255)
                                                                                                             aagTAATGCCGTAGTTTGTAAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAg                                          <  2:46140/148‑1 (MQ=255)
                                                                                                              tGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGt                                         <  2:845671/150‑1 (MQ=255)
                                                                                                                 aaTGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTct                                      <  2:1028608/150‑1 (MQ=255)
                                                                                                                   tGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTctac                                    <  1:1012159/150‑1 (MQ=255)
                                                                                                                       gTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTctactact                                <  2:1664251/150‑1 (MQ=255)
                                                                                                                         aGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAg                              <  1:957158/150‑1 (MQ=255)
                                                                                                                              gTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTcc                         >  1:942921/1‑150 (MQ=255)
                                                                                                                               ttAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTccc                        <  1:433408/150‑1 (MQ=255)
                                                                                                                                tAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGATGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCca                       >  2:1006971/1‑150 (MQ=255)
                                                                                                                                tAAATTTATTAACGGGAGCGTAACGCTCACGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCca                       >  1:1542609/1‑150 (MQ=255)
                                                                                                                                  aaTTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCcaca                     >  1:221747/1‑150 (MQ=255)
                                                                                                                                   aaaaaacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACAt                    >  1:1544177/8‑150 (MQ=255)
                                                                                                                                   aTTTATTAACGGGAGCGTAACGCTCCCGTTTTTTTTTGGAAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGGTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTCCTAGGCTCCCACAt                    >  2:821213/1‑150 (MQ=255)
                                                                                                                                      aaacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTTGTAGTTCTACTACTAGTCTCCCACATAGa                 >  2:916078/5‑150 (MQ=255)
                                                                                                                                      aaacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGa                 >  2:892746/5‑150 (MQ=255)
                                                                                                                                      aaacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATAAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGa                 >  1:1350497/5‑150 (MQ=255)
                                                                                                                                       attaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGATCTACTACTAGTCTCCCACATAGat                >  2:1207167/1‑150 (MQ=255)
                                                                                                                                         taACGGGAGCGTAACGCTCCCTTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGatat              >  1:330463/1‑150 (MQ=255)
                                                                                                                                         taACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGatat              >  1:691934/1‑150 (MQ=255)
                                                                                                                                         cgaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGatat              >  2:1392647/3‑150 (MQ=255)
                                                                                                                                          aaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATAtt             >  1:1553929/1‑150 (MQ=255)
                                                                                                                                          aaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATAtt             >  1:109644/1‑150 (MQ=255)
                                                                                                                                                aGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAg       >  1:747141/1‑150 (MQ=255)
                                                                                                                                                  cGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCt     >  2:795236/1‑150 (MQ=255)
                                                                                                                                                     aaCGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCtttt  <  1:1313188/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGATATTCCTTAGCTTTT  >  NC_000913/1196259‑1196557

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: