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

ATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTGA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCC  >  NC_000913/4624343‑4624825
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
aTGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCg                                                                                                                                                                                                                                            <  1:55597/250‑1 (MQ=255)
aTGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGAGCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCg                                                                                                                                                                                                                                            <  2:240367/250‑1 (MQ=255)
           gtgtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATGGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGc                                                                                                                                                                                                                                <  2:90514/251‑1 (MQ=255)
               cAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GGCGGTACGCATCAGCgggg                                                                                                                                                                                                                            >  2:3470/1‑249 (MQ=255)
                atgcgaaGATAGCCGAGCGGTCCTTTTTGCCCCGCAATTGGCTCGACAATCTGGTGGTTAAAAGAGCCGCCTGTTGTTTTACCGTCGGGCGCCAGCCGGTCGCGCGCTTCGACTCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCTAGGATACGGCTTTCATCCGGTAACTGAGCGAGGATCGCCGGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACTCATCAGCGGTGc                                                                                                                                                                                                                           <  2:199715/245‑1 (MQ=255)
                          aGCCGCGCGGTTCGTTTTTCCCCGCAATTGGCCCGACAATCGGCTGGTAAAAATAGCCGCCTGCTTTTTTCCCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATTCACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGc                                                                                                                                                                                                                 <  2:212766/251‑1 (MQ=255)
                                       gTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGAcgcg                                                                                                                                                                                                    >  1:255882/1‑251 (MQ=255)
                                           ttACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCa                                                                                                                                                                                                 <  1:291754/250‑1 (MQ=255)
                                              cccGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGt                                                                                                                                                                                             <  1:290414/251‑1 (MQ=255)
                                                   aaTTGGCTCGCCAATCTGCTGGTTATAATCGCCGCCTGCTTTTTTCCCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGTGCGATAAGATCGCCTTGTTCGTCATACACCCCCGCGTCGCGGCTACGGCTTTCATCCTTTACCTGATCGAGCATCGCCTTAATGC‑GTTTTTCATCCGGTAA‑GTCGGTCCGCATCAGCGGTTCAACTTTCAGCGTCACCttcctcgccacgttgcggg                                                                                                                                                                                        >  2:54577/1‑233 (MQ=255)
                                                    aTTGGCTCGACAATATGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGGCGCGCACTTCGACGCTTTCGCCAAAACGTGCGATAAGATCGCCTTGTTCGTCATACACCCCCGCGTAGGGGGTACGGCTTTCATCCGTTAACTGATCGAGGCACGCCTGAATGC‑GTTTTTTATTCCGTAA‑GTCGGTACGCCTCCGCGGGGCAACGTTCAGCGTCCCCTGACGCGCCAGGGCGCGGGc                                                                                                                                                                                       >  2:313603/1‑251 (MQ=255)
                                                       ggCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGcc                                                                                                                                                                                      >  1:116622/1‑249 (MQ=255)
                                                             aCAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGttcttc                                                                                                                                                                              >  1:274673/1‑251 (MQ=255)
                                                              cAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTcc                                                                                                                                                                             >  1:60068/1‑251 (MQ=255)
                                                                    gctgGTTAAAATAGCCTCCTGCTTTTTTACCATCGCGCGCCAGCCGTTCGCGCTCTTCGACGCTTTCGCCAGAACGTGCGATAAGTTCGCCTTTTTCTTCATACACACCCTCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGTGTGC‑GTTTTTCATCCG‑TGAAGTCGGTACGCATCTGCGGTGCTACGTTCAGCGTCCCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAg                                                                                                                                                                           <  1:20118/247‑1 (MQ=255)
                                                                     ctgGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCtgtg                                                                                                                                                                      >  1:54981/1‑251 (MQ=255)
                                                                                   cgccTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCCGCTGTGGATTACGCTGTCGc                                                                                                                                                        >  2:71531/1‑251 (MQ=255)
                                                                                    gccTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCCCCGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCt                                                                                                                                                       >  2:140889/1‑251 (MQ=255)
                                                                                    gccTGCTTTTTTAACGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGACAGAACGTGCGATAAGATCGCCTTGTTCGTAATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGATTGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAAAGTTCAGCGTCACCTGCCGCGCCAGGGTGCAGGCCACGTCTTCCAGCTGAGGATTACGCTTTCGCt                                                                                                                                                       >  2:124282/1‑251 (MQ=255)
                                                                                     ccTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTTCTGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCt                                                                                                                                                       <  1:324253/250‑1 (MQ=255)
                                                                                               tacctgcGAGCGCCAGCCGGTCGCGCACTTCGACGCTTGCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATGCACACCCGCGTCGAGGATACGGCTTTCATACGTTAACTGATCGAGTATCGCCTGAATGC‑GTTTTTCAGCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTg                                                                                                                                             <  2:116788/244‑1 (MQ=255)
                                                                                                                cGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACTCACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTGTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTTCATCGTTCAGCGGCACCTGACGCGCCTGGGTGCGGGCCAGTTCGTCCAGCTGTGTATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTcc                                                                                                                           <  1:344242/251‑1 (MQ=255)
                                                                                                                        cACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGACTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATc                                                                                                                    >  2:18519/1‑250 (MQ=255)
                                                                                                                             cGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGcc                                                                                                              >  2:7809/1‑251 (MQ=255)
                                                                                                                             cGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCGTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCGGTCTCTTGTGGTTTTGACTAAACCATGGCGCTCCCTGCATCAGCGcc                                                                                                              <  1:200834/251‑1 (MQ=255)
                                                                                                                                    ttCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCGCCTGACGCGCCAGGGTGCGGGCCAGTTCTGCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTaa                                                                                                         <  1:314471/249‑1 (MQ=255)
                                                                                                                                      cGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGTATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACaa                                                                                                      <  2:54981/250‑1 (MQ=255)
                                                                                                                                        ccAGAACGTGAGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGc                                                                                                   <  2:252440/251‑1 (MQ=255)
                                                                                                                                                           cGCCTAGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GACGGTACGAATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCCCTAACAAGGCAAGACAGAACAGGACAAt                                                                                 >  1:259715/1‑250 (MQ=255)
                                                                                                                                                                 gTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGcc                                                                         <  1:177081/251‑1 (MQ=255)
                                                                                                                                                                  ttCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGCTCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGCCCCCTAACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTCCGCTGTCGCTGGTGGTTTTGACTAAACCCTGACGCTCCCTGCATCAGCGCCACTAACACGGCAAGACAGACCAGGACAATCCC‑TGccc                                                                        >  2:256309/1‑251 (MQ=255)
                                                                                                                                                                        aTACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGc                                                                   >  1:324901/1‑250 (MQ=255)
                                                                                                                                                                           cacacCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGcc                                                               <  1:186951/251‑1 (MQ=255)
                                                                                                                                                                            acacCCGCGTCGAGGATACGGATTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCg                                                              <  1:240905/251‑1 (MQ=255)
                                                                                                                                                                                 cgcgTCGAGGATACGGCTTTCATCCGTTAACTTATCGAGGATCGCCTGATTGC‑GTTTTTCATCCGGTAT‑GTCGGTACGCATCAGCGGTTCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTTTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAAtt                                                         <  1:269665/251‑1 (MQ=255)
                                                                                                                                                                                       gAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCGTCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCAGGAATTTCAGtt                                                   >  2:317709/1‑251 (MQ=255)
                                                                                                                                                                                         ggagaCGGCTTTCATCCGTTAACTGAGCGAGGGTCGCCTGAAGGC‑GTTTTTCATCCG‑TGAAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGGCATCTGGCGCGCCAGGGTGCGTGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCAGCAGCGCCACTAAGAATGCTAGACAGAACAGGACATTCAC‑TGCCCGATGCAGCCGGAATTTCAGtttt                                                 <  2:215641/247‑1 (MQ=255)
                                                                                                                                                                                              cGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGAACGCCTCAGCGGTGCCACGTTCAGCGTCCCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCCTGACGCTCCCTGCATCAGCGCCACTCACAAGGCAAGACCGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTTTTGTGCg                                            >  2:144645/1‑251 (MQ=255)
                                                                                                                                                                                               ggCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTTTTGTGCg                                            >  2:281919/1‑250 (MQ=255)
                                                                                                                                                                                                       tCCGTTAACTGATCGCGGAACGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCtt                                   >  2:87872/1‑251 (MQ=255)
                                                                                                                                                                                                       tCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGGTTTGTGCGAGCCATCtt                                   >  1:129382/1‑251 (MQ=255)
                                                                                                                                                                                                        ccgatAACTGATCGAGGATCGCTAGAATGC‑GGTTTTCATCCGGTAA‑GTCGGTACGCAGCCGCGGTGCAACGTTCAGCGTCATCTGACGCGCCAGGGTGCGGGCCAGGTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCAGCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTc                                  <  2:274673/247‑1 (MQ=255)
                                                                                                                                                                                                                gtagcGAGGATCGCCTGAATGCTGTTTTT‑ATCCGGTAA‑GTCGGTACGCATCAGCGTTGCAACGTTCCGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCGCCCtt                           <  2:297480/246‑1 (MQ=255)
                                                                                                                                                                                                                    cGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCTCCCTTTGaaa                      <  1:260042/251‑1 (MQ=255)
                                                                                                                                                                                                                      aGGATCGCCTGAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAAt                    >  1:325160/1‑251 (MQ=255)
                                                                                                                                                                                                                                gAATGC‑GTTTTTCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTGTTGTGCGAGCCATCTTCCTCCCTTTGAAAATTTGAGACtt           <  2:67354/250‑1 (MQ=255)
                                                                                                                                                                                                                                        tttttCTTCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGaa                                                                                          >  1:96565/1‑165 (MQ=255)
                                                                                                                                                                                                                                        tttttCATCCGGTAA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGaa                                                                                          <  2:96565/165‑1 (MQ=255)
                                                                                                                                                                                                                                           ttCATCCGGTAC‑GTCGGTACGCATCCGCGGAGCATCGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGCGGATAACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCACGACAGAAAAGGACACT‑ACATGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCACCCCTTTGTAAATTTGAGACTTAATGttacc  >  2:233126/1‑246 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
ATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGC‑GTTTTTCATCCGGTGA‑GTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCAC‑TGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCC  >  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: ‑