Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 73,413 A→G V370A (GTG→GCG)  thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,4130AG89.1% 116.9 / ‑2.6 46V370A (GTG→GCG) thiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (36/5);  minor base C (1/4);  total (37/9)
Fisher's exact test for biased strand distribution p-value = 3.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCACCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGC  >  NC_000913/73168‑73653
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
cTTTTAAGCGCGACCCGCCTTCAATGCCCAGCGACTGACACAAAATACTGTAGAGGGCGGTGATATCGCGCATCGGGTTTTCAAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATACCGTCAACAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACCATCACCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTg                                                                                                                                                                                                                                               >  1:182538/1‑250 (MQ=255)
                      aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCcgcg                                                                                                                                                                                                                       >  1:336400/5‑251 (MQ=255)
                      aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGAGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCGCATCACCAGCGCCTGACCC‑GCCAGCCTTTTCTGCcgcg                                                                                                                                                                                                                       >  1:356498/5‑251 (MQ=255)
                            cAGCGACTGAACCAACATACTGGAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGCCAATCACAATGCCGTCCGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGGAGCCAGCACACTCCCCGGCAAGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCCGACCC‑GCCAGCCTTTTCTGCCGCGCCCgac                                                                                                                                                                                                                 >  2:240693/1‑249 (MQ=255)
                                gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTc                                                                                                                                                                                                             >  1:262059/1‑251 (MQ=255)
                                                       ggcggTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGcatc                                                                                                                                                                                       >  2:124723/1‑250 (MQ=255)
                                                         cggTGATATCGCGCGTCGGGTTCTCCAGTACTTTCAGAGCATAAGGGATCGCCATTAAAGCATTGGTGAAAATCAAACTGCCGTCAGCATATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGACAGCACAATCCCCGGCATGGCGAGGATACACCA‑GCCGCTCATCTCCAGAGCCTGACCC‑GCCAGCATTTTCTTCCGCGCCCGAAGTTCAAGACTGCGCCATCGCCGCATCATg                                                                                                                                                                                    >  2:243883/1‑251 (MQ=255)
                                                            tGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATACCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATACCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTc                                                                                                                                                                                 >  2:185244/1‑251 (MQ=255)
                                                              atatCGCCCATCGGGTTTGCGAGTACTTGCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGACAGCAGATTGCGGCAGCCCAAGGG‑TGTTGTTGACAAGAAAAAAAAAGCCGTTAGCCAGCACAACCCCCGGCTTGGCGAGGAT‑CAACATGCCGCGCATCTCCAGCCCCTGACCC‑GCCATCATTTGCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg                                                                                                                                                                               <  1:32135/251‑1 (MQ=255)
                                                                     gcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTAcc                                                                                                                                                                        >  2:230326/1‑250 (MQ=255)
                                                                        tCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGcc                                                                                                                                                                                                                          >  1:109276/1‑198 (MQ=255)
                                                                         cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTACCGCATTGGTGAACATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGCT‑AAACATGCCGCTCCTCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCCGCAGAAGTTCGCGACTGCTCCCTAGCAGCATCATGGACAGCACAACGccc                                                                                                                                                                     >  2:176587/1‑248 (MQ=255)
                                                                          gggTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGcaca                                                                                                                                                                    >  1:191797/1‑250 (MQ=255)
                                                                              tttCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑Gcc                                                                                                                                                                                                                                       <  1:23811/179‑1 (MQ=255)
                                                                              tttCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑Gcc                                                                                                                                                                                                                                       >  2:23811/1‑179 (MQ=255)
                                                                                      tcTTTAAGCGCAAAGGGGACCGCCATTAACGCATTGGTAAAAAACACAATGCCGTCAGCAGATTGCGGCAGGC‑AATGGGTTTTGTTGAGCAGAAAAAAGAAGCCGGAAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTgccg                                                                                                                                                        <  1:195088/249‑1 (MQ=255)
                                                                                               gcATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAAt                                                                                                                                              >  2:347080/1‑251 (MQ=255)
                                                                                                cATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCCGCGCAAt                                                                                                                                              >  2:65863/1‑250 (MQ=255)
                                                                                                   aaGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGACAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGCCCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGCCTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCCGCGCAATACGc                                                                                                                                          >  2:22228/1‑251 (MQ=255)
                                                                                                               ttAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCcac                                                                                                                              >  2:358699/1‑251 (MQ=255)
                                                                                                                 aaCGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCACCGAGGTCcac                                                                                                                              >  2:214889/1‑249 (MQ=255)
                                                                                                                    gCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAgcgc                                                                                                                         >  2:168279/1‑251 (MQ=255)
                                                                                                                       ttGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCt                                                                                                                                                                                                                                                <  2:5255/130‑1 (MQ=255)
                                                                                                                       ttGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCt                                                                                                                                                                                                                                                >  1:5255/1‑130 (MQ=255)
                                                                                                                        tGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGc                                                                                                                     >  1:40136/1‑251 (MQ=255)
                                                                                                                             aaaaTCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑GGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCCTAGCAGCATCATGGTCCGCACTACACACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCcacac                                                                                                                >  2:139360/1‑250 (MQ=255)
                                                                                                                                aTCACAATGCCGTCAGCAGATTGCGGCAGGCCCATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCCATACGCAACGAGGTCCACAGCGCCTGCCACAGCAc                                                                                                             >  2:116534/1‑251 (MQ=255)
                                                                                                                                       tGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTtgt                                                                                                      >  1:194401/1‑251 (MQ=255)
                                                                                                                                              agcagATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGc                                                                                                >  2:352094/1‑250 (MQ=255)
                                                                                                                                                 agATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGCCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGCGCCAGCACtt                                                                                            >  2:266532/1‑251 (MQ=255)
                                                                                                                                                   aTTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTcc                                                                                          >  1:242189/1‑251 (MQ=255)
                                                                                                                                                           aGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTg                                                                                  >  1:366765/1‑251 (MQ=255)
                                                                                                                                                             gccaatgg‑tgtTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGc                                                                                >  1:344261/8‑251 (MQ=255)
                                                                                                                                                                  tgg‑tgtTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCAGGCATGGCGAGCAT‑CAACATGCCGATCAACTCCAGCCCCTGACCC‑GCCAGCCTTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCAGGTGTGCCAGCACTTCCGGCAACAGGCGAttt                                                                           >  1:158223/3‑251 (MQ=255)
                                                                                                                                                                   gg‑tgtTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGAttt                                                                           >  1:128752/2‑250 (MQ=255)
                                                                                                                                                                       gttgttGAGCAAAAAAAAGACGCCGGTAGCCAGCACAACCCCCGGAAGGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCCCCTGCCCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACTCCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTagcc                                                                       <  1:64040/251‑4 (MQ=255)
                                                                                                                                                                               gaaataaaaaaaggcgggtgaccAGCAAAACCCCCGGAGTGGCGAGAAC‑CAACAGCCCCCTCATCCCCAGCCCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATcgacg                                                               <  2:365809/230‑1 (MQ=255)
                                                                                                                                                                                cAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATcgacg                                                               >  1:230807/1‑250 (MQ=255)
                                                                                                                                                                                 agaaaaaaGAAGCCGGTAGCCACCACAACCCCCGGCAGGGCGAGGAC‑AAACATGCCCCTCATCCCCAGCGCCTGCCCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGAt                                                             <  1:309716/248‑1 (MQ=255)
                                                                                                                                                                                  gTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATc                                                            >  1:72593/1‑251 (MQ=255)
                                                                                                                                                                                         aaagCCGGTCGCCAGCACAACCCCCGGCAGGGCGGGGAC‑AAACATCCCGCTCATCCCCAGCCCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCa                                                     <  1:374877/250‑1 (MQ=255)
                                                                                                                                                                                               ggTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACg                                               >  2:98132/1‑251 (MQ=255)
                                                                                                                                                                                                               cccGGCATGGCGCGGAA‑CAACCTGCCGCTCATCACACGCGCCTGACCA‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGACAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCAACGGTTGTGCCAGCACTTCCGGCACCTGGCGATTTAGCCAATCGACGATCACCGCCAGCAACGGTGGCAGCAAcagc                                 >  2:374539/1‑249 (MQ=255)
                                                                                                                                                                                                                 cGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAAcagcagc                              >  1:293835/1‑250 (MQ=255)
                                                                                                                                                                                                                        ggggagaaaaaaCATGCCCCTCACCCCCAGCGCC‑GACCCCGCCAGCATTTCCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCAGGGTCAGCACTACGCACAATACCCCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACGGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCa                       <  1:124723/241‑1 (MQ=255)
                                                                                                                                                                                                                                             cATCTCCAGCGCCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGc  >  1:169790/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCACCTGACCC‑GCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGC  >  NC_000913/73168‑73653

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