Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 89,699 T→C Y22Y (TAT→TAC mraZ → DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091389,6990TC98.2% 199.4 / ‑6.0 57Y22Y (TAT→TACmraZDNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (33/23);  minor base A (1/0);  total (34/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCG  >  NC_000913/89454‑89940
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
cATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCggg                                                                                                                                                                                                                                               >  1:380029/1‑250 (MQ=255)
 aTAACGTTAATTTGCTGCGCAAGTTGGACGTAAAATAAAAAACGCTGAGATTAGCCGTAAACATTGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTGGACAAGATTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTGCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGaa                                                                                                                                                                                                                                             <  2:421785/251‑1 (MQ=255)
   aaCGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAAc                                                                                                                                                                                                                                            <  1:207538/250‑1 (MQ=255)
      gTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAgct                                                                                                                                                                                                                                        >  1:467289/1‑251 (MQ=255)
       ttAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGGGGGAATAAGGGGTGAGGATGGCCAGGTACCGGGAGCAACGGTAGGCAATCTCGCAAGCAACGGGCGGTTATCCATGTCTACCCGTTACCGGGAACAgcgg                                                                                                                                                                                                                                       >  2:222934/1‑249 (MQ=255)
          aTTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCtt                                                                                                                                                                                                                                    >  1:448887/1‑251 (MQ=255)
                 gcGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAATCTTTTCCTTCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAAc                                                                                                                                                                                                                              <  1:90189/249‑1 (MQ=255)
                        ttGGACGTAAAATAAACAACGCTGATATTAGCCGTAACCATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCCAAGTGGGAATAAGGGCTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCACCCCGTTACCGGGCAAAGCGGCTTGGGAACGCTGCCgg                                                                                                                                                                                                                      >  2:28254/1‑251 (MQ=255)
                         tGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAACTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGCGAACCGCTGCTTGAGAACCCTGCCGGt                                                                                                                                                                                                                     >  2:390942/1‑251 (MQ=255)
                                aaaaTAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATg                                                                                                                                                                                                               <  1:361702/250‑1 (MQ=255)
                                     aaaCAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTg                                                                                                                                                                                                          <  1:331827/250‑1 (MQ=255)
                                               gATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGCATAAGGGGTGAGGCTGGCATGTTCCGGGGGGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGAc                                                                                                                                                                                                >  1:542389/1‑250 (MQ=255)
                                                 tatTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTTTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAGTTGTGAGGCAAGGTGGGAATAAGGGGTTAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCACTCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTCCCGGGCACAGATGCTTGAGACAGCTGCCAGTCAAATGGCTTGCCCCATTGACAtt                                                                                                                                                                                             >  2:89352/1‑251 (MQ=255)
                                                  atTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACAtt                                                                                                                                                                                             <  1:79696/250‑1 (MQ=255)
                                                       ccGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCAc                                                                                                                                                                                       >  1:484140/1‑251 (MQ=255)
                                                         gTAAACATCGTGTTGTTTACCTCGGTAGGCCTTGGGACTGGCTTGCTAAGCGTTTCCTCCGCTCCGTAAACTCCGTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGTAAAGTGCGCTTATCAGGGCCTACCCGTTTCCGGGAACAGCTGTTTGAGAACGCCGCAGGTCAAAGGGTTTGCACCATTGACATTCATAAccc                                                                                                                                                                                     <  2:31424/251‑1 (MQ=255)
                                                         gTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCAccc                                                                                                                                                                                     <  1:525989/251‑1 (MQ=255)
                                                             aCATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCCACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGGGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATccaccgttc                                                                                                                                                                                 >  2:316991/1‑243 (MQ=255)
                                                                 cGGGTTTTTTACCTCGGTATGCCTTGAGACTGGCTTGACAAGCTTTTCCTCCGCACCGTAAACTCCCTTCCGTGGGCAATTGTGGGCCAAAGAGAGAACAAGGGGTGACGCTGGCATGTTCCGGCGAGCAACGGTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTCCCTACCCGTTACCGGTAACAGCTGCTTGCGAACGCTGCCCGTCCAATGTATTGCACCCTTGACATTCCTCACCCGTGCctgc                                                                                                                                                                             >  1:552882/1‑251 (MQ=255)
                                                                             cTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATATCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATGCATCACCCGTGCCTGCTGCTTTAccccc                                                                                                                                                                 >  1:213481/1‑251 (MQ=255)
                                                                                ggTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTTCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCcctgc                                                                                                                                                              >  2:283320/1‑251 (MQ=255)
                                                                                                ggCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTAt                                                                                                                                              <  2:448887/251‑1 (MQ=255)
                                                                                                    tGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCg                                                                                                                                            <  1:91977/249‑1 (MQ=255)
                                                                                                    tGACAAGCTTTTCCTCCGATCCGTAAACTAACTTCCGTGGGACATTGTGGCGCACAGTGGGAATACGGGGACAGGCTGGCATGTTCCGGGGAGCAACGCTAGTCACTCTCGCCAGCAAAGGGCGCTTATCAGTGCCTACCCGTTAACGGGAACAGCTGCTTGAGAACGCTTCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGGTTGGGAAATTATCgtg                                                                                                                                          >  2:246295/1‑249 (MQ=255)
                                                                                                           cTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCACATGGTTTGCa                                                                                                                                                                                                        >  1:298039/1‑182 (MQ=255)
                                                                                                           cTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCa                                                                                                                                                                                                        <  2:298039/182‑1 (MQ=255)
                                                                                                                          gTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCgg                                                                                                                    >  1:32889/1‑250 (MQ=255)
                                                                                                                           tAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCg                                                                                                                     >  1:564597/1‑249 (MQ=255)
                                                                                                                              aCTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCAt                                                                                                                >  1:375311/1‑251 (MQ=255)
                                                                                                                                tCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATg                                                                                                               >  1:295708/1‑250 (MQ=255)
                                                                                                                                         tGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATc                                                                                                                                                                                         >  1:150239/1‑167 (MQ=255)
                                                                                                                                         tGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATc                                                                                                                                                                                         <  2:150239/167‑1 (MQ=255)
                                                                                                                                            gAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGATAGTCAATCTCGACAGCAAAGGGCGCTTAGCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAgc                                                                                                  <  2:127325/251‑1 (MQ=255)
                                                                                                                                              aatggtgGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGGCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGt                                                                                                                       <  2:499859/224‑1 (MQ=255)
                                                                                                                                              aaTTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGt                                                                                                                       >  1:499859/1‑228 (MQ=255)
                                                                                                                                                          aaGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAg                                                                                                                                          >  1:277591/1‑197 (MQ=255)
                                                                                                                                                          aaGTGGGAAGAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACAGTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAg                                                                                                                                          <  2:277591/197‑1 (MQ=255)
                                                                                                                                                            gggggAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGGCAAACTCGACAGCAAAGGGCGCTGATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGTTTGAGCGCCGTGTGCAGCGCCt                                                                                  <  2:378418/249‑1 (MQ=255)
                                                                                                                                                            gagggAATCAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGCCAAATGGTTTGAACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCt                                                                                  >  2:496423/3‑251 (MQ=255)
                                                                                                                                                               ggAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACt                                                                               <  1:395884/251‑1 (MQ=255)
                                                                                                                                                                gAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTg                                                                              >  1:313094/1‑251 (MQ=255)
                                                                                                                                                                gAATAAGGGGCGAAGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCACAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGACCGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAACAATTATCACGTCTTTCGAGCATGAACCCGGTTGTGCGCCGTGTGCAGCGCCTACTg                                                                              >  2:258258/1‑251 (MQ=255)
                                                                                                                                                                        ggTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAAAAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGGGCAGCGCCTACTGTTAGGTc                                                                       >  2:33501/1‑250 (MQ=255)
                                                                                                                                                                              gCTGGCATGTGCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAg                                                                <  2:564597/251‑1 (MQ=255)
                                                                                                                                                                                tGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCg                                                              >  1:422440/1‑251 (MQ=255)
                                                                                                                                                                                      tGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGt                                                         >  1:207329/1‑251 (MQ=255)
                                                                                                                                                                                               gAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatggat                                                >  1:112201/1‑251 (MQ=255)
                                                                                                                                                                                                         aGTCAAGCGCGACAGCAAAGGGCGCTTATCAGTGCCTTCCCGTTACAGGGAACAGCTGCTTTAGAACGCTGCCGGGCAAATGGTTTGCACCATTGACATTCAGCACCCGTGCCTGCTGCTTTACCCCCTGCCTTAATGGGAAATTAGCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTTTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATTTCGCAGGTc                                      <  2:88548/251‑1 (MQ=255)
                                                                                                                                                                                                          gccAATCTCGACAGCAAAGGGCGCTTAACTGTGTCTACCCTGTTCCTGGAACAGCTGCTTGAGAACGCTTCCGGTCAAATGGATTGCACCATTGACATTCATCACCCGTGCCTGCTGGTTTACCCCCTGCCGGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCg                                     <  2:479848/249‑1 (MQ=255)
                                                                                                                                                                                                            cAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAACAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTATTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGAt                                   >  2:233097/1‑251 (MQ=255)
                                                                                                                                                                                                                ctcGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGCGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGTTGGCGCCGGTCGATTGtt                               >  2:426496/1‑251 (MQ=255)
                                                                                                                                                                                                                    aCAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTACTGCTTTACCCCCTGCCTGAATGGAAAATTATCGCGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATc                           >  1:460132/1‑251 (MQ=255)
                                                                                                                                                                                                                                gcTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTg               <  1:390942/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ttATCAGTGCCTACCCGTTACCGGGACCAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCTTGCCTGAATGGGAAATTATCGAGCACAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGATTGTCAGATGGATGGCGCAGGTCGATTGGTAATCGCGCAAGTACTGCg             >  2:425987/1‑251 (MQ=255)
                                                                                                                                                                                                                                         tGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACAt      >  1:330600/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ccccGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGcc   <  2:380029/248‑1 (MQ=255)
                                                                                                                                                                                                                                              aCCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCg  <  2:482500/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCG  >  NC_000913/89454‑89940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: