Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 175,410 C→A 5.5% P102T (CCG→ACG)  clcA → chloride:H(+) antiporter ClcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913175,4100CA 5.5% 253.2 / 3.2 74P102T (CCG→ACG) clcAchloride:H(+) antiporter ClcA
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base A (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

CAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGA  >  NC_000913/175140‑175621
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                   
cAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGCTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTTCTGGCGGTGTTTGGCTACGTTTTGGTGGCCACAAACGCGcc                                                                                                                                                                                                                    >  2:475592/1‑272 (MQ=255)
          gcCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCCACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGGTCTGGCGTTGTTTGGCTACTTTTTGGTGCGCAAATACGCGcccgggagcggg                                                                                                                                                                                                          >  2:960016/1‑262 (MQ=255)
                                                  cgATAATACCCCGTTAGCCTTTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGTTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTAGGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTTTTTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTGGGCTACTTTTTGGTGCGCAAATACGCGCGGGAAGCAGGTGGTTCGGGGATCCCGAAAATTGAAGGGGCGCTGGAAGATc                                                                                                                                                                  <  1:778419/272‑1 (MQ=255)
                                                     tAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATAATGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGAGGTTTTGCTACTTTTTGGTGAGCAAATACGCGCCGGaaa                                                                                                                                                                                                               >  1:529418/1‑223 (MQ=255)
                                                        aatcccGTTAGGCATTTTGTTGATTGCGGCAGGCGTCGGCAGGCTTGGTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGGCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGGTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTc                                                                                                                                                            <  1:475593/269‑1 (MQ=255)
                                                        aaCCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTc                                                                                                                                                            >  1:829937/1‑272 (MQ=255)
                                                          cctggTTGGCCTTTTTGTTTTTGGCGGCGTTTGTGGGGACGGTTGTTGGGGTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTccc                                                                                                                                                          <  2:829936/268‑1 (MQ=255)
                                                                  gCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGTGCGCTGGAAGAACAACGGTCCGTTTGcgg                                                                                                                                                  >  2:345620/1‑270 (MQ=255)
                                                                   ccATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTATTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCtg                                                                                                                                                  <  1:691964/271‑1 (MQ=255)
                                                                                        tcgtcgGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGaa                                                                                                                                                                                                                >  1:516377/1‑188 (MQ=255)
                                                                                         cgtcgGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAg                                                                                                                                                                                                               <  2:516376/188‑1 (MQ=255)
                                                                                                     tgttgGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATGAACGTCCCGTTCGCGGGTGGCTTGTAATTGCGGTGGAGTCATTTgggcggt                                                                                                               >  2:766955/1‑267 (MQ=255)
                                                                                                          ggCTTGCAGCGGTTGTTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATTTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTgggg                                                                                                          <  1:960017/272‑1 (MQ=255)
                                                                                                               gCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTTCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCGGGTGGCGTGGATTGGCGGGGAAGTTTTTTGGCGGGCTcgggggggc                                                                                                     >  2:947214/1‑267 (MQ=255)
                                                                                                                 aGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCaa                                                                                                                                                                >  1:880347/1‑211 (MQ=255)
                                                                                                                 aGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTc                                                                                                                                                      >  2:945829/1‑221 (MQ=255)
                                                                                                                 aGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGACGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGa                                                                                                                                                                                   >  1:800910/1‑192 (MQ=255)
                                                                                                                  gCGGTTGCTTTTGACAAAGGTGTCGCCTTGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCg                                                                                                                                                     <  1:945830/221‑1 (MQ=255)
                                                                                                                  gCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAAc                                                                                                                                                               <  2:880346/211‑1 (MQ=255)
                                                                                                                  gCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGACGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGaa                                                                                                                                                                                  <  2:800909/192‑1 (MQ=255)
                                                                                                                    ggTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTccc                                                                                                                                                          <  1:15026/214‑1 (MQ=255)
                                                                                                                       ttgttttGACAAAGTTGTCGCCTTGTTGCAGAACCAACGGATGGGTGCTGTGGTACATACTTCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  1:175545/269‑1 (MQ=255)
                                                                                                                       tGCTTTTGACAAAGGTGTTGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  1:766956/272‑1 (MQ=255)
                                                                                                                       tGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTGCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  1:947215/272‑1 (MQ=255)
                                                                                                                        gCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGGATTGCCGGGGAAGTTCTTTGGCGGGGTGGGGAAACGTGGCGGAgg                                                                                            >  1:1020652/1‑272 (MQ=255)
                                                                                                                                       ttCGCCTTGTTGCAGAACCAACGGATGGGGGCGCTGGTACATACTTCTGATAATTATCCGCTTCTGTTAACCGTCTCTTTTCTCTGTTCGGCGGTGCTGGGGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGgcg                                                                             <  1:380015/271‑1 (MQ=255)
                                                                                                                                                              tATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTc                                                                                                    >  1:534814/1‑226 (MQ=255)
                                                                                                                                                                tacgtataggggTATATATTGCTGTTATTTTTCCGCTTCTGTTATCCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATcggc                                                    <  2:381165/262‑1 (MQ=255)
                                                                                                                                                                   gggCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGa                                                                                                                                                                                   >  1:321235/1‑142 (MQ=255)
                                                                                                                                                                   gggAGCTGGTACATACTGCTGATAATTATCCGCTTATGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTCTTTGTTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCGGGGGACACTCGGCGGAGGCATGGTGTTGGGGAGCGAGGGGCCAACCGTGCAAATAGGCGGt                                                 >  1:355506/1‑272 (MQ=255)
                                                                                                                                                                    ggCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGaa                                                                                                                                                                                  <  2:321234/142‑1 (MQ=255)
                                                                                                                                                                     gcgcTGGTACATACTGCTGATAATTTTCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGGTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGTGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGGTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTaa                                               <  1:808163/272‑1 (MQ=255)
                                                                                                                                                                     gcgcTGGTACATACTGCTGATAATTATCCGCGTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTaa                                               <  1:345621/272‑1 (MQ=255)
                                                                                                                                                                      acagggTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTaa                                               <  1:238060/267‑1 (MQ=255)
                                                                                                                                                                        cgtgtgtttattgctgATGTTTTTCTGTTTTTGGTAACTTTCTCTTTTCTCTGGTCGGCGGGGCTGGCGAGGTTTGGCTGCGTTGTGGTGCGCCAATACGCGCCGGAAGCAGGTGGTGCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGGTCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACAt                                            <  2:638564/261‑1 (MQ=255)
                                                                                                                                                                        cTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTa                                                                                                                                        >  1:520487/1‑180 (MQ=255)
                                                                                                                                                                         tGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTa                                                                                                                                        >  1:836215/1‑179 (MQ=255)
                                                                                                                                                                         tGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTAt                                                                                                                                       <  2:520486/180‑1 (MQ=255)
                                                                                                                                                                          ggTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTAt                                                                                                                                       <  2:836214/179‑1 (MQ=255)
                                                                                                                                                                           gTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGGGTATTGCCGGTGAAGGTCTTTGGCGGGCTGGGGACACTCGGGGGAGGCATGGTGGTGGGGCGCGCAGGGACAACCGTGCAGATCGGCGGTAACATTgg                                         >  1:775013/1‑272 (MQ=255)
                                                                                                                                                                                    gctgATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCtggt                                                                                                                                                >  1:922270/1‑160 (MQ=255)
                                                                                                                                                                                     ctgATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCtggtg                                                                                                                                               <  2:922269/160‑1 (MQ=255)
                                                                                                                                                                                      tgATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAgg                                                                                                                                                                                >  2:141069/1‑126 (MQ=255)
                                                                                                                                                                                      tgATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  2:948468/209‑1 (MQ=255)
                                                                                                                                                                                      tgATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             >  1:948469/1‑209 (MQ=255)
                                                                                                                                                                                      tgATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAGGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTg                               >  1:514448/1‑271 (MQ=255)
                                                                                                                                                                                       gATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAggg                                                                                                                                                                               <  1:141069/126‑1 (MQ=255)
                                                                                                                                                                                        aTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACgcgc                                                                                                                                                                                                                     >  2:577681/1‑87 (MQ=255)
                                                                                                                                                                                        aTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTcc                                                                                                                                                           >  2:121888/1‑145 (MQ=255)
                                                                                                                                                                                        aTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGcc                                                                                                                                   >  2:972882/1‑169 (MQ=255)
                                                                                                                                                                                        aTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTc                                                                                                                         >  1:156024/1‑179 (MQ=255)
                                                                                                                                                                                        aTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGa                                                                                              >  2:1026218/1‑206 (MQ=255)
                                                                                                                                                                                        aTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGa                                                                                              >  2:1118711/1‑206 (MQ=255)
                                                                                                                                                                                        aTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATcgg                                                     >  1:381166/1‑247 (MQ=255)
                                                                                                                                                                                        aTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGAGCGGCGGTAACa                                             >  1:638565/1‑255 (MQ=255)
                                                                                                                                                                                        aTAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACCCGGCGGAGTCGTGGTGTTGGGGCGCGAAGGGCCACACGGGCAGGACGGCGGGAAAGTAGTGCGGATGGTGCtg                            >  2:429519/1‑271 (MQ=255)
                                                                                                                                                                                         tagttATCCGCTTCTGTTAACCGTCGCTTGTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  1:36068/203‑1 (MQ=255)
                                                                                                                                                                                         tAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGcc                                                                                                                                                                                                                    <  1:577682/87‑1 (MQ=255)
                                                                                                                                                                                         tAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTccc                                                                                                                                                          <  1:121888/145‑1 (MQ=255)
                                                                                                                                                                                         tAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCg                                                                                                                                  <  1:972883/169‑1 (MQ=255)
                                                                                                                                                                                         tAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCt                                                                                                                        <  2:156024/179‑1 (MQ=255)
                                                                                                                                                                                         tAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  1:1026219/206‑1 (MQ=255)
                                                                                                                                                                                         tAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  1:1118712/206‑1 (MQ=255)
                                                                                                                                                                                         tAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGGGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGgcg                                                                             <  1:729538/222‑1 (MQ=255)
                                                                                                                                                                                         tAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTGTGGCTACTTTTTGGGGCGCAAATACGCGGCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTaa                                               <  1:322979/252‑1 (MQ=255)
                                                                                                                                                                                         tAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCGGGCGATGTTTGGCTACTGTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  1:1009065/206‑1 (MQ=255)
                                                                                                                                                                                                  gCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGAAAATACGCGACGGAAGCAGGTGGTTCGGGGATCTCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGa                                                                                              >  2:568032/1‑196 (MQ=255)
                                                                                                                                                                                                  gCTTCTGTTAACCGTCGCTTTTCTCTGTTAGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGc                                                                                                                                                    >  1:293646/1‑142 (MQ=255)
                                                                                                                                                                                                   cTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGAAAATACGCGACGGAAGCAGGTGGTTCGGGGATCTCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  1:568033/196‑1 (MQ=255)
                                                                                                                                                                                                   cTTCTGTTAACCGTCGCTTTTCTCTGTTAGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCt                                                                                                                                                   <  2:293645/142‑1 (MQ=255)
                                                                                                                                                                                                                  gctttctgttGTTCGGCGGTGCTGGCGTTGTTTGGCTACTTTTTGGTGACAAAAGACGGGCCGGAAGCAGGTGGGTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCGACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACg   <  1:429520/262‑1 (MQ=255)
                                                                                                                                                                                                                  gctgtttcttGTTCGTGGTTGCTGGGGATGTTTGGTTACTTGTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGa  <  2:775012/263‑1 (MQ=255)
                                                                                                                                                                                                                                     tGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAggg                                                                                                                                                                               >  2:111488/1‑80 (MQ=255)
                                                                                                                                                                                                                                      gCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAgggg                                                                                                                                                                              <  1:111488/80‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGa                                                                                              >  2:187904/1‑148 (MQ=255)
                                                                                                                                                                                                                                                   ggCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAg                                                                                             <  1:187904/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaTACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACa                                             >  2:1045263/1‑178 (MQ=255)
                                                                                                                                                                                                                                                                      aaTACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACAt                                            <  1:1045264/178‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                   
CAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGCTGGCAGCGGTTGCTTTTGACAAAGGTGTCGCCTGGTTGCAGAACCAACGTATGGGGGCGCTGGTACATACTGCTGATAATTATCCGCTTCTGTTAACCGTCGCTTTTCTCTGTTCGGCGGTGCTGGCGATGTTTGGCTACTTTTTGGTGCGCAAATACGCGCCGGAAGCAGGTGGTTCGGGGATCCCGGAAATTGAAGGGGCGCTGGAAGATCAACGTCCCGTTCGCTGGTGGCGTGTATTGCCGGTGAAGTTCTTTGGCGGGCTGGGGACACTCGGCGGAGGCATGGTGTTGGGGCGCGAAGGGCCAACCGTGCAGATCGGCGGTAACATTGGCCGTATGGTGCTTGATATTTTCCGCCTGAAAGGTGACGA  >  NC_000913/175140‑175621

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