Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,769,474 G→A G134E (GGG→GAG)  ydiK → putative transporter YdiK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,769,4740GA93.8% 162.2 / 0.2 64G134E (GGG→GAG) ydiKputative transporter YdiK
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (34/26);  total (36/28)
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.49e-01

ATTGTTACGTTTGCAAAAGATCATGTTTGGCCGCCGCTCTCTCGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGC  >  NC_000913/1769226‑1769717
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
aTTGTTACGTTTGCAAAAGATCATGTTTGGCCGCCGCTCTCTCGGCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATggggg                                                                                                                                                                                                                                                   >  1:56378/1‑251 (MQ=255)
       cGTTTACAACAGATCAAGTTTGGCCGCCGCTATCTCGCCGTTCTGGTGAAGACGCTGTTTCTAGAGAAGGGATTTATCATCACTCGTGCTTTGCTGGTTAACAGTATCGTCGAAGACCGCGGCCCGCTCATTAAAGACATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCCGGCTGGCACAACTTGCTGGATATGGAGGGGACGGc                                                                                                                                                                                                                                            >  2:615992/1‑251 (MQ=255)
                aaGATCATGTTCGGCCGCCGCTCTCTCGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCCTCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCCCAACTTGCTGGATATGGAGGGGACGGCGATCATggg                                                                                                                                                                                                                                   >  2:98007/1‑250 (MQ=255)
                   aTCATGTTTGGCCGCCGCTCTCTCGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGaa                                                                                                                                                                                                                                <  1:167690/251‑1 (MQ=255)
                      aTGTTTGGCCGCCGCTCTCTCGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATg                                                                                                                                                                                                                                     <  1:391030/243‑1 (MQ=255)
                      aTGTTTGGCCGCCGCTCTCTCGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATg                                                                                                                                                                                                                                     >  2:391030/1‑243 (MQ=255)
                        gTTTGGCCGCCGCTCTCTCGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTcc                                                                                                                                                                                                                           <  1:196762/251‑1 (MQ=255)
                                gccgcTCTCTCGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGGCGGCGATCATGGCGAAAGTCCGCCCTtat                                                                                                                                                                                                                   >  2:242968/1‑251 (MQ=255)
                                acctctctctAGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCACTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTACAGCGGAGACATGACGTTACCCGATCTGGCGTGGCTTATTACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTTGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGACCGCCCtttt                                                                                                                                                                                                                   >  2:123550/5‑249 (MQ=255)
                                       tctcGCCGTTCTGGTGATGACGCTGTTATTAGTGCTGGTGTTTATCATCACTCTTGCATTGCTGGTTAACAGTATCGTCGCCGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCAGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGGCGGCGATCATGGCGCAAGTCCGCCCTTATATTGGc                                                                                                                                                                                                             >  2:600492/1‑250 (MQ=255)
                                        ctcGCAGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATACCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGACAGTCCGCCCTTATATTGGcac                                                                                                                                                                                                           >  2:83280/1‑251 (MQ=255)
                                         tcGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcacc                                                                                                                                                                                                          >  2:255238/1‑251 (MQ=255)
                                                   ggTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCAGTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacca                                                                                                                                                                                                   <  2:210216/248‑1 (MQ=255)
                                                       atgaCGGTGTTATTAGTGATGGTGTTTATCCTCCCTATTGCTTTGCTGGTTAACAGTATCGTCGATGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACACATTGCTGGATATGGAGGGGACGGCGGTCATGGCGAGAGAACGCACTTAGATTGGCACCACCACCCCCAGGtt                                                                                                                                                                                            >  2:337091/1‑251 (MQ=255)
                                                       atgaCGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGtt                                                                                                                                                                                            <  1:168972/251‑1 (MQ=255)
                                                                     aGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacca                                                                                                                                                                                                   >  1:623627/1‑230 (MQ=255)
                                                                     aGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccaccacca                                                                                                                                                                                                   <  2:623627/230‑1 (MQ=255)
                                                                           ggtatttatacacccTATTGCTTTGCTGGTTAACAGTGTCGTCGACGGAAGCGGGTCTCTAATTAAAGTCATTTCCAGCGGTGACATGTCGTTACCCGCTCTGGCGTTGCTTAATACCATTCCGGTGATTGGCGAGGAGCTGTATGCAGGCTGGCACAACTTGCTGGCTATGGAGGGGACGGCGATCATGGCGAAAGTACGCCCTTATATTGGCACCACCACCACCGGGTTCGTTGGGCGGGTGgcgcgt                                                                                                                                                                         <  2:386599/238‑3 (MQ=255)
                                                                           ggTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGGCGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGCAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGGAAGTCCGCCCTTATATTGGCAACACCACCACCTGGTTCGTTGGGCCGGCGgcgcct                                                                                                                                                                         >  2:546981/1‑248 (MQ=255)
                                                                               tttATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCAGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCAGCAGTGGCTTAATACCATTCCGCTGATTGGCGCGCAGCTGTATGCAAGCTGGCACAACTTGCTGGCTATGCAGGGGACGGCGACCATGGCGAAAGTCCGACATTACATTCGCCACACCCCCACCTGGTTCATTAGGCAGGCGCCGCTTATCgg                                                                                                                                                                    >  1:616170/1‑251 (MQ=255)
                                                                               tttATCATCACTATTGCTGTGCTGGTTAACAGTATCGTCGACTGCCGCGGCCAGCTAATTAACGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCACGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATCTGGAGGGGACGGCGATCATGGAGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCCATGGGCAGGCCGCGCATATCgg                                                                                                                                                                    >  2:195441/1‑251 (MQ=255)
                                                                                 tATCATCCCTATTGCTTTGCTGGTTAACCGTATCGTCGACGGCAGCGGCCCGCTATTTAAAGCCATTTACAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCAAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCAATCATGGCGAAATGCCGCCCTTATATTGGCCCCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgc                                                                                                                                                                  >  2:130073/1‑251 (MQ=255)
                                                                                 tATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCggg                                                                                                                                                                   <  2:251363/250‑1 (MQ=255)
                                                                                  atcatcCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgcg                                                                                                                                                                 >  1:276006/1‑251 (MQ=255)
                                                                                      tcCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCttt                                                                                                                                                             >  2:17660/1‑251 (MQ=255)
                                                                                       cccTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTGGGGCAGGCGGCGCATATCGGGCGCTTTa                                                                                                                                                            <  2:65149/251‑1 (MQ=255)
                                                                                           aTTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATg                                                                                                                                                          <  1:329094/249‑1 (MQ=255)
                                                                                           aTTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGt                                                                                                                                                        <  1:242968/251‑1 (MQ=255)
                                                                                                  ttctgtttAACAGTATCGTCGACGGCAGCGGCCCGTTAACTTACGCCTTTTCCAGCGGTGACATGACGTTACCCGAGATGGCGTGGCTTGATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCAGCCACAACTTGCTGGATATGGAGGGTACGGCGATCATGGCGAAATTCCGCCCTTATATTGGCACCACCACCACCTGGTGCGTTGGGCAGTCCGGGCATATCGGGCGCTTTATGGTGCATtg                                                                                                                                                  <  2:374087/244‑1 (MQ=255)
                                                                                                    ctGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTgc                                                                                                                                               >  1:17763/1‑251 (MQ=255)
                                                                                                             aGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACTACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATAGGACGCTTTATGGTGCATTGTGCGCTGATGCt                                                                                                                                      >  1:20262/1‑251 (MQ=255)
                                                                                                                  cggcgACGGCATCGGCCAGCTAATGAAAGCTATGTTCAGCGGGGTCATGTGGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGGGATTGGCGAGAAGCGGTATGCAGGCTGGCACAACGTGCTGGATATGGAGGGGACGGCGAGCATGGCGAAAGTTCGCCCTTATATCGTCTCCACCTCGATCTGGTTCGTTGGGCAGGCGGCGCATTTCGGGCGCTTTTTGGTGCATTGTGCGCTGATGCGTctct                                                                                                                                 <  2:323504/248‑1 (MQ=255)
                                                                                                                      gACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTc                                                                                                                               >  1:491532/1‑249 (MQ=255)
                                                                                                                          gCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGcc                                                                                                                         >  1:306811/1‑251 (MQ=255)
                                                                                                                          gCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTAACCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGcc                                                                                                                         >  1:171653/1‑251 (MQ=255)
                                                                                                                               ggCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCctgct                                                                                                                    <  1:92850/251‑1 (MQ=255)
                                                                                                                                  ccGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACTGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGt                                                                                                                  >  1:375654/1‑250 (MQ=255)
                                                                                                                                    gCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGtcct                                                                                                               >  2:3794/1‑248 (MQ=255)
                                                                                                                                         ttAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGCTTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGCACGGCGATCATGGCGAAAGTCCTCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCAGGCGCTTTTTGGTGCATTGTGCGCTGATGCTTCTCTTCATTGCACTGCTGTACTGgcgc                                                                                                          >  1:606407/1‑251 (MQ=255)
                                                                                                                                                        gCGGGGACATGACGTTGCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGca                                                                                           <  2:350481/251‑1 (MQ=255)
                                                                                                                                                          ggTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATCTGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTATCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGcac                                                                                          >  1:551075/1‑250 (MQ=255)
                                                                                                                                                          ggTGAAATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGACGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGTGCGCTTTATGGTGCATTGTGCGCTGTTGCTTCTCTTCAGTGCCCTGCTGGACTGGCGCGGGGAACCGGTGGcac                                                                                          >  2:444505/1‑250 (MQ=255)
                                                                                                                                                                  gACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTc                                                                                 >  2:52166/1‑251 (MQ=255)
                                                                                                                                                                    cGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGc                                                                               <  1:241482/251‑1 (MQ=255)
                                                                                                                                                                     gTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGcc                                                                              >  1:10715/1‑251 (MQ=255)
                                                                                                                                                                     gTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGcc                                                                              >  1:155016/1‑251 (MQ=255)
                                                                                                                                                                         cccGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTGCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATGCGCCAtt                                                                           <  1:83280/250‑1 (MQ=255)
                                                                                                                                                                         cccGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCAttt                                                                          >  1:116339/1‑251 (MQ=255)
                                                                                                                                                                                 ggCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAAccc                                                                  <  1:194827/251‑1 (MQ=255)
                                                                                                                                                                                 ggCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAAccc                                                                  <  1:533762/251‑1 (MQ=255)
                                                                                                                                                                                              aCCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGt                                                     <  1:98007/251‑1 (MQ=255)
                                                                                                                                                                                                cATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTc                                                   <  1:546981/251‑1 (MQ=255)
                                                                                                                                                                                                    ccGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGGCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGGACTGGCGCGGTGAACAGGAGGCCAAAGGCATTCGCCATTTTGCACACCGTCCGGCGGGCGTTCGCgg                                               >  2:14943/1‑251 (MQ=255)
                                                                                                                                                                                                      ggTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCgggg                                             >  2:422288/1‑249 (MQ=255)
                                                                                                                                                                                                         gATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCATGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATg                                          <  1:255506/251‑1 (MQ=255)
                                                                                                                                                                                                          aTTGGCGCGAAGCTGTATGCAGGCTGGCACAACGTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCTCCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGAGGCTTCGCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCAATCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATgc                                         <  1:445211/251‑1 (MQ=255)
                                                                                                                                                                                                            tGGCGCGAAGCTTTATGCAGGCTGGCACAACTTTCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTTCGCCCTGATATTGGCACCACCACCACCTGGTGCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCGGATGCGTCTCTTCAGTGCCCTGCTGTACTGGCGCGTTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTACGCGGAGATgcc                                        <  2:281492/250‑1 (MQ=255)
                                                                                                                                                                                                             ggCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCTCGGTGAACAGGTGGCACTAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATgccg                                       <  1:585566/250‑1 (MQ=255)
                                                                                                                                                                                                                     gCTGTATGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATAAAGGCGAAAGTCCGCCCTTATATTGGCACCCCCACCACCTGGATCGTTGGGCGGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGACCTACTGTACTGGCGCGGTGAACAGGTGGCACAAGGAATCCGCCATTTTGCAACCCGTCTGGCAGGCGTTCACGGAGATGCCGCCGtgctgc                              >  1:395719/1‑251 (MQ=255)
                                                                                                                                                                                                                          aTGCAGGCTGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTggcg                         >  1:473154/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tGGCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGCAAGCCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGc                 >  2:431952/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ggCACAACTTGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGc                 <  1:166783/250‑1 (MQ=255)
                                                                                                                                                                                                                                            tGCTGGATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCcgcg         <  1:94560/249‑1 (MQ=255)
                                                                                                                                                                                                                                            tGCTGGATATGGAGGGGACCGAGATCATGGCGAAAGTCCTCCCTTAGATTGGTACCCCCACCACCTGGTTCGTTGGGCAGTTGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCGGCTCGTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACATGTCATTCGCCATTTTGCGACCCGTCGGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGTTGGCGGCACAGGCTAGCcgcg         <  2:606407/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 gATATGGAGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGGCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGGGCGGTGgc  >  2:433880/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
ATTGTTACGTTTGCAAAAGATCATGTTTGGCCGCCGCTCTCTCGCCGTTCTGGTGATGACGCTGTTATTAGTGATGGTGTTTATCATCCCTATTGCTTTGCTGGTTAACAGTATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGC  >  NC_000913/1769226‑1769717

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