Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 74,084 A→G H146H (CAT→CAC thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091374,0840AG98.0% 165.0 / NA 51H146H (CAT→CACthiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (29/21);  minor base T (0/1);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCC  >  NC_000913/73854‑74324
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
gCCTGAGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCtt                                                                                                                                                                                                                               <  1:201884/250‑1 (MQ=255)
           gccgcccAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCACAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACTTGCGCCAGCAAAATGCCTGGCAGTGCGTAgg                                                                                                                                                                                                                   <  2:230807/247‑1 (MQ=255)
                   gCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAg                                                                                                                                                                                                           <  1:328184/251‑1 (MQ=255)
                        aGCACGGTGGAGAAGCTGGCGAAACAGAGCATAACGATAAGGGCAGCAACCGGCGGGATTTGTCGCCGGAACCACGGCAATTCGACGAAGCGGAAAAAATGCCAGCCGCGCCTCCCAAGCTGGGCGGCAAGCTGACGCTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGTTTAAAACATACGTGCGCCAGCAAAATACCTTGCAGACCGTAGGGCGAAAAGGTccc                                                                                                                                                                                                      >  1:304956/1‑250 (MQ=255)
                          cACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACt                                                                                                                                                                                                    <  1:166398/251‑1 (MQ=255)
                                ggCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGAc                                                                                                                                                                                              >  1:286446/1‑251 (MQ=255)
                                cgCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGAc                                                                                                                                                                                              <  2:307526/250‑1 (MQ=255)
                                         ggCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGAtt                                                                                                                                                                                     <  1:352094/251‑1 (MQ=255)
                                         ggCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGAAACCACGGCCCTTCGACGAAGCGGACAAACTGCCAGCCGCGCCTCCCAAGCTGGGCGGCAAGCTGCCGTTGTTCGCCGGGAATGTTTTCCCGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGAtt                                                                                                                                                                                     >  2:248104/1‑251 (MQ=255)
                                             aaaCAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCa                                                                                                                                                                                 <  1:141999/251‑1 (MQ=255)
                                              aaCAGAGCATAACGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGTCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGGTTTCCAGTGCCTGGAGTAATAAGCGGCCCGCCATCGGGAGATTAAAAAAAACGTGGCCCAGCAAAAACCCTTGCAGGCCGAAGGGCGAAAAGGGCCACTCCAGACCGAGCGATTGGc                                                                                                                                                                                  >  2:82385/1‑249 (MQ=255)
                                                           gATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgccagcc                                                                                                                                                                   >  1:43436/1‑251 (MQ=255)
                                                            aTAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGAGGAAAAAATGCCAGCCCCGGCTCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgccagcc                                                                                                                                                                   >  1:233400/1‑250 (MQ=255)
                                                             tAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgccagcc                                                                                                                                                                   >  1:97443/1‑249 (MQ=255)
                                                                                 aTTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGcc                                                                                                                                             >  2:270901/1‑250 (MQ=255)
                                                                                  tgtGTCGCAGTAGCCACGGCCATTATACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCACGCTGACGGTGTTCGCCGGGAATGTTGTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAAGACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCTAAAAGGGCCACTCCAGTCCGAGCGATTGGCAGAGTGTTGCCAGGCATCCCTTGCGTCCATAGACGCt                                                                                                                                             <  2:366765/248‑1 (MQ=255)
                                                                                         ccGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGc                                                                                                                                     >  1:248932/1‑251 (MQ=255)
                                                                                            tAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCg                                                                                                                                   >  1:318418/1‑250 (MQ=255)
                                                                                               ccACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGaaa                                                                                                                                >  1:358858/1‑250 (MQ=255)
                                                                                                   ggCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGc                                                                                                                           >  1:267932/1‑251 (MQ=255)
                                                                                                         tCGACGAAGCGGAAAAAATGCCAGCCGCGCATTCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCa                                                                                                                      <  1:370315/250‑1 (MQ=255)
                                                                                                          cgacgaAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTGGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATCCCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAACGCCGAAAACAGCAACCagc                                                                                                                    <  2:233400/251‑3 (MQ=255)
                                                                                                              gaAGAGGAAAAAATGCCAGGCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTAGTTCGCCGGGAATGATTTCCAGTGCATGGAGTAATAAGCGGCTAGCCATCGGGAGAATAAAAAAAACGTGCGCCAGCAAAATACCTAGCAGGCCGTAGGGCGAAAAGGTCGACTCCAGACCGAGCGATTGGCAGAGTGTTGACAGACAGCCACGGAGGACAGAGACGCACAGAATGCCAAAAACAGCAACCAacaacg                                                                                                                >  2:298185/1‑248 (MQ=255)
                                                                                                                aGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCAGCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGcagcaagccccggg                                                                                                              <  2:336400/251‑14 (MQ=255)
                                                                                                                                 ccgcGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGc                                                                                             <  1:248104/248‑1 (MQ=255)
                                                                                                                                   gcGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGcgc                                                                                           >  1:57033/2‑249 (MQ=255)
                                                                                                                                      cATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGcgcac                                                                                         >  1:172024/1‑250 (MQ=255)
                                                                                                                                                tGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAAc                                                                                >  1:109588/1‑249 (MQ=255)
                                                                                                                                                 gggCGGCTAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGTAGGATTAAGTTCATTTCGCACAGACGCAACa                                                                               <  2:366450/249‑1 (MQ=255)
                                                                                                                                                 gggCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCTAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACa                                                                               <  1:374539/249‑1 (MQ=255)
                                                                                                                                                       cAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGc                                                                       <  2:344261/251‑1 (MQ=255)
                                                                                                                                                              acgttGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCg                                                                <  2:120003/247‑1 (MQ=255)
                                                                                                                                                               cgttGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCgg                                                               <  1:256900/248‑1 (MQ=255)
                                                                                                                                                               cgttGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCgg                                                               >  2:77573/4‑251 (MQ=255)
                                                                                                                                                                 ttGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGa                                                              >  2:21201/2‑250 (MQ=255)
                                                                                                                                                                   gtgcgcagtgaaTGTTTTCCAGTGACTGGAGGAATAAGCGGCTCGCCACCGGTAGATTAAAAAATACGTGCGCCAGAAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCAGTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAg                                                           <  2:144475/240‑1 (MQ=255)
                                                                                                                                                                   gTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAg                                                           <  1:270901/251‑1 (MQ=255)
                                                                                                                                                                         cGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTg                                                     >  1:17789/1‑251 (MQ=255)
                                                                                                                                                                          gggaaTGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGc                                                    >  2:332832/5‑251 (MQ=255)
                                                                                                                                                                          gggaaTGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGACGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGc                                                    >  2:332815/5‑251 (MQ=255)
                                                                                                                                                                           ggaaTGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCg                                                   >  2:132625/4‑251 (MQ=255)
                                                                                                                                                                                             gggAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGGCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAg                                 <  1:168526/250‑1 (MQ=255)
                                                                                                                                                                                                 gTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGact                             >  2:1829/1‑249 (MQ=255)
                                                                                                                                                                                                     aaaGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTCCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATcg                         <  2:331689/250‑1 (MQ=255)
                                                                                                                                                                                                        tcGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGCGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAACGCGCCTGCGATAGAGCGCGCGGGCGCGGAATATCGCgg                      >  2:54168/2‑251 (MQ=255)
                                                                                                                                                                                                                 cATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTCGGGCGCAAAGGTCCACTCCAGACCGCGCGCTTGGCCGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGCATGCCGACAACCGCAACCAGCACCGGGAGGATTCAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGAc               >  1:176052/1‑249 (MQ=255)
                                                                                                                                                                                                                       taGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAACACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGCCGCAACAGCGCCAGCCCCCCCGGAAAGCGCCTGCGATCGAGCGCGCGGGCGAGGAATATCGCGGGTCCGCCAGAGAGcc       >  2:339722/2‑250 (MQ=255)
                                                                                                                                                                                                                       taGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTg                                                                                              <  2:313960/163‑1 (MQ=255)
                                                                                                                                                                                                                       taGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTg                                                                                              >  1:313960/2‑164 (MQ=255)
                                                                                                                                                                                                                        aGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGc        >  1:189904/1‑249 (MQ=255)
                                                                                                                                                                                                                            tAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCCGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTACGACAGAGAGCAGCGcc  >  2:97185/1‑251 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
GCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCC  >  NC_000913/73854‑74324

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