Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 73,984 2 bp→CG coding (537‑538/1611 nt) thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,9840TC95.7% 147.1 / NA 47S180G (AGC→GGC) thiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/24);  minor base A (2/0);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_00091373,9850AG97.9% 133.3 / NA 47R179R (CGT→CGCthiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/24);  minor base C (1/0);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

ACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGG  >  NC_000913/73744‑74216
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                       
gcatcACCAACCCGAGGCAACACGCCAGCTGGATCAGCGCCAGCATCGCCGCGCGGGCATGTTCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGTTTGCCGGCGGACCTCCGCCCAGCTATAGCACGGTGGCGAAGCTGGCGAAACAGAGCGTAAAGATAAGTTCAGCAACCGGCGTGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCa                                                                                                                                                                                                                                <  1:243883/247‑1 (MQ=255)
gcAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCCGCATAGCAGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGAGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCCACGAAGCGGAAAAAATGCCAGCAGCGCATccc                                                                                                                                                                                                                                 >  1:154657/2‑250 (MQ=255)
    cACCAACCCGAGGCAACCCACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGAGGAAAAAATGCCCGCCGCGCATCCCCAGcc                                                                                                                                                                                                                            >  2:168526/1‑250 (MQ=255)
              aGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGGGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCAGCGCCTCCCCAGCTGggcggcc                                                                                                                                                                                                                    >  2:93826/1‑248 (MQ=255)
                  aacacaCCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCCGCCGCGCCTCCCAAGCTGGGCGGCAAgctg                                                                                                                                                                                                               >  2:328184/3‑247 (MQ=255)
                       acCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACgttg                                                                                                                                                                                                         >  1:70996/1‑248 (MQ=255)
                        ccATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGCAACTCAGCGCCGGATAGATTGCCAGCTCGATAGTGGTGGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGGAACCACGGCCATTCGACGAAGCGGAAAAAAGGCCGGCCGCGCATCCCAAGCTGGGCGGCAGGCTGACGTTGt                                                                                                                                                                                                        <  2:219680/251‑1 (MQ=255)
                              ggatcagcGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGcc                                                                                                                                                                                                   <  1:142916/246‑1 (MQ=255)
                                      gcCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGGTGTTCGCCGGGAATGtt                                                                                                                                                                                          >  2:166398/1‑251 (MQ=255)
                                       ccAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGttt                                                                                                                                                                                         <  2:264699/251‑1 (MQ=255)
                                          gCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCCTTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCCAGCTGACGATGTTCGCCGGGGATGTTTTcc                                                                                                                                                                                      >  2:318492/1‑251 (MQ=255)
                                            atcgCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAg                                                                                                                                                                                    <  2:293835/248‑1 (MQ=255)
                                              cgCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTg                                                                                                                                                                                  <  1:214889/250‑1 (MQ=255)
                                                                  aGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTc                                                                                                                                                              <  1:88793/251‑1 (MQ=255)
                                                                            cAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCTGCTCGCCATCggt                                                                                                                                                     <  1:65863/250‑2 (MQ=255)
                                                                            cAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCggta                                                                                                                                                    <  1:116534/251‑3 (MQ=255)
                                                                              gcgcCTGATAGATTGCCAGCTCGATGGTGGTTGCCTGCGGGCCGCCGCCCAGCGATAGCGCGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGGTGTTCGCCGGGAATGTTTTCCAGTGTGTGGAGTAATAAGCGGCTCgtcagctgtaga                                                                                                                                                  <  2:341432/251‑12 (MQ=255)
                                                                              gcgcCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGa                                                                                                                                                  <  1:22228/251‑1 (MQ=255)
                                                                                 ccTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGAtt                                                                                                                                                >  1:117806/1‑250 (MQ=255)
                                                                                           tGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTATATTaaaaaatacg                                                                                                                                      <  1:302536/250‑5 (MQ=255)
                                                                                           tGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAACGGAAAAAATGACAGCCGCGCATCCCCAGCTGGGAGGCAAGCTGACGTTGTTCTCCGGGAATGCTTTCCAGTGCCTGGAGTAATCAACGGCTCGCCCTCGGGAGATTaaaaaatacgg                                                                                                                                     >  2:2166/1‑246 (MQ=255)
                                                                                                              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)
                                                                                                                                                                          aTAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaatacgtgcgcca                                                                                                                               <  2:357881/178‑12 (MQ=255)
                                                                                                                                                                          aTAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaatacgtgcgcca                                                                                                                               >  1:357881/1‑167 (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)
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                       
ACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGG  >  NC_000913/73744‑74216

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