Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,579,610 C→T E77E (GAG→GAA symE ← nucleoid‑associated protein SymE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,579,6100CT97.4% 113.7 / NA 39E77E (GAG→GAAsymEnucleoid‑associated protein SymE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/21);  minor base A (1/0);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCT  >  NC_000913/4579363‑4579847
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
ccTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCTCCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACGGCATCAGTTCGCTCTCTTcg                                                                                                                                                                                                                                            <  1:305914/251‑1 (MQ=255)
   ctATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTcggcg                                                                                                                                                                                                                                         >  2:384216/1‑251 (MQ=255)
                  cAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCg                                                                                                                                                                                                                          <  2:120782/251‑1 (MQ=255)
                             aaCAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAAt                                                                                                                                                                                                               <  1:121300/251‑1 (MQ=255)
                               cAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACAGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATAc                                                                                                                                                                                                             >  1:197031/1‑251 (MQ=255)
                               cAATAGATTTTTAATGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAACCGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATAc                                                                                                                                                                                                             >  1:364979/1‑251 (MQ=255)
                                     attattagtGCAAAATCAGACTGTTAAATATTTACTGAGGGGTAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACTCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCTGCTCCTGCTTTTGTTTACGCGCCGACAGTTTGGACACCTGGCGTAGCGACTTCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCtt                                                                                                                                                                                                       <  1:223553/243‑1 (MQ=255)
                                            cTGCAAAATCAGACTGGTGAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGCGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCCATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGTGGACAATACCGCCTTCCCTCAc                                                                                                                                                                                                >  2:376245/1‑251 (MQ=255)
                                                          tGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTAccg                                                                                                                                                                                  >  1:332847/1‑251 (MQ=255)
                                                           ggTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTAccgc                                                                                                                                                                                 <  1:402401/251‑1 (MQ=255)
                                                                              ggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCACTCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGGCAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGGGAGGACAATACAGCCTTCCATCACTTTGACCTCTACCGCCGTGCCCGTGGCCAAACcg                                                                                                                                                              >  2:375988/1‑251 (MQ=255)
                                                                                        tATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTcc                                                                                                                                                     >  1:222761/1‑250 (MQ=255)
                                                                                         aTTGAGTCAGTGGAAGGCTCTCGGTGGTTAACCGTGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGTCGGCGGGTGGTTGGGCGGGGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCat                                                                                                                                                   >  2:151376/1‑250 (MQ=255)
                                                                                         aTTGAGTCAGTGGAAGGCCCCCGGGGGTTACCCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGAAATCACTCCAATAAACGCCTCCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAg                                                                                                                                                   <  1:123401/251‑1 (MQ=255)
                                                                                                 aGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCCGTGGCAAAACCGGCGGCTTCCAGCCACTGAc                                                                                                                                           >  1:110122/1‑251 (MQ=255)
                                                                                                    ggAAGGCTCCCGTTGGTTAACCGGGAGTAAACGTGGTTAGGCGACGTTCTGTTGACCGGCAACCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCTCCGACAGTTTGCACACCGGGCGTAGCGACTGAATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCCCTTTGCCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCtt                                                                                                                                        <  2:207512/251‑1 (MQ=255)
                                                                                                      aaGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATc                                                                                                                                                                                       <  2:10095/202‑1 (MQ=255)
                                                                                                      aaGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATc                                                                                                                                                                                       >  1:10095/1‑202 (MQ=255)
                                                                                                        ggCTCCCGGTGGTTAACCGGGATTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCCCTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAg                                                                                                                                    <  1:624/251‑1 (MQ=255)
                                                                                                                                aaaCGCTGTTACGCGACTTTCTGTTTCCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACCCCGGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCCCTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTg                                                                                                                                             <  1:423864/218‑1 (MQ=255)
                                                                                                                                aaaCGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACCCCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTg                                                                                                                                             <  1:264848/218‑1 (MQ=255)
                                                                                                                                aaaCGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAACCCGGCGGCTTCCAGCCACTg                                                                                                                                             >  2:264848/1‑218 (MQ=255)
                                                                                                                                aaaCGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTg                                                                                                                                             >  2:423864/1‑218 (MQ=255)
                                                                                                                                          aCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAGGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCAGCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGGGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTa                                                                                                   <  2:130177/250‑1 (MQ=255)
                                                                                                                                                   tctttttACCGGCAATCACCCCAAAAAACGCCTGCACCTGCTTTTTTTTCCGCTCCGACAGTTGGCACACCTGGCGTAGCGACTGCATCATTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGTCTGTAATCCGGGTAGCGACTCGc                                                                                          <  1:230421/246‑1 (MQ=255)
                                                                                                                                                      ttttcccGGCAATCCCTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGTCAGTGTGCGCACCTGGCGTAGCGACTGCATCAGTTCCCTCTCTTCGGCTGCGGGTGGTTGGGCGGTGAGGACAATACAGCTTTCCATCTCTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCAt                                                                                        <  1:307507/244‑1 (MQ=255)
                                                                                                                                                      gtgtACCGGAAATCACTCCAAGAAACCCCTGCACCTGCGTTGGTTTACGCGCCGCAAGTTTGCCCACCAGGCGTTGCGACGGCATGAGTACTCTCTCTACGGCGGCGGGTGGTTGGGCGGTGGGGAAAATACAGCCTTCCATCTCTTTGCCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCGTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCAt                                                                                        <  1:391125/246‑1 (MQ=255)
                                                                                                                                                      gTTTACCGGCCATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTCGCGTAGCGACTGCCTCACTTCGCTCTCATCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATATACCGCCGTGCCAGTCGCAAAACCGGCGGCTTCCAGCCCCTGACCTTTCAGGGAGATGGCGGGAATTAGGCTGAAATCCGGGTAGTGACTCTCATa                                                                                       >  1:287302/1‑250 (MQ=255)
                                                                                                                                                               cAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACcggcgggg                                                                             >  2:347399/1‑249 (MQ=255)
                                                                                                                                                                      ccaaaaaaccccttcacccTGGTTTTGTTTCCGGGCCGACAGTTGGCACGCCGGGCGGGGCGACTGCATAAGTTCGCTCTCTTCGGGGGCGGGTGGTGGGGCGGTGAGGCAAAAACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGAc                                                                      <  2:112184/234‑1 (MQ=255)
                                                                                                                                                                          aTAACCGCCTGCACCTGTTTTTGTTTACGGTCCGCAAGTTTGCCCACCGGGCGTAGCTCCGCCTTCAGTTCGCTCTTTCCGGCGGCGGTTGGTTGGGCTGTGAGGACTATACAGCCTCCCATCCCTTTGACATTTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGTTTCGCATAACCGACGGTGACATGTCGGtt                                                                  <  2:222761/251‑1 (MQ=255)
                                                                                                                                                                                   tGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCCGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGGGACATGACGGTTATTTGCCgg                                                         >  1:22040/1‑251 (MQ=255)
                                                                                                                                                                                        cTGCTTTTGTTTACGCGCCGACCGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGCATACGGCTGTAATCCGGGGAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCgggg                                                       >  1:130790/1‑248 (MQ=255)
                                                                                                                                                                                                             cAGTTTGCACACCGGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGGGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTtgtg                               <  1:294924/251‑1 (MQ=255)
                                                                                                                                                                                                                  tgcccacCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGTTTGGGGGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGTCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCggg                                                        <  1:431739/217‑1 (MQ=255)
                                                                                                                                                                                                                  tGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCATTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGAGGGCATACGGCTGTAATCCGGGTAGCGACTCCCATAACCGACGGTGGCATGACGGTTATTTGGCggg                                                        >  2:431739/1‑221 (MQ=255)
                                                                                                                                                                                                                                  ggACTGCAGCATTTCGCTCTCTTCGGCGTGGGGGGTTTGGGCGGTGAGGACAATACAGCCTCCCACCACGTTGACATCGACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACGTCTGCTTCGAACGGTTGTGCAATAGAATGCGTgtcagtca          <  2:197031/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aTCAGTTCGCTCTCTTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTCGCGACTCGCATAACCGAAGGTGGCCTGCCGGTTATTTGCCGTGGAGACTTCTGCTTCGAAAGGTTGTTCAATAGAATGCGTGTCCGTCATAACTGc   >  2:133008/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ccatttCGCTCTCTGCGGCGGCGGGTGGTGGGGCGGTGAGGACAATACAGCCTCCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCt  <  2:293233/247‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
CCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCT  >  NC_000913/4579363‑4579847

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