Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,925,158 T→C 6.0% N17N (AAT→AAC holE → DNA polymerase III, theta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,925,1580TC 6.0% 206.3 / 3.0 69N17N (AAT→AACholEDNA polymerase III, theta subunit
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base C (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCACCGTTTGGCTTCGGTCAATCTG  >  NC_000913/1925014‑1925305
                                                                                                                                                |                                                                                                                                                   
ttGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTc                                                                                                                                                  >  6:8317/1‑148 (MQ=255)
   tCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGAtt                                                                                                                                              <  6:124992/149‑1 (MQ=255)
     tAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTg                                                                                                                                            >  3:6996/1‑149 (MQ=255)
       aTGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATg                                                                                                                                                    <  2:387611/139‑1 (MQ=255)
       aTGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATg                                                                                                                                                    >  1:387611/1‑139 (MQ=255)
       aTGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGATCACTCTTAGTAACTTTTGCCCGCTAATGATGATGAGATTAAGAATGCTTAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGAttt                                                                                                                                             <  3:43852/146‑1 (MQ=255)
        tGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTggcg                                                                                                                                         <  1:234434/149‑1 (MQ=255)
          tCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGt                                                                                                                                                   <  7:377403/137‑1 (MQ=255)
          tCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGt                                                                                                                                                   >  8:377403/1‑137 (MQ=255)
                  gAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTggc                                                                                                                                          <  1:267358/138‑1 (MQ=255)
                  gAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTggc                                                                                                                                          >  2:267358/1‑138 (MQ=255)
                   aGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGt                                                                                                                              <  1:240211/149‑1 (MQ=255)
                        catcaATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAACGTCGATTTGGCGGCGGCCGGGGTGGCAt                                                                                                                         >  6:134250/1‑149 (MQ=255)
                          tcaATTCCGGCTTGCCATCCAGGCTAACTCTTAGTAAGTTTTGCCGGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCGACACTGGATCAAACAGAAATGGAGAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAtt                                                                                                                        >  6:130740/1‑148 (MQ=255)
                           caATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAttt                                                                                                                       >  1:167016/1‑148 (MQ=255)
                            aaTTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTaa                                                                                                                     <  3:62229/149‑1 (MQ=255)
                            aaTTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTaa                                                                                                                     <  8:83368/149‑1 (MQ=255)
                             aTTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTaaa                                                                                                                    <  2:167016/149‑1 (MQ=255)
                             aTTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTa                                                                                                                      >  3:354515/1‑147 (MQ=255)
                              ttCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGACGGCGGCCGGGGGGGCATTTAAAg                                                                                                                   <  3:184610/149‑1 (MQ=255)
                                 cGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAAc                                                                                                                <  7:676/149‑1 (MQ=255)
                                   gCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGc                                                                                                              >  7:233693/1‑149 (MQ=255)
                                        ccATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACaa                                                                                                         >  3:11969/1‑149 (MQ=255)
                                           tCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAatat                                                                                                      >  5:167442/1‑149 (MQ=255)
                                           tCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAatat                                                                                                      >  2:131087/1‑149 (MQ=255)
                                               ggCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATGAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGACGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGACg                                                                                                  >  5:20238/1‑149 (MQ=255)
                                                gCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCg                                                                                                  <  2:39099/148‑1 (MQ=255)
                                                    aCTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGAt                                                                                             <  1:131087/149‑1 (MQ=255)
                                                        ttAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTTAACTGGATCAAACAGAAATGGATAAAGTGAATGTCTATTTGGCGGCGGCAGGGGTGGCATTTACAGAACGCTACAAAATGCCGTTGATAGCt                                                                                         >  2:37405/1‑149 (MQ=255)
                                                        ttAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTAc                                                                                                           >  2:23635/1‑131 (MQ=255)
                                                        ttAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTAc                                                                                                           <  1:23635/131‑1 (MQ=255)
                                                        ttAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGATGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAAGGTCGATTTGTCGGCGGCAGGGGTGGCATTTAAATAACGCTACAATATGCCGGTGATctct                                                                                         >  2:6447/1‑149 (MQ=255)
                                                           gTAACTTTTACCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCCATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACTCTACAATATGCCGGTGATc                                                                                            >  7:93769/1‑143 (MQ=255)
                                                                ttttGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGt                                                                                 >  3:43774/1‑149 (MQ=255)
                                                                      ccGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAATGACCGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACg                                                                           >  8:346261/1‑149 (MQ=255)
                                                                      ccGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACg                                                                           >  1:314154/1‑149 (MQ=255)
                                                                      ccGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACAGGATCAAACAGAAATGGATAAAGTGAAAGTCGATTTGGCGGAGGCCGGGGTGGCAATTAAAGAACGCTACAAAATGCCGGTGATCGCTGAAGCGGTTGAACg                                                                           >  1:115833/1‑149 (MQ=255)
                                                                       cgcgAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGt                                                                          <  4:354515/149‑1 (MQ=255)
                                                                       cgcgAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGt                                                                          <  4:6996/149‑1 (MQ=255)
                                                                                tgaGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCt                                                                 <  4:43774/149‑1 (MQ=255)
                                                                                    gagaTTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAAc                                                             <  6:167442/149‑1 (MQ=255)
                                                                                      gaTTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACAt                                                           >  5:191844/1‑149 (MQ=255)
                                                                                      gaTTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAACGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACAt                                                           >  6:18659/1‑149 (MQ=255)
                                                                                       aTTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATATaaa                                                                                                                    >  6:131709/1‑91 (MQ=255)
                                                                                        ttAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACAttt                                                         <  2:314154/149‑1 (MQ=255)
                                                                                          aaGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAACGGCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTg                                                        <  5:134250/148‑1 (MQ=255)
                                                                                                  tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTc                                                                                                                                                  <  3:326304/50‑1 (MQ=255)
                                                                                                  tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTc                                                                                                                                                  >  4:326304/1‑50 (MQ=255)
                                                                                                      gaaTCTGGCTAATCTGGATGAAACAGAATTGGATAAAGTGATTGTGGATTTGGCGGCGGCCGGGGTGTCATTTAAAGAACGCTACTATTTGCCGGTGATCGCTGAAGCGGTTGAAGGTGAACTGCCTGATCATTTGCGCAGGTGGTTTc                                           <  2:115833/149‑1 (MQ=255)
                                                                                                      gaaTCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTc                                           <  5:8317/149‑1 (MQ=255)
                                                                                                             gCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGAGCAGCTGGTTTCGCGAGCg                                    >  4:624/1‑149 (MQ=255)
                                                                                                             gCTAAACTGGATCAAACAGAAATGGATAAAGTGAACGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGc                                     <  5:18659/148‑1 (MQ=255)
                                                                                                               tAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGc                                  <  3:12656/149‑1 (MQ=255)
                                                                                                                 aaCTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGc                                  >  4:12656/1‑147 (MQ=255)
                                                                                                                 aaCTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCtt                                <  6:191844/149‑1 (MQ=255)
                                                                                                                      gATCATACAGAAATGGATATAGTGAATGTCGATTTGGCGGGGGCCGGTGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTg                            <  6:20238/148‑1 (MQ=255)
                                                                                                                       aTCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCGGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGc                           >  7:95974/1‑148 (MQ=255)
                                                                                                                           aaCAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCAcc                      >  4:231737/1‑149 (MQ=255)
                                                                                                                           aaCAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCAcc                      >  3:199275/1‑149 (MQ=255)
                                                                                                                            aCAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCACCg                     >  1:160472/1‑149 (MQ=255)
                                                                                                                                    ggATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCACCGTTTGGCtt             >  7:97010/1‑149 (MQ=255)
                                                                                                                                     gATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGAGCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCACCGTTTGGCTTc            <  4:11969/149‑1 (MQ=255)
                                                                                                                                          aGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACAttt                                                         <  7:60443/99‑1 (MQ=255)
                                                                                                                                          aGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACAttt                                                         >  8:60443/1‑99 (MQ=255)
                                                                                                                                          aGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCACCGTTTGGCTTCGGTc        >  6:92021/1‑148 (MQ=255)
                                                                                                                                           gTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTAt                              >  2:356119/1‑125 (MQ=255)
                                                                                                                                           gTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTAt                              <  1:356119/125‑1 (MQ=255)
                                                                                                                                              aaTGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCACCGTTTGGCTTCGGTCAATc    <  8:97010/148‑1 (MQ=255)
                                                                                                                                               aTGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCACCGTTTGGCTTCGGTCAATCTg  >  3:320316/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TTGTCTAATGTCTGAGTCGAGGATCATCAATTCCGGCTTGCCATCCTGGCTCACTCTTAGTAACTTTTGCCCGCGAATGATGAGGAGATTAAGAATGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACAGCCTGAACATTTGCGCAGCTGGTTTCGCGAGCGGCTTATTGCCCACCGTTTGGCTTCGGTCAATCTG  >  NC_000913/1925014‑1925305

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: