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,4130AG97.5% 131.5 / NA 40V370A (GTG→GCG) thiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (16/23);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAG  >  NC_000913/73121‑73702
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                   
agCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGAACACTTTTAAGCGCGAACAGCCTTCAATGCCCAGCGACTGACACAACACACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGCTCGCCATTAACGCCTTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTAAcc                                                                                                                                                                                                                                                                                              >  2:369768/2‑300 (MQ=255)
                 gTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGcc                                                                                                                                                                                                                                                                             >  1:651372/1‑300 (MQ=255)
                                              acTTTTAAGCGCTACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGTTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATagc                                                                                                                                                                                                                                                <  2:268067/300‑1 (MQ=255)
                                              acTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGNTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCACGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATag                                                                                                                                                                                                                                                 <  1:495832/299‑1 (MQ=255)
                                              acTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATag                                                                                                                                                                                                                                                 <  1:675511/299‑1 (MQ=255)
                                              acTTTTAAGCGCGACCAGCCTTCAAGGCCCAGCGACTGACACAACATACCTGAGCGGGCGGTGATATCGCGCAGCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAGTGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAACCCCCGGCATGTCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCAACCGAAGTTCGCGACTGCTCCATagc                                                                                                                                                                                                                                                <  2:43719/300‑1 (MQ=255)
                                              acTTTGAAGCGCGACCACCCGGCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATagc                                                                                                                                                                                                                                                <  1:236493/300‑1 (MQ=255)
                                                ttggaagctcgACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATagc                                                                                                                                                                                                                                                <  1:236518/288‑1 (MQ=255)
                                                 tttAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAACGGGCGTTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTTAGCAGATTGCGGCAGGCCAATGGTGTTTTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATagc                                                                                                                                                                                                                                                <  2:397045/297‑1 (MQ=255)
                                                          gACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTc                                                                                                                                                                                                                                    <  1:625981/300‑1 (MQ=255)
                                                               gCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCTCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCAc                                                                                                                                                                                                                               <  2:392874/300‑1 (MQ=255)
                                                               gCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCAc                                                                                                                                                                                                                               <  1:601179/300‑1 (MQ=255)
                                                                  ttCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTAc                                                                                                                                                                                                                            <  2:147564/300‑1 (MQ=255)
                                                                     aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc                                                                                                                                                                                                                          <  1:51929/295‑1 (MQ=255)
                                                                     aatgcccAGCGACTGACACAACATACTGTAGCGGACGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACACATTGGTGAAAATCACAATGCCGTCAGCAAATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCACCTCCAACGCCTGACTCGCCAGCATTTTCTGCCGCGCCCGAAGGTCGCGACTGCTCCATagc                                                                                                                                                                                                                                                <  2:43758/273‑1 (MQ=255)
                                                                               gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACAcc                                                                                                                                                                                                                <  2:666038/299‑1 (MQ=255)
                                                                                            atacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCa                                                                                                                                                                                                   >  2:399275/4‑300 (MQ=255)
                                                                                            atacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCa                                                                                                                                                                                                   >  1:390416/4‑300 (MQ=255)
                                                                                                           gTGATATCGCGCATCGNGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGGTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCANAATACACCTGCCGCCAGCNCAATACGCAACGAGGTc                                                                                                                                                                                    >  1:589072/1‑300 (MQ=255)
                                                                                                           gTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCACGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTc                                                                                                                                                                                    >  1:248187/1‑300 (MQ=255)
                                                                                                             gATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCca                                                                                                                                                                                  >  1:262819/1‑300 (MQ=255)
                                                                                                                 tCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAgc                                                                                                                                                                              <  2:633517/300‑1 (MQ=255)
                                                                                                                     gcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCt                                                                                                                                                                          >  1:401063/1‑300 (MQ=255)
                                                                                                                               tttcagtaCTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCAc                                                                                                                                                                <  1:85608/293‑1 (MQ=255)
                                                                                                                               ttCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCAc                                                                                                                                                                <  1:85977/300‑1 (MQ=255)
                                                                                                                                        tttCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGc                                                                                                                                                       >  1:650394/1‑300 (MQ=255)
                                                                                                                                             gcgcATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGATTGTGCCAGCa                                                                                                                                                  >  1:277194/1‑300 (MQ=255)
                                                                                                                                                                                                  gATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGANGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGgc                                                                                           >  2:546316/1‑300 (MQ=255)
                                                                                                                                                                                                      gcggCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGgcagca                                                                                       <  2:580623/298‑1 (MQ=255)
                                                                                                                                                                                                                aatggtgtTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGAAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAAcagcagc                                                                               <  1:399275/294‑1 (MQ=255)
                                                                                                                                                                                                                  tggtgtTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAgcg                                                                           <  1:168227/298‑1 (MQ=255)
                                                                                                                                                                                                                    gtgtTGTTGAGCAGTANAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCAGCGACGATNCACCGCCAGCAACGGTGGNCAGCAAcagcagggagc                                                                           >  1:574829/1‑295 (MQ=255)
                                                                                                                                                                                                                                  tAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGNAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAacac                                                           <  2:502307/300‑1 (MQ=255)
                                                                                                                                                                                                                                  tAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAacac                                                           <  1:288685/300‑1 (MQ=255)
                                                                                                                                                                                                                                              ggTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCaa                                                 <  1:199281/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ggCGAGGATCGACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCAAAGCACCNGTTGTGCCAGCACTTCCGGCAAATGGCGATTTAGCCCATCGACGATCAACGCAAGCAACGGTGGCAGCAACAGCAGCAGCGCCGGCACAATTAACACCGTGTAGCAGATGCGGATATGCAGACGATCGTcc                         >  2:213350/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                      ggCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAAGTCCACAGCGCCTGCCACAGCACCAGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCAACGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCATAATTAACACCGTGTAGCAAATGCGGATATGCAGACGATCGTcc                         >  2:558104/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                            gATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCAGCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGTAGACGATCGTCCGGGTcgc                  >  2:447367/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                            aTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGcagcag  <  1:49952/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              ctcCAGCGCCTGACCCGCCAGTATTTTCTGCCGCGCCAGAAGTTCGCGACTACTCCATAGCCGCATCATGGTAAGCACTACGCACAATACACCTGCCGACAGCGCAATACGCAACGAGGTCCACAGCGCCTGCAACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCGATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCAAGcaca                                                                   >  2:72583/1‑232 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                   
GGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAG  >  NC_000913/73121‑73702

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: