Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,622 C→T V56V (GTG→GTA 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,6220CT97.0% 202.9 / ‑2.3 68V56V (GTG→GTAytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (32/33);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTT  >  NC_000913/4624376‑4624863
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
cGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGACCGCTTAAGCGGTGCAACGTTCCGCGTTACCt                                                                                                                                                                                                                                               >  2:418253/1‑251 (MQ=255)
      gTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGAcgcg                                                                                                                                                                                                                                         <  1:304369/251‑1 (MQ=255)
                             aaaTCTGTTGGGTAAAATAGCCGCCCGCTTTTTTACCGTCGGGCGCCAGCCGGTCGCGCACTTTGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTcc                                                                                                                                                                                                                  <  2:132392/250‑1 (MQ=255)
                                   gctgGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAATGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAg                                                                                                                                                                                                                <  2:47189/247‑1 (MQ=255)
                                   gctgGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAg                                                                                                                                                                                                                <  1:409850/247‑1 (MQ=255)
                                            aaTAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACgt                                                                                                                                                                                                   <  1:129227/251‑2 (MQ=255)
                                                 ccgcacgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATAGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTAGGTACGCATCAGCGGTGCCACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCAGGGGATTACGTTGTCg                                                                                                                                                                                              >  2:430795/7‑251 (MQ=255)
                                                   gcccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTTCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCt                                                                                                                                                                                            <  1:301412/247‑1 (MQ=255)
                                                         ttttttACCGTCGAGCGGCAGCCGGTCGCGCACTTCGCCGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGTATGCGATTTTCATCCGGCGAGTCGGGACGCCTCAGCGGTGAACCGTTCAGCGTTACCTGACGAGCCTAGGTGCGAGCAAGTTCTTCCAGCTGTGGTTTTCGTTGTCGCtggtggg                                                                                                                                                                                      >  2:263503/1‑250 (MQ=255)
                                                              tACCGTCGGGCGCCAGCCGGGCGCGCACTTCGTCGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGAGCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGttttt                                                                                                                                                                                  <  2:424130/250‑2 (MQ=255)
                                                                         gcCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGATAGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCCGTTCTTCCCGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATg                                                                                                                                                                       >  1:375644/1‑250 (MQ=255)
                                                                           cAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATAAGGCTTTCATACGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTAGATGGAGGTTTTGACTAAACCATGa                                                                                                                                                                      >  1:320636/1‑249 (MQ=255)
                                                                               cGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTcc                                                                                                                                                                <  1:439504/251‑1 (MQ=255)
                                                                               cGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGGTTTGACTAAACCATGACGCTcc                                                                                                                                                                >  2:228040/1‑251 (MQ=255)
                                                                                     cgcACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCAt                                                                                                                                                          <  1:240000/251‑1 (MQ=255)
                                                                                        aCTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAg                                                                                                                                                       <  1:276184/251‑1 (MQ=255)
                                                                                             gcagCTTTCGCCAGAACGGGCGATAAGATCGCCGTGTTCGTCAGTCACGCCCGCGTCGAGGATACGGTTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCa                                                                                                                                                  <  2:198195/248‑1 (MQ=255)
                                                                                                gCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCAAGGGTGCGGGCCAGCTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGCCGCTACCTGCATCAGCGCCACt                                                                                                                                                >  1:454325/1‑250 (MQ=255)
                                                                                                 cTTTCGCCAGAACGGGCGCTAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTACACCATGACGCTCCCTGCATCAGCGCCACTac                                                                                                                                              >  2:314930/1‑250 (MQ=255)
                                                                                                        ccgaacaggCGATAAGATCGCCTTTTTCCTCATCCATGCCCGCGTCGAGGATACGGCTTTCATCCGTTATCTGATCGAGGATCGCTTGTATGTGTTGTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTGTCCTGTCGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGGCGCTGGTGGTTTTGACGAAACCATGACGCTCCCTTCATCAGCGCCACTAACAAGGCa                                                                                                                                       <  2:171556/243‑1 (MQ=255)
                                                                                                           aaCGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGa                                                                                                                                    <  1:241672/251‑1 (MQ=255)
                                                                                                               ggCGATAAGATCGCCTTGGTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGa                                                                                                                                <  1:227361/250‑1 (MQ=255)
                                                                                                                  gATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGGTACCTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAAACATGTCGCTCCCTGCATCAGCGCACCTAACAAGGCTAGaacgaac                                                                                                                              >  2:384575/1‑244 (MQ=255)
                                                                                                                            ccTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACCATccc                                                                                                                   >  2:266524/1‑249 (MQ=255)
                                                                                                                               tGTTCGTCATACACGCCCGCGGCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTTCGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACt                                                                                                                  <  2:341866/249‑1 (MQ=255)
                                                                                                                                   cGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGCGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCg                                                                                                             >  1:431250/1‑250 (MQ=255)
                                                                                                                                    gTCATACACGCACGCGTCGAGGATAAGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGGGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCAAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCCTGACGCACCCTGCATCAACGCCACTAACAAGGCAAGACAGAACAGGACAATCACAGCCCGAt                                                                                                           >  1:262884/1‑251 (MQ=255)
                                                                                                                                      cATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGc                                                                                                         >  1:290144/1‑251 (MQ=255)
                                                                                                                                          cacgcccGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGcc                                                                                                     >  1:287289/5‑251 (MQ=255)
                                                                                                                                                       gggATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGttt                                                                                        <  2:237649/250‑1 (MQ=255)
                                                                                                                                                       aGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGttt                                                                                        >  1:352309/1‑251 (MQ=255)
                                                                                                                                                        ggATACGGCGTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTGTCACCCGGCGAGTCGGTACGCAGCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGtttt                                                                                       <  2:431250/251‑1 (MQ=255)
                                                                                                                                                          aTACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTtgt                                                                                     >  2:249962/1‑251 (MQ=255)
                                                                                                                                                              ggCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGa                                                                                 <  2:352309/251‑1 (MQ=255)
                                                                                                                                                                cTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATAACTGCCCGATGCAGCCGACATTTCAGTTTTGTGCGAGc                                                                               >  1:221253/1‑251 (MQ=255)
                                                                                                                                                                 tttCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGcc                                                                              >  1:247783/1‑251 (MQ=255)
                                                                                                                                                                         gTTAACTGATCGAGGATCGCCTTAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCa                                                                      <  1:388097/251‑1 (MQ=255)
                                                                                                                                                                            aaCTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGCGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGAc                                                                                                                                                                                >  2:28187/1‑142 (MQ=255)
                                                                                                                                                                            aaCTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGAc                                                                                                                                                                                <  1:28187/142‑1 (MQ=255)
                                                                                                                                                                                gATCGAGGAGCCTCTGAATGCGTTTTTCATCCGGCGAGCCGGTACGCATCAGCGGTGCAACGTTCAGCGTTCCCTGACGCGCGAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGGTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGTCGGAAGTTTAGGTTTGTGCGAGCCATCTTCCTCCCTTTg                                                               <  1:405491/251‑1 (MQ=255)
                                                                                                                                                                                    gagtgtCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCCGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGaaaa                                                           <  2:433673/246‑1 (MQ=255)
                                                                                                                                                                                    gAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGaaaa                                                           <  1:191121/251‑1 (MQ=255)
                                                                                                                                                                                     aGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCCCCCTTTGAAAAt                                                          >  2:169848/1‑251 (MQ=255)
                                                                                                                                                                                               gAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACtt                                                 >  2:197682/1‑250 (MQ=255)
                                                                                                                                                                                               gAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGAAAATCAATGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCGTCTTCCACCCTTTGAAAATTTGAGCCtt                                                 >  2:82185/1‑250 (MQ=255)
                                                                                                                                                                                                 aTGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAAt                                              <  1:249962/251‑1 (MQ=255)
                                                                                                                                                                                                        tttCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGcc                                        <  1:159467/250‑1 (MQ=255)
                                                                                                                                                                                                         ttCATCCGGCGAGTCGGGACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGTCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCTCCCGTTGAAAATTTGAGACTTAATGTTGCCa                                       <  1:106489/250‑1 (MQ=255)
                                                                                                                                                                                                          tCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGcc                                        >  1:7914/1‑248 (MQ=255)
                                                                                                                                                                                                               cggcgattCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTgag                                                     <  2:329487/223‑1 (MQ=255)
                                                                                                                                                                                                               cggcgAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTgag                                                     >  1:329487/5‑230 (MQ=255)
                                                                                                                                                                                                                ggcgAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCaa                                                                                                                                      >  1:407624/4‑148 (MQ=255)
                                                                                                                                                                                                                ggcgAGTCGGGACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCaa                                                                                                                                      <  2:407624/145‑1 (MQ=255)
                                                                                                                                                                                                                    aGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAg                                                                                <  2:237378/198‑1 (MQ=255)
                                                                                                                                                                                                                    aGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAg                                                                                >  1:237378/1‑198 (MQ=255)
                                                                                                                                                                                                                      gcGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCGGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACGAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATAc                         <  2:437045/250‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTCACAAGGCAAGACAGAACAGGACCATCCCTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGATAATTTGAGACTTAATGTTGCCAGAAGCAATGGAtcc                         >  2:125658/1‑248 (MQ=255)
                                                                                                                                                                                                                       cGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCAGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACACGGCAAGACAGAACAGGACCATCACTGCCCGATGCAGCCGGAATTTCAGTTTCGTGCGAGCCATCTTCCACCCTTTGAAACTTTGAGACTTAATGTTGCCAGAAGCAAAGGATAcc                        >  2:71185/1‑250 (MQ=255)
                                                                                                                                                                                                                           aCGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAgg                     <  1:108130/250‑1 (MQ=255)
                                                                                                                                                                                                                             gccacAGCGGTGCAACGGTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAg                  <  2:117005/247‑1 (MQ=255)
                                                                                                                                                                                                                             gCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGACGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAGTTTGAGACTTAATATTGACAGAAGCAATGGATACAAgga                    >  2:263466/1‑248 (MQ=255)
                                                                                                                                                                                                                             gCATCAGCGGTGCAACGTTCAGCGTTACCCGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGGTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCGTGGAAAATTTGAGACTTAATGTTGCCAGAAGCAAGGGATACAAGGTAg                  <  2:454325/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aGCGGTGCAACGTTCAGCGTTCCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCa             <  1:116975/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATg           <  1:247142/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGCTTTGACTAAACCATGCCGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGATACTTAATGTTGCCAGAAGCAATGGATCAAAGGTAGCCTCATGc          >  2:97224/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGACGCAATGGATACAAGGTAGCCTCATGCGtt       >  1:250214/1‑250 (MQ=255)
                                                                                                                                                                                                                                             gTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTAtttt  <  1:125658/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gTTCACCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTAtttt  >  1:130108/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
CGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTT  >  NC_000913/4624376‑4624863

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