Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,973,627 T→G K566N (AAA→AAC cheA ← fused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,973,6270TG97.8% 162.2 / NA 46K566N (AAA→AACcheAfused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (27/18);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGAT  >  NC_000913/1973404‑1973854
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
gCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTcgacg                                                                                                                                                                                                                      >  2:1026423/1‑239 (MQ=255)
        gagaCATCAACAAGTGGTCCCACGCTGCCGTCGCCAAGAGTGGGCGCACGAGAACGGCCGGGGACTTTGCGATTTTTTTTTTCATGGTTTTTAACCACCACCTTGTGTTGACAAATTAATGACTCCACCAGCAAGGCGTGGCGGCACCCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCa                                                                                                                                                                                                               <  2:774458/238‑1 (MQ=255)
                aacaaTCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGAACGTGAACAAGTGCCACAGTTCGACGACGGGGAGATATTCa                                                                                                                                                                                                       >  2:681179/1‑238 (MQ=255)
                       aGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGc                                                                                                                                                                                                 >  2:28990/1‑237 (MQ=255)
                                             gAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACCGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGccgccg                                                                                                                                                                         >  1:625672/1‑239 (MQ=255)
                                                                                     aCTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGAt                                                                                                                                  >  2:1194862/1‑238 (MQ=255)
                                                                                     aCTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGGGAt                                                                                                                                  >  1:469836/1‑238 (MQ=255)
                                                                                                  taacccACACCCTGGTGTTGACAAATTAATGGTACCCACAGCAAGGCGTGGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCAtt                                                                                                                    <  1:1056928/235‑1 (MQ=255)
                                                                                                     cccACACCCGTGTGTGCACCAATAAATTGCCCCCCCAGAAGGGGGCGGGGGCGCCCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACg                                                                                                                                               <  2:710819/208‑1 (MQ=255)
                                                                                                         caACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGc                                                                                                             >  1:737617/1‑239 (MQ=255)
                                                                                                           tcttGGTTTTGACAAATAAATGACCCCACCAGCAGGGGGTGGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACg                                                                                                                                               <  2:817500/200‑1 (MQ=255)
                                                                                                            ccTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGtggatgga                                                                                                                                                            >  1:601646/1‑189 (MQ=255)
                                                                                                               ggTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGt                                                                                                                                     >  1:710819/1‑209 (MQ=255)
                                                                                                                      aCCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTcgacg                                                                                                                                                                                                                      >  1:144145/1‑121 (MQ=255)
                                                                                                                      aCCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGt                                                                                                                                     >  1:817500/1‑202 (MQ=255)
                                                                                                                      aCCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGGGATTCCATAACAGCATTCAGCGGCAGAATGaaa                                                                                                    >  1:841475/1‑235 (MQ=255)
                                                                                                                                   tCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAAcgcg                                                                                   >  1:1016311/1‑239 (MQ=255)
                                                                                                                                               gCGTGGCGGCACCCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATg                                                                         <  2:1000800/237‑1 (MQ=255)
                                                                                                                                                  tAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTc                                                                    >  1:878901/1‑239 (MQ=255)
                                                                                                                                                   aGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCg                                                                   >  1:1000800/1‑239 (MQ=255)
                                                                                                                                                        cgACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCca                                                                                                                                                                                                                                <  1:507576/77‑1 (MQ=255)
                                                                                                                                                        cgACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCca                                                                                                                                                                                                                                >  2:507576/1‑77 (MQ=255)
                                                                                                                                                          aCCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTcggcg                                                                                                                                                                                                                      >  2:674315/1‑85 (MQ=255)
                                                                                                                                                          aCCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTcggcg                                                                                                                                                                                                                      <  1:674315/85‑1 (MQ=255)
                                                                                                                                                          aCCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGaaa                                                                                                    >  1:222217/1‑199 (MQ=255)
                                                                                                                                                          aCCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGaaa                                                                                                    <  2:222217/199‑1 (MQ=255)
                                                                                                                                                                       cataTCACCACAATTCCCTGGGGTGCTTAGGTTTTCGGGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCCCCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCgg                                               <  2:429797/236‑1 (MQ=255)
                                                                                                                                                                           tccacacaATTCCCTGGGTGGCTTCTGTTTTCGCGCCCGCGACGTTGAACCCTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCa                                             <  2:589956/233‑1 (MQ=255)
                                                                                                                                                                              ccacAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGACAGCGTCATCGGCAGTaaa                                        >  1:984085/1‑239 (MQ=255)
                                                                                                                                                                                         tGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATAtt                                                                                                                                                                                                         >  2:964936/1‑67 (MQ=255)
                                                                                                                                                                                         tGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATAtt                                                                                                                                                                                                         <  1:964936/67‑1 (MQ=255)
                                                                                                                                                                                         tGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGTAAAATGCGGATCGt                             <  2:1016311/239‑1 (MQ=255)
                                                                                                                                                                                         tGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGGAAAATGCggcacgt                             >  2:1070148/1‑234 (MQ=255)
                                                                                                                                                                                            gTGGCTTCGGTTTTCGCGCCCGCGCCGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCCCCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGTAAAATGCGGATCGTAg                           <  1:517722/238‑1 (MQ=255)
                                                                                                                                                                                            gTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTc                                                                                                                                >  1:566589/1‑137 (MQ=255)
                                                                                                                                                                                            gTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTc                                                                                                                                <  2:566589/137‑1 (MQ=255)
                                                                                                                                                                                            gTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGTAAAATGCGGATCGTAg                           <  1:1194862/238‑1 (MQ=255)
                                                                                                                                                                                            gTGGCTTCGGTTTGCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGTAAAATGCGGATCGTAg                           <  1:722214/238‑1 (MQ=255)
                                                                                                                                                                                               gCTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCca                                                                                                                                                                                                                                >  2:826056/1‑38 (MQ=255)
                                                                                                                                                                                                cTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGTAAAATGCGGATCGTAGTGCCa                      >  1:733932/1‑239 (MQ=255)
                                                                                                                                                                                                cTTCGGTTTTCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGGAAAAAGCGGATCCGAGTGCCa                      >  2:357625/1‑239 (MQ=255)
                                                                                                                                                                                                         tCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACa                                                                           <  1:545158/177‑1 (MQ=255)
                                                                                                                                                                                                         tCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACa                                                                           >  2:545158/1‑177 (MQ=255)
                                                                                                                                                                                                         tCGCGCCCGCGACGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCCGCGGCACCGGCAGTAACATGCGGATCGTAGTGCCCGTACCCt               >  2:304133/1‑237 (MQ=255)
                                                                                                                                                                                                                cgcgACGTTGAACACGTTCCCAAGTTCGACGATGGGCAGATATTCCCCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACt      <  2:737617/239‑1 (MQ=255)
                                                                                                                                                                                                                     cGTTGAACACGTTCCACAGTTCGACGATGGGCAGATATTCCCCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGAt  <  2:878901/238‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
GCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCTTCACGGGGTTGCAGTGATTCCATAACAGCATTCAGCGGCAGAATGAAAACTTCATCCGCAACGCGTACGGACATGCCGTCGAGGATGGCCAGCGTCAGCGGCAGTAAAATGCGGATCGTAGTGCCAGTACCCTGCTTCGACTGGAT  >  NC_000913/1973404‑1973854

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