Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,968,073 T→C 5.0% A143A (GCA→GCG) ‡ cheB ← fused chemotaxis regulator: protein‑glutamate methylesterase in two‑component regulatory system with CheA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,968,0730TC 5.0% 341.1 / 4.2 96A143A (GCA→GCG) ‡cheBfused chemotaxis regulator: protein‑glutamate methylesterase in two‑component regulatory system with CheA
Reads supporting (aligned to +/- strand):  ref base T (46/45);  new base C (2/3);  total (48/48)
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

AATTAACAGTGCCGGGCTGGAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGC  >  NC_000913/1967941‑1968214
                                                                                                                                    |                                                                                                                                             
aaTTAACAGTGCCGGGCTGGAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAAt                                                                                                                                      <  2:1052817/142‑1 (MQ=255)
aaTTAACAGTGCCGGGCTGGAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAAt                                                                                                                                      >  1:1052817/1‑142 (MQ=255)
     aCAGTGCCGGGCTGGAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACaa                                                                                                                                       >  4:234393/1‑136 (MQ=255)
     aCAGTGCCGGGCTGGAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACaa                                                                                                                                       <  3:234393/136‑1 (MQ=255)
          gCCGGGCTGGAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCa                                                                                                                       >  2:532062/1‑147 (MQ=255)
          gCCGGGCTGGAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCa                                                                                                                       <  1:532062/147‑1 (MQ=255)
            cGGGCTGGAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAgg                                                                                                                    <  2:116756/148‑1 (MQ=255)
                  ggAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGcc                                                                                                              >  4:277101/1‑148 (MQ=255)
                           cAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGcc                                                                                                     <  3:277101/148‑1 (MQ=255)
                             aCGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTAtgtg                                                                                                                            >  2:158180/1‑123 (MQ=255)
                             aCGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTAtgtg                                                                                                                            <  1:158180/123‑1 (MQ=255)
                             aCGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGt                                                                                                   <  1:714086/148‑1 (MQ=255)
                             aCGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGt                                                                                                   >  2:714086/1‑148 (MQ=255)
                                     aGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCt                                                                                           >  4:714332/1‑148 (MQ=255)
                                      gTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCtt                                                                                          >  3:376229/1‑148 (MQ=255)
                                            gACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAAt                                                                                                                                      >  4:590833/1‑98 (MQ=255)
                                            gACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAAt                                                                                                                                      <  3:590833/98‑1 (MQ=255)
                                             aCGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATgg                                                                                                                                    <  1:891374/99‑1 (MQ=255)
                                             aCGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATgg                                                                                                                                    >  2:891374/1‑99 (MQ=255)
                                             aCGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCa                                                                                   <  1:572816/148‑1 (MQ=255)
                                             aCGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCa                                                                                   >  2:572816/1‑148 (MQ=255)
                                               gAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAAt                                                                                 <  4:376229/148‑1 (MQ=255)
                                                aaTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCa                                                                                                                       >  1:227388/1‑109 (MQ=255)
                                                aaTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCa                                                                                                                       <  2:227388/109‑1 (MQ=255)
                                                aaTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCg                                                                                                                >  3:315681/1‑116 (MQ=255)
                                                aaTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCg                                                                                                                <  4:315681/116‑1 (MQ=255)
                                                aaTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCtttt                                                                                        >  2:913998/1‑140 (MQ=255)
                                                aaTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCtttt                                                                                        <  1:913998/140‑1 (MQ=255)
                                                 aTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTa                                                                                                  <  3:471509/129‑1 (MQ=255)
                                                 aTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTa                                                                                                  >  4:471509/1‑129 (MQ=255)
                                                 aTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTa                                                                                                  >  4:58972/1‑129 (MQ=255)
                                                 aTTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTa                                                                                                  <  3:58972/129‑1 (MQ=255)
                                                        cAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTtg                                                                                                          <  1:431750/114‑1 (MQ=255)
                                                        cAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTtg                                                                                                          >  2:431750/1‑114 (MQ=255)
                                                        cAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGtt                                                                        >  2:803891/1‑148 (MQ=255)
                                                        cAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGtt                                                                        <  3:714332/148‑1 (MQ=255)
                                                                cccGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCa                                                                >  1:211595/1‑148 (MQ=255)
                                                                 ccGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTa                                                                                                                                <  4:290246/83‑1 (MQ=255)
                                                                 ccGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTa                                                                                                                                >  3:290246/1‑83 (MQ=255)
                                                                  cGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCg                                                                                                    <  1:871434/110‑1 (MQ=255)
                                                                  cGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCg                                                                                                    >  2:871434/1‑110 (MQ=255)
                                                                    ttGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCttt                                                                                         <  3:801012/119‑1 (MQ=255)
                                                                    ttGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCttt                                                                                         >  4:801012/1‑119 (MQ=255)
                                                                     tGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATa                                                           <  1:803891/148‑1 (MQ=255)
                                                                       aaGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATAcc                                                         >  3:543736/1‑148 (MQ=255)
                                                                        aGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTa                                                                                                  >  3:142392/1‑106 (MQ=255)
                                                                        aGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTa                                                                                                  <  4:142392/106‑1 (MQ=255)
                                                                         gCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCtt                                                       <  4:543736/148‑1 (MQ=255)
                                                                         gCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCtt                                                       >  1:986196/1‑148 (MQ=255)
                                                                             cAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCtt                                                       <  1:296498/144‑1 (MQ=255)
                                                                             cAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCtt                                                       >  2:296498/1‑144 (MQ=255)
                                                                              aaTCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCa                                                                <  3:238093/134‑1 (MQ=255)
                                                                              aaTCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCa                                                                >  4:238093/1‑134 (MQ=255)
                                                                                tCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAAt                                                >  1:303229/1‑148 (MQ=255)
                                                                                   cAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAAt                                                                                 >  1:41092/1‑112 (MQ=255)
                                                                                   cAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAAt                                                                                 <  2:41092/112‑1 (MQ=255)
                                                                                    aaTCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGtgctg                                                                                                                         <  3:416941/71‑1 (MQ=255)
                                                                                    aaTCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGtgctg                                                                                                                         >  4:416941/1‑71 (MQ=255)
                                                                                     aTCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATAcc                                             <  2:356147/146‑1 (MQ=255)
                                                                                     aTCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATAcc                                             >  1:356147/1‑146 (MQ=255)
                                                                                          tttttCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGc                                      <  2:303229/148‑1 (MQ=255)
                                                                                           ttttCAGAACTCAACAACGGCGCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTc                                                                                       >  1:686267/1‑98 (MQ=255)
                                                                                           ttttCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTc                                                                                       <  2:686267/98‑1 (MQ=255)
                                                                                           ttttCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCa                                                             <  3:501033/124‑1 (MQ=255)
                                                                                           ttttCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCa                                                             >  4:501033/1‑124 (MQ=255)
                                                                                           ttttCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCg                                     >  2:470567/1‑148 (MQ=255)
                                                                                            tttCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCgg                                    <  2:211595/148‑1 (MQ=255)
                                                                                                aGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTg                                <  2:986196/148‑1 (MQ=255)
                                                                                                 gAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGTCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTa                                                                                                  >  2:37245/1‑81 (MQ=255)
                                                                                                 gAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTa                                                                                                  <  1:37245/81‑1 (MQ=255)
                                                                                                 gAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGc                                                              >  1:557200/1‑117 (MQ=255)
                                                                                                 gAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGc                                                              <  2:557200/117‑1 (MQ=255)
                                                                                                    cTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTc                                                                                       <  2:574077/89‑1 (MQ=255)
                                                                                                    cTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTc                                                                                       >  1:574077/1‑89 (MQ=255)
                                                                                                      caacaaCGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAAt                                                                                 >  2:672853/1‑93 (MQ=255)
                                                                                                      caacaaCGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAAt                                                                                 <  1:672853/93‑1 (MQ=255)
                                                                                                      caacaaCGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTgg                               >  1:1085897/1‑143 (MQ=255)
                                                                                                      caacaaCGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTgg                               <  2:1085897/143‑1 (MQ=255)
                                                                                                          aacGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCtt                                                       <  1:513430/115‑1 (MQ=255)
                                                                                                          aacGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCtt                                                       >  2:513430/1‑115 (MQ=255)
                                                                                                          aacGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAAt                      <  1:470567/148‑1 (MQ=255)
                                                                                                            cGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACaaa                        <  4:693480/144‑1 (MQ=255)
                                                                                                            cGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACaaa                        >  3:693480/1‑144 (MQ=255)
                                                                                                                  cGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTa                   <  1:960473/143‑1 (MQ=255)
                                                                                                                  cGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTa                   >  2:960473/1‑143 (MQ=255)
                                                                                                                       tCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTCCGGTTTGGTGACAAAATCTATCGCCCCCAg         >  1:731719/1‑148 (MQ=255)
                                                                                                                            attGTCGGCACCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGttt                                 <  2:514367/118‑1 (MQ=255)
                                                                                                                            attGTCGGCACCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGttt                                 >  1:514367/2‑119 (MQ=255)
                                                                                                                             ttGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATa                                                           <  4:139867/92‑1 (MQ=255)
                                                                                                                             ttGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATa                                                           >  3:139867/1‑92 (MQ=255)
                                                                                                                             ttGTCGGCACCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACaaa                        <  3:362134/127‑1 (MQ=255)
                                                                                                                              tGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAgc  >  4:692360/1‑148 (MQ=255)
                                                                                                                              tGTCGGCACCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTcc     >  1:675121/1‑145 (MQ=255)
                                                                                                                              tGTCGGCACCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTcc     <  2:675121/145‑1 (MQ=255)
                                                                                                                                 cGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGtt                                  >  3:444968/1‑113 (MQ=255)
                                                                                                                                 cGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGtt                                  <  4:444968/113‑1 (MQ=255)
                                                                                                                                    tGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGa                                                  <  2:183650/94‑1 (MQ=255)
                                                                                                                                    tGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGa                                                  >  1:183650/1‑94 (MQ=255)
                                                                                                                                    |                                                                                                                                             
AATTAACAGTGCCGGGCTGGAAAGCGGCAACGGTTGCAGTACGTGACGAATTGCCTCAGTTCCACCCGTTGAAGCACCAATCGCAATCAGTTTTTCAGAACTCAACAACGGCCCCGCCTTCAGCGTTGTCGGTGCCGACAATGGCTTATGTGCTGCAAGGCTCGCCTTTGCTGCCGTACGCACCTTTTCAGCAATCATTTCGTTATACGCCAGCATACCTTCGCGAATACCCAGTTGCGGTTTGGTGACAAAATCTATCGCCCCCAGCTCCAGC  >  NC_000913/1967941‑1968214

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