Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,318,639 A→C F413C (TTT→TGT)  rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,318,6390AC98.3% 203.5 / NA 58F413C (TTT→TGT) rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (29/28);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACAT  >  NC_000913/2318406‑2318871
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        
tACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACacaa                                                                                                                                                                                                                                        >  2:110898/1‑236 (MQ=255)
     aTAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGAc                                                                                                                                                                                                                                   <  1:76673/236‑1 (MQ=255)
       aaTGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGCATTGTTGCTGTTCAGGAAATCAACACAATTGa                                                                                                                                                                                                                                    >  1:432250/1‑233 (MQ=255)
           ccATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACctgctg                                                                                                                                                                                                                             <  1:110898/236‑1 (MQ=255)
               aCAGCGTCGTGCGCATCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACATTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCca                                                                                                                                                                                                                         <  2:685076/236‑1 (MQ=255)
                   cggcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGTCATCAACACAATTGACCATCTGCCCGCAgc                                                                                                                                                                                                                     >  2:96062/1‑235 (MQ=255)
                   cggcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACctgctg                                                                                                                                                                                                                             <  1:155401/228‑1 (MQ=255)
                   cggcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACctgctg                                                                                                                                                                                                                             >  2:155401/1‑228 (MQ=255)
                   cggcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGa                                                                                                                                                                                                                     >  1:648249/1‑236 (MQ=255)
                     gcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATa                                                                                                                                                                                                                   >  1:4928/1‑236 (MQ=255)
                      cgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATaa                                                                                                                                                                                                                  >  1:578907/1‑236 (MQ=255)
                                       tgcaggTTGCAAGGAACTTCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCAGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTgg                                                                                                                                                                                                 <  2:257314/232‑1 (MQ=255)
                                         acGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGt                                                                                                                                                                                                           <  1:372459/224‑1 (MQ=255)
                                         acGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGt                                                                                                                                                                                                           >  2:372459/1‑224 (MQ=255)
                                          cGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCg                                                                                                                                                                                               <  1:173511/235‑1 (MQ=255)
                                              ggCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGAc                                                                                                                                                                                                                                   <  2:663112/195‑1 (MQ=255)
                                              ggCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGAc                                                                                                                                                                                                                                   >  1:663112/1‑195 (MQ=255)
                                                    gAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCg                                                                                                                                                                                    <  1:423311/236‑1 (MQ=255)
                                                     aaCATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCg                                                                                                                                                                                    <  1:645130/235‑1 (MQ=255)
                                                                       cTCGCCTGTTCCGCTGCTTGGGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGt                                                                                                                                                                 <  1:96062/236‑1 (MQ=255)
                                                                            cTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATgcg                                                                                                                                                            >  1:512376/1‑236 (MQ=255)
                                                                            cTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATgcg                                                                                                                                                            <  1:196021/236‑1 (MQ=255)
                                                                               ttCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCCGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCCTGCGTAAGCACATTGAGATAGGTATGCTCCa                                                                                                                                                         >  2:377229/1‑236 (MQ=255)
                                                                                    ctgctTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGTTAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATCTTGCGATAGGTATCCGCCAGTTc                                                                                                                                                     >  1:317241/1‑235 (MQ=255)
                                                                                            gCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCg                                                                                                                                                    <  2:497841/228‑1 (MQ=255)
                                                                                            gCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCg                                                                                                                                                    >  1:497841/1‑228 (MQ=255)
                                                                                                 ctcCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTa                                                                                                                                             <  2:759045/230‑1 (MQ=255)
                                                                                                 ctcCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTa                                                                                                                                             >  1:759045/1‑230 (MQ=255)
                                                                                                     ttgcaCGTCTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCCCACAAAGGGCCACGTTTTCATTCCGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGAc                                                                                                                                   <  2:322340/232‑1 (MQ=255)
                                                                                                     tGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGAc                                                                                                                                   <  2:512376/236‑1 (MQ=255)
                                                                                                          cGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGCTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGGGCCAGTTCGTTACTTAAAATATTGACGCCAt                                                                                                                              >  2:696339/1‑236 (MQ=255)
                                                                                                                      cTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCa                                                                                                                  >  2:631714/1‑236 (MQ=255)
                                                                                                                        tCACGCGCAAAGAAACATCCACCAGCACACAAATGGCCAAGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCaaa                                                                                                                <  2:648249/236‑1 (MQ=255)
                                                                                                                            gcgcgAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCat                                                                                                                                                                             >  1:349436/1‑171 (MQ=255)
                                                                                                                            gcgcgAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCat                                                                                                                                                                             <  2:349436/171‑1 (MQ=255)
                                                                                                                                                    acaAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGACTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATAACACCTGGCTCCGAGGCGa                                                                                    >  2:704582/1‑236 (MQ=255)
                                                                                                                                                        aTGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATa                                                                                                       <  1:169073/213‑1 (MQ=255)
                                                                                                                                                        aTGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATa                                                                                                       >  2:169073/1‑213 (MQ=255)
                                                                                                                                                           gCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCtt                                                                             >  1:576526/1‑236 (MQ=255)
                                                                                                                                                                    tCATTACGATAGCTCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGaa                                                                    <  2:766806/236‑1 (MQ=255)
                                                                                                                                                                    tCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGaa                                                                    <  1:156720/236‑1 (MQ=255)
                                                                                                                                                                                    aaTGGACGAAGCTGATTGTCAAATTGAGATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGttcttc                                                     <  2:576526/235‑1 (MQ=255)
                                                                                                                                                                                     aTGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCAt                                                            <  2:332013/227‑1 (MQ=255)
                                                                                                                                                                                     aTGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCAt                                                            >  1:332013/1‑227 (MQ=255)
                                                                                                                                                                                        gACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTc                                                <  1:378131/236‑1 (MQ=255)
                                                                                                                                                                                                            tGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTccgcc                            >  2:630677/1‑236 (MQ=255)
                                                                                                                                                                                                                    ttgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGaa                                                                    >  1:529499/1‑188 (MQ=255)
                                                                                                                                                                                                                    ttgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGaa                                                                    <  2:529499/188‑1 (MQ=255)
                                                                                                                                                                                                                    ttgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGaa                    <  2:4928/236‑1 (MQ=255)
                                                                                                                                                                                                                      gCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGTGCGTCGCTTTCCGCCGGATTGgaaa                  >  2:395560/1‑236 (MQ=255)
                                                                                                                                                                                                                      gCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTccgccg                           >  2:395186/1‑227 (MQ=255)
                                                                                                                                                                                                                      gCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTccgccg                           <  1:395186/227‑1 (MQ=255)
                                                                                                                                                                                                                      gCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGaaaa                  >  1:159996/1‑236 (MQ=255)
                                                                                                                                                                                                                        tggcaaggcaaTCAACCAAATTGACCGGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATa                <  2:317241/226‑1 (MQ=255)
                                                                                                                                                                                                                         ggTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCATGGCGTCGCTTTCCGCCGGAATGAAAATAc               >  1:602822/1‑236 (MQ=255)
                                                                                                                                                                                                                         ggTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGATTTCCGCCGGCATGAACATAc               >  2:748027/1‑236 (MQ=255)
                                                                                                                                                                                                                         ggTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCACGGCGTCGCTTTCCGCCGGACTGAAAATaa               >  1:332002/1‑235 (MQ=255)
                                                                                                                                                                                                                            cAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTc            <  2:578907/236‑1 (MQ=255)
                                                                                                                                                                                                                               gACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCt         >  2:236058/1‑236 (MQ=255)
                                                                                                                                                                                                                                    caacacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTa     >  2:163704/1‑235 (MQ=255)
                                                                                                                                                                                                                                    caacacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCAGCCGGAATGAAAATActcggctcgaac    >  2:758756/1‑225 (MQ=255)
                                                                                                                                                                                                                                     aacacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACa   >  1:65330/1‑236 (MQ=255)
                                                                                                                                                                                                                                      acacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACAt  >  1:438985/1‑236 (MQ=255)
                                                                                                                                                                                                                                      acacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACAt  <  1:630677/236‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        
TACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATATTGACGCCATCAGCGGTACGCAAAATGCAGATACCCACTGGCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACAT  >  NC_000913/2318406‑2318871

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