Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,587 G→A S68L (TCA→TTA)  ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,5870GA95.9% 145.8 / ‑6.0 49S68L (TCA→TTA) ytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (24/23);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

ATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCC  >  NC_000913/4624343‑4624825
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
aTGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCg                                                                                                                                                                                                                                            <  1:55597/250‑1 (MQ=255)
aTGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGAGCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCg                                                                                                                                                                                                                                            <  2:240367/250‑1 (MQ=255)
           gtgtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATGGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGc                                                                                                                                                                                                                                <  2:90514/251‑1 (MQ=255)
               cAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGGCGGTACGCATCAGCgggg                                                                                                                                                                                                                            >  2:3470/1‑249 (MQ=255)
                atgcgaaGATAGCCGAGCGGTCCTTTTTGCCCCGCAATTGGCTCGACAATCTGGTGGTTAAAAGAGCCGCCTGTTGTTTTACCGTCGGGCGCCAGCCGGTCGCGCGCTTCGACTCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCTAGGATACGGCTTTCATCCGGTAACTGAGCGAGGATCGCCGGAATGCGTTTTTCATCCGGTAAGTCGGTACTCATCAGCGGTGc                                                                                                                                                                                                                           <  2:199715/245‑1 (MQ=255)
                          aGCCGCGCGGTTCGTTTTTCCCCGCAATTGGCCCGACAATCGGCTGGTAAAAATAGCCGCCTGCTTTTTTCCCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATTCACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGc                                                                                                                                                                                                                 <  2:212766/251‑1 (MQ=255)
                                       gTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGAcgcg                                                                                                                                                                                                    >  1:255882/1‑251 (MQ=255)
                                           ttACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCa                                                                                                                                                                                                 <  1:291754/250‑1 (MQ=255)
                                              cccGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGt                                                                                                                                                                                             <  1:290414/251‑1 (MQ=255)
                                                   aaTTGGCTCGCCAATCTGCTGGTTATAATCGCCGCCTGCTTTTTTCCCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGTGCGATAAGATCGCCTTGTTCGTCATACACCCCCGCGTCGCGGCTACGGCTTTCATCCTTTACCTGATCGAGCATCGCCTTAATGCGTTTTTCATCCGGTAAGTCGGTCCGCATCAGCGGTTCAACTTTCAGCGTCACCttcctcgccacgttgcggg                                                                                                                                                                                        >  2:54577/1‑233 (MQ=255)
                                                    aTTGGCTCGACAATATGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGGCGCGCACTTCGACGCTTTCGCCAAAACGTGCGATAAGATCGCCTTGTTCGTCATACACCCCCGCGTAGGGGGTACGGCTTTCATCCGTTAACTGATCGAGGCACGCCTGAATGCGTTTTTTATTCCGTAAGTCGGTACGCCTCCGCGGGGCAACGTTCAGCGTCCCCTGACGCGCCAGGGCGCGGGc                                                                                                                                                                                       >  2:313603/1‑251 (MQ=255)
                                                       ggCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGcc                                                                                                                                                                                      >  1:116622/1‑249 (MQ=255)
                                                             aCAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGttcttc                                                                                                                                                                              >  1:274673/1‑251 (MQ=255)
                                                              cAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTcc                                                                                                                                                                             >  1:60068/1‑251 (MQ=255)
                                                                    gctgGTTAAAATAGCCTCCTGCTTTTTTACCATCGCGCGCCAGCCGTTCGCGCTCTTCGACGCTTTCGCCAGAACGTGCGATAAGTTCGCCTTTTTCTTCATACACACCCTCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGTGTGCGTTTTTCATCCGTGAAGTCGGTACGCATCTGCGGTGCTACGTTCAGCGTCCCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAg                                                                                                                                                                           <  1:20118/247‑1 (MQ=255)
                                                                     ctgGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCtgtg                                                                                                                                                                      >  1:54981/1‑251 (MQ=255)
                                                                                   cgccTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCCGCTGTGGATTACGCTGTCGc                                                                                                                                                        >  2:71531/1‑251 (MQ=255)
                                                                                    gccTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCCCCGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCt                                                                                                                                                       >  2:140889/1‑251 (MQ=255)
                                                                                    gccTGCTTTTTTAACGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGACAGAACGTGCGATAAGATCGCCTTGTTCGTAATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGATTGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAAAGTTCAGCGTCACCTGCCGCGCCAGGGTGCAGGCCACGTCTTCCAGCTGAGGATTACGCTTTCGCt                                                                                                                                                       >  2:124282/1‑251 (MQ=255)
                                                                                     ccTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTTCTGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCt                                                                                                                                                       <  1:324253/250‑1 (MQ=255)
                                                                                               tacctgcGAGCGCCAGCCGGTCGCGCACTTCGACGCTTGCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATGCACACCCGCGTCGAGGATACGGCTTTCATACGTTAACTGATCGAGTATCGCCTGAATGCGTTTTTCAGCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTg                                                                                                                                             <  2:116788/244‑1 (MQ=255)
                                                                                                                cGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACTCACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTGTTCATCCGGTAAGTCGGTACGCATCAGCGGTTCATCGTTCAGCGGCACCTGACGCGCCTGGGTGCGGGCCAGTTCGTCCAGCTGTGTATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTcc                                                                                                                           <  1:344242/251‑1 (MQ=255)
                                                                                                                        cACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGACTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATc                                                                                                                    >  2:18519/1‑250 (MQ=255)
                                                                                                                             cGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGcc                                                                                                              >  2:7809/1‑251 (MQ=255)
                                                                                                                             cGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCGTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCGGTCTCTTGTGGTTTTGACTAAACCATGGCGCTCCCTGCATCAGCGcc                                                                                                              <  1:200834/251‑1 (MQ=255)
                                                                                                                                    ttCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCGCCTGACGCGCCAGGGTGCGGGCCAGTTCTGCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTaa                                                                                                         <  1:314471/249‑1 (MQ=255)
                                                                                                                                      cGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGTATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACaa                                                                                                      <  2:54981/250‑1 (MQ=255)
                                                                                                                                        ccAGAACGTGAGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGc                                                                                                   <  2:252440/251‑1 (MQ=255)
                                                                                                                                                           cGCCTAGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGACGGTACGAATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCCCTAACAAGGCAAGACAGAACAGGACAAt                                                                                 >  1:259715/1‑250 (MQ=255)
                                                                                                                                                                 gTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGcc                                                                         <  1:177081/251‑1 (MQ=255)
                                                                                                                                                                  ttCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGCTCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGCCCCCTAACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTCCGCTGTCGCTGGTGGTTTTGACTAAACCCTGACGCTCCCTGCATCAGCGCCACTAACACGGCAAGACAGACCAGGACAATCCCTGccc                                                                        >  2:256309/1‑251 (MQ=255)
                                                                                                                                                                        aTACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGc                                                                   >  1:324901/1‑250 (MQ=255)
                                                                                                                                                                           cacacCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGcc                                                               <  1:186951/251‑1 (MQ=255)
                                                                                                                                                                            acacCCGCGTCGAGGATACGGATTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCg                                                              <  1:240905/251‑1 (MQ=255)
                                                                                                                                                                                 cgcgTCGAGGATACGGCTTTCATCCGTTAACTTATCGAGGATCGCCTGATTGCGTTTTTCATCCGGTATGTCGGTACGCATCAGCGGTTCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTTTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAAtt                                                         <  1:269665/251‑1 (MQ=255)
                                                                                                                                                                                       gAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCGTCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCAGGAATTTCAGtt                                                   >  2:317709/1‑251 (MQ=255)
                                                                                                                                                                                         ggagaCGGCTTTCATCCGTTAACTGAGCGAGGGTCGCCTGAAGGCGTTTTTCATCCGTGAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGGCATCTGGCGCGCCAGGGTGCGTGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCAGCAGCGCCACTAAGAATGCTAGACAGAACAGGACATTCACTGCCCGATGCAGCCGGAATTTCAGtttt                                                 <  2:215641/247‑1 (MQ=255)
                                                                                                                                                                                              cGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGAACGCCTCAGCGGTGCCACGTTCAGCGTCCCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCCTGACGCTCCCTGCATCAGCGCCACTCACAAGGCAAGACCGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCg                                            >  2:144645/1‑251 (MQ=255)
                                                                                                                                                                                               ggCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCg                                            >  2:281919/1‑250 (MQ=255)
                                                                                                                                                                                                       tCCGTTAACTGATCGCGGAACGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCtt                                   >  2:87872/1‑251 (MQ=255)
                                                                                                                                                                                                       tCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGGTTTGTGCGAGCCATCtt                                   >  1:129382/1‑251 (MQ=255)
                                                                                                                                                                                                        ccgatAACTGATCGAGGATCGCTAGAATGCGGTTTTCATCCGGTAAGTCGGTACGCAGCCGCGGTGCAACGTTCAGCGTCATCTGACGCGCCAGGGTGCGGGCCAGGTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCAGCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTc                                  <  2:274673/247‑1 (MQ=255)
                                                                                                                                                                                                                gtagcGAGGATCGCCTGAATGCTGTTTTTATCCGGTAAGTCGGTACGCATCAGCGTTGCAACGTTCCGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCGCCCtt                           <  2:297480/246‑1 (MQ=255)
                                                                                                                                                                                                                    cGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCTCCCTTTGaaa                      <  1:260042/251‑1 (MQ=255)
                                                                                                                                                                                                                      aGGATCGCCTGAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAAt                    >  1:325160/1‑251 (MQ=255)
                                                                                                                                                                                                                                gAATGCGTTTTTCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTGTTGTGCGAGCCATCTTCCTCCCTTTGAAAATTTGAGACtt           <  2:67354/250‑1 (MQ=255)
                                                                                                                                                                                                                                        tttttCTTCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGaa                                                                                          >  1:96565/1‑165 (MQ=255)
                                                                                                                                                                                                                                        tttttCATCCGGTAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGaa                                                                                          <  2:96565/165‑1 (MQ=255)
                                                                                                                                                                                                                                           ttCATCCGGTACGTCGGTACGCATCCGCGGAGCATCGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGCGGATAACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCACGACAGAAAAGGACACTACATGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCACCCCTTTGTAAATTTGAGACTTAATGttacc  >  2:233126/1‑246 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
ATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCC  >  NC_000913/4624343‑4624825

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