Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 674,992 C→A intergenic (‑209/‑26) leuS ← / → ybeL leucine‑‑tRNA ligase/DUF1451 domain‑containing protein YbeL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913674,9920CA98.3% 209.0 / NA 58intergenic (‑209/‑26)leuS/ybeLleucine‑‑tRNA ligase/DUF1451 domain‑containing protein YbeL
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (32/25);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ATCTCTTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACGATGAACAAGGTTGCTCAATATTACCGTGAACTGGTTGCGTCACTGAGCGAACGCCTGCGCAATGGCGAACGTGATATCGACGCACTGGTGGAACAGGCGCGCGAGCGCGTAATAAAAACAGGGGAGTTAACGCGAACCGAGGTCGATGAGCTGACGCGAGCTGTCAGACGTGACCTGGAAGAGTTCGCCATGAGCTATGAAGAGAGCCTGAAAGAAGAA  >  NC_000913/674755‑675236
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
aTCTCTTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGAccc                                                                                                                                                                                                                                          <  2:347016/250‑1 (MQ=255)
    ctTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTTTATTATTATAGTAAGTTAACGACCCg                                                                                                                                                                                                                                         <  1:60217/247‑1 (MQ=255)
        cGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATg                                                                                                                                                                                                                                  >  1:171472/1‑250 (MQ=255)
            cGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACg                                                                                                                                                                                                                             >  1:26779/1‑251 (MQ=255)
             ggTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGg                                                                                                                                                                                                                            >  1:394608‑M1/1‑250 (MQ=255)
                aTTGCTCTTGCATGGCAGCCAGGGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCATTTCGCGCACTCGTCGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCAGGGAGATGAAACGgcat                                                                                                                                                                                                                         >  1:342505‑M1/1‑247 (MQ=255)
                  tGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatcc                                                                                                                                                                                                                       >  2:4581‑M1/1‑245 (MQ=255)
                    ctctTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccgg                                                                                                                                                                                                                     >  1:204950‑M1/1‑243 (MQ=255)
                      ctTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccgg                                                                                                                                                                                                                     <  1:356153‑M1/249‑9 (MQ=255)
                              gCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgcc                                                                                                                                                                                                           <  2:131226‑M1/251‑19 (MQ=255)
                               cAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccg                                                                                                                                                                                                          <  1:267361‑M1/251‑20 (MQ=255)
                                  ccAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGACAAAACAGCCTTTCGCGCACTCCACGTTGAAATGATGCCCGGATTACTCATACATTTATTTACAGAGTATGGGGGCGTATTAGCAAAGCAAGGAACAAAGAACGGCTATTATTATAGTAAGTTAACGAACCGGGAGATGAAACGgcatctggtgatggtgccgaa                                                                                                                                                                                                        >  1:67318‑M1/1‑229 (MQ=255)
                                   cAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaat                                                                                                                                                                                                      >  1:158537‑M1/1‑228 (MQ=255)
                                    aGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTCCTCATACATTAATTTACAGCGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaat                                                                                                                                                                                                      >  2:427344‑M1/1‑227 (MQ=255)
                                     gTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTTAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGGCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatcctgtgagcgtgccgaaatat                                                                                                                                                                                                    <  2:158537‑M1/251‑26 (MQ=255)
                                       ggTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgg                                                                                                                                                                                                  >  1:77498‑M1/1‑224 (MQ=255)
                                          ccTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggat                                                                                                                                                                                               <  1:489759‑M1/251‑31 (MQ=255)
                                             gTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgt                                                                                                                                                                                             >  1:381793‑M1/1‑218 (MQ=255)
                                                ttCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgttccgaaatatgggatgtatt                                                                                                                                                                                          <  2:161157‑M1/250‑36 (MQ=255)
                                                  cAATACGGCTACAAATGTAGCGTTGAGGTGNTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattcc                                                                                                                                                                                        <  2:292369‑M1/250‑38 (MQ=255)
                                                   aaTACGGCTACAAATGTCGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccgg                                                                                                                                                                                      <  2:383001‑M1/251‑40 (MQ=255)
                                                   aaTACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccgg                                                                                                                                                                                      >  2:465388‑M1/1‑212 (MQ=255)
                                                          cTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattcaggcacgat                                                                                                                                                                                >  1:556310‑M1/1‑205 (MQ=255)
                                                              aaaTGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgggccnaaatatgggatggattccggcacgataaga                                                                                                                                                                            >  2:284222‑M1/1‑201 (MQ=255)
                                                                                       aGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgaagc                                                                                                                                                    >  2:74402‑M1/1‑176 (MQ=255)
                                                                                         aTCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagac                                                                                                                                                <  1:274289‑M1/251‑78 (MQ=255)
                                                                                            cGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAACGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactc                                                                                                                                              >  1:505101‑M1/1‑171 (MQ=255)
                                                                                              catagcagagCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatc                                                                                                                                           >  1:389182‑M1/1‑169 (MQ=255)
                                                                                                 agcatagcCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatggattccggcacgataagaagggattatttacggcgctgacgcagactcatcaac                                                                                                                                        >  2:222552‑M1/1‑166 (MQ=255)
                                                                                                      agcCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagc                                                                                                                                   >  2:59156‑M1/1‑161 (MQ=255)
                                                                                                         ccAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagc                                                                                                                                <  2:247712‑M1/251‑94 (MQ=255)
                                                                                                            aaCGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaa                                                                                                                              >  1:291796‑M1/1‑155 (MQ=255)
                                                                                                                    cgTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaac                                                                                                                      >  1:3790‑M1/1‑147 (MQ=255)
                                                                                                                       cAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatg                                                                                                                  >  1:251723‑M1/1‑144 (MQ=255)
                                                                                                                                cTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagc                                                                                                          >  1:2755‑M1/1‑135 (MQ=255)
                                                                                                                                    cgcgCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgc                                                                                                                 >  1:87699‑M1/1‑131 (MQ=255)
                                                                                                                                    cgcgCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgc                                                                                                                 <  2:87699‑M1/239‑109 (MQ=255)
                                                                                                                                       gcACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccac                                                                                                    >  1:148495‑M1/1‑128 (MQ=255)
                                                                                                                                          cTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTGACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTAGAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagact                                                                                                                                               <  2:439107‑M1/203‑79 (MQ=255)
                                                                                                                                                gTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagtt                                                                                         >  1:390664‑M1/1‑119 (MQ=255)
                                                                                                                                                  tGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttg                                                                                        <  1:317718‑M1/250‑134 (MQ=255)
                                                                                                                                                     aaTGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccgg                                                                                    >  1:52758‑M1/1‑114 (MQ=255)
                                                                                                                                                               ggATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTACCGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcccccacagcggccagttcccggtacgtgcgt                                                                           >  2:44772‑M1/1‑104 (MQ=255)
                                                                                                                                                                    actcaTACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACACGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgatacgaagggattatttacgtcgctgacgcagactcatcaacacagcagcaacaccaaacaatgccgtccgcacccacagcggccagatgacggcacgtccgtatgggg                                                                     >  2:22091‑M1/3‑99 (MQ=255)
                                                                                                                                                                       cATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtga                                                                  <  1:32993‑M1/251‑156 (MQ=255)
                                                                                                                                                                          aCATTAATTTACAGAGTTTGTGGGCGTATTANCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtc                                                               <  2:251723‑M1/251‑159 (MQ=255)
                                                                                                                                                                          aCATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtaggtc                                                               >  2:24057‑M1/1‑93 (MQ=255)
                                                                                                                                                                          aCATTAATTGACAGAGTTTGTGTGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagt                                                                <  2:52758‑M1/250‑158 (MQ=255)
                                                                                                                                                                                  ttACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggaggtc                                                       >  1:344632‑M1/1‑85 (MQ=255)
                                                                                                                                                                                    aCAGAGTTTGTTGGCGTATTAGCAAAGCAAGGAGCAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagagattatttacgtcgctgacgctgactcatcaacacagcagcaacaccaaacaatgccgtcagctcccacagcgggcagttgccggtacgtgcgtatggttggagtccggtggtcg                                                      <  2:60056‑M1/250‑168 (MQ=255)
                                                                                                                                                                                      agagTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcg                                                   <  2:129353‑M1/251‑171 (MQ=255)
                                                                                                                                                                                             gtgGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgtt                                            <  1:4581‑M1/251‑178 (MQ=255)
                                                                                                                                                                                                     aTTAGCAAAGCAAGGAACAAAGAACGTCTATTATGATAGTAAGTTAACGACCAGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccggcagcacccacagcggccagttgccggtacgtgcgtagggtgtgagtccggtggtcggcgtcacgttagtggtta                                    <  2:171472‑M1/251‑186 (MQ=255)
                                                                                                                                                                                                                         aGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgacctgc                >  1:427625‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                                                  attcttatagtaAGGTAATGACCCGGGAGATGAAAAGgctgccggtgagcgtgccgatatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagacgcagcaacacagcagcaaaaccaaacaatgccgtcagcgcccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaactcctcgcgggtgaattgcgggatcaac       <  2:554898‑M1/247‑215 (MQ=255)
                                                                                                                                                                                                                                  attattatAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcgggatcat        <  1:509532‑M1/250‑214 (MQ=255)
                                                                                                                                                                                                                                  attattatAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgcgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcgggatcat        <  2:381793‑M1/250‑214 (MQ=255)
                                                                                                                                                                                                                                     attatAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcgggatcatcgc     <  1:311243‑M1/250‑217 (MQ=255)
                                                                                                                                                                                                                                       tataGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtagatcccggtggccgccgtacgataatggaaaacacctcgcgggcgaacggcgggagcagcgcctg  >  2:166474‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
ATCTCTTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACGATGAACAAGGTTGCTCAATATTACCGTGAACTGGTTGCGTCACTGAGCGAACGCCTGCGCAATGGCGAACGTGATATCGACGCACTGGTGGAACAGGCGCGCGAGCGCGTAATAAAAACAGGGGAGTTAACGCGAACCGAGGTCGATGAGCTGACGCGAGCTGTCAGACGTGACCTGGAAGAGTTCGCCATGAGCTATGAAGAGAGCCTGAAAGAAGAA  >  NC_000913/674755‑675236

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