Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,639,623 T→C 8.2% intergenic (‑171/‑140) czcR‑III ← / → PP_1439 response regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,639,6230TC 8.2% 240.1 / 8.4 85intergenic (‑171/‑140)czcR‑III/PP_1439response regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/36);  new base C (3/4);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

GGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGA  >  NC_002947/1639474‑1639772
                                                                                                                                                     |                                                                                                                                                     
ggattgtaTTTTTGACCCTGTAGTGTGCTCTGGGTCAAGTTTCCCAATGTTACGTCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGt                                                                                                                                                       <  1:1368001/143‑1 (MQ=255)
   gACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAgg                                                                                                                                                    >  1:1533254/1‑150 (MQ=255)
      gATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAAc                                                                                                                                                 <  2:1533254/150‑1 (MQ=255)
      gATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTACGCGAGTAAACCCGCGCCTGTAGGAACg                                                                                                                                                >  1:672454/1‑150 (MQ=255)
       aTTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACg                                                                                                                                                >  2:1664543/1‑150 (MQ=255)
       aTTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACg                                                                                                                                                >  2:1665064/1‑150 (MQ=255)
           ttGACCCTGATGTCACTACAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGAATTGTAATTCTGCTGTTATCTGGCTGTCTTTTCTGGCCTCTTTGCGAGTAAACCCGCGCCTGTAGGAACgggca                                                                                                                                            >  1:314722/1‑148 (MQ=255)
                ccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGt                                                                                                                                                       >  2:1119747/1‑134 (MQ=255)
                ccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGt                                                                                                                                                       <  1:1119747/134‑1 (MQ=255)
                ccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTAccc                                                                                                                                       <  1:72001/150‑1 (MQ=255)
                 cTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcg                                                                                                                                      >  2:1455216/1‑150 (MQ=255)
                    tAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGa                                                                                                                                   >  2:1006162/1‑150 (MQ=255)
                       tGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACggg                                                                                                                                              <  1:643540/136‑1 (MQ=255)
                       tGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACggg                                                                                                                                              >  2:643540/1‑136 (MQ=255)
                        gACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGaa                                                                                                                                  <  1:1012946/147‑1 (MQ=255)
                        gACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGaa                                                                                                                                  <  1:1013343/147‑1 (MQ=255)
                        gACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGaa                                                                                                                                  >  2:1013343/1‑147 (MQ=255)
                        gACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGaa                                                                                                                                  >  2:1012946/1‑147 (MQ=255)
                            tcacagTCAAGCGCCCAAATGTTACGGAATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTTGCTGTCTGTTCGGGCCTCTGCGCGAGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGccc                                                                                                                           <  2:771667/145‑1 (MQ=21)
                                  tCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGc                                                                                                                             <  2:787221/142‑1 (MQ=255)
                                  tCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGc                                                                                                                             >  1:787221/1‑142 (MQ=255)
                                  tCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTCAGCGGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAg                                                                                                                     <  1:510074/150‑1 (MQ=255)
                                   cAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACTCGCGCTTGTAGGAACGGGTTTACCCGCGTAAGGGCCCGAGCAGt                                                                                                                    >  2:422202/1‑150 (MQ=255)
                                    aaGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                                   <  1:753326/150‑1 (MQ=255)
                                        gcCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAgacagcc                                                                                                               <  1:1142970/150‑7 (MQ=18)
                                         ccccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacac                                                                                                              <  1:1438216/150‑1 (MQ=255)
                                          cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGCTAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTGCGCGAGGAAACCCGCGCCTTTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacaca                                                                                                             >  2:709550/1‑150 (MQ=255)
                                            cAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAg                                                                                                           <  2:1688320/150‑1 (MQ=255)
                                            cAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAg                                                                                                           >  2:490779/1‑150 (MQ=255)
                                            cAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAg                                                                                                           >  2:843086/1‑150 (MQ=255)
                                             aaTGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAgg                                                                                                          <  1:843086/150‑1 (MQ=255)
                                                 ttACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCcac                                                                                                      >  1:867587/1‑150 (MQ=255)
                                                 ttACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCcac                                                                                                      >  1:651154/1‑150 (MQ=255)
                                                   aCGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                                   <  1:75374/135‑1 (MQ=255)
                                                   aCGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                                   >  2:75374/1‑135 (MQ=255)
                                                   aCGGCATGAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTACGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACAtt                                                                                                   <  2:672454/150‑1 (MQ=255)
                                                    cGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACAtt                                                                                                   <  1:1455216/150‑1 (MQ=255)
                                                       cATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGGTTACTCGCGTAGGGGCCCGTAAGGCACAGACAAGGCCACCTTGCa                                                                                                >  2:827338/1‑150 (MQ=255)
                                                          taaaaaaGGCCTGCGGTGCAGATAACGGTTTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGa                                                                                             <  2:1656756/149‑1 (MQ=255)
                                                             aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAg                                                                                                                                                     <  1:787544/91‑1 (MQ=255)
                                                             aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAg                                                                                                                                                     >  2:787544/1‑91 (MQ=255)
                                                             aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTg                                                                                          >  2:122065/1‑150 (MQ=255)
                                                             aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTg                                                                                          >  2:581185/1‑150 (MQ=255)
                                                               ttttccTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGGGGCCCGAAAAGTCCACACAAGGCCACATTGCAGGACTGTa                                                                                        >  1:1037626/5‑150 (MQ=255)
                                                                    cTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATccc                                                                                   >  2:1032294/1‑150 (MQ=255)
                                                                       cGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGc                                                                                <  2:651154/150‑1 (MQ=255)
                                                                              gATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAgaca                                                                                                                  >  1:382374/1‑106 (MQ=14)
                                                                              gATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAgaca                                                                                                                  <  2:382374/109‑4 (MQ=14)
                                                                                  aCGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGcc                                                                     <  1:1720819/150‑1 (MQ=255)
                                                                                   cGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGAGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCGCCTTGCCg                                                                    >  2:470301/1‑150 (MQ=255)
                                                                                      aTTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTa                                                                 <  1:1032294/150‑1 (MQ=255)
                                                                                           tAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagccag                                                            <  2:867587/150‑1 (MQ=255)
                                                                                           tAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACACGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACGTTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagccag                                                            >  2:1326164/1‑150 (MQ=255)
                                                                                               tCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCt                                                        >  2:1929341/1‑150 (MQ=255)
                                                                                               tCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGACTGTAGTAACGGGTTTACCAGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCt                                                        >  2:228364/1‑150 (MQ=255)
                                                                                                ctgctgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTc                                                       <  1:581185/150‑1 (MQ=255)
                                                                                                ctgctgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTc                                                       <  1:470301/150‑1 (MQ=255)
                                                                                                     gTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAg                                                  <  1:1929341/150‑1 (MQ=255)
                                                                                                     gTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTACCCAGCCCTCGCTAg                                                  >  1:305217/1‑150 (MQ=255)
                                                                                                      ttATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAgg                                                 >  2:70908/1‑150 (MQ=255)
                                                                                                         tCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATg                                              <  2:305217/150‑1 (MQ=255)
                                                                                                           tGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCTCTAGGATGGc                                            <  1:1664543/150‑1 (MQ=255)
                                                                                                           tGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGc                                            <  1:1665064/150‑1 (MQ=255)
                                                                                                             gtggactgtTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagcca                                                             >  1:1835750/6‑131 (MQ=14)
                                                                                                             gtggactgtTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagcca                                                             <  2:1835750/126‑1 (MQ=14)
                                                                                                                   tgtTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTg                                    >  2:988076/1‑150 (MQ=21)
                                                                                                                   tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagc                                                               >  1:486121/1‑123 (MQ=255)
                                                                                                                   tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagc                                                               <  2:486121/123‑1 (MQ=255)
                                                                                                                    gtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTGCCCGCGAAGAGGCCCGAACAGGCCACACAAGGCCACATTGCAGGACTGTAATCCCGGCCACACCTGGCCGTTAGGCAGCCCTCCCTAGGAGGGCCGTGTCTGa                                   >  2:69604/1‑150 (MQ=255)
                                                                                                                     ttCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgtgt                                       >  2:1837663/1‑145 (MQ=255)
                                                                                                                     ttCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgtgt                                       <  1:1837663/145‑1 (MQ=255)
                                                                                                                      tCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATc                                 <  1:490779/150‑1 (MQ=255)
                                                                                                                        gggCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGa                                   <  1:1519510/146‑1 (MQ=255)
                                                                                                                        gggCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGa                                   >  2:1519510/1‑146 (MQ=255)
                                                                                                                        gggCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCgg                               >  1:368062/1‑150 (MQ=255)
                                                                                                                         ggCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCt                                                                              <  2:1590590/102‑1 (MQ=255)
                                                                                                                         ggCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCt                                                                              >  1:1590590/1‑102 (MQ=255)
                                                                                                                          gCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGtt                             >  2:724313/1‑150 (MQ=255)
                                                                                                                            ctctTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTGCCCGCGAAGAGGACCGAACAGCCCTCAACGTTGCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttat                           <  1:1006162/150‑1 (MQ=255)
                                                                                                                             tctTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttatt                          >  1:1212816/1‑150 (MQ=255)
                                                                                                                               ttCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAt                                  <  1:1644917/140‑1 (MQ=255)
                                                                                                                               ttCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAt                                  >  2:1644917/1‑140 (MQ=255)
                                                                                                                                 cgcgAGTAAACCCGCGCTTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGa                      <  1:827338/150‑1 (MQ=255)
                                                                                                                                 cgcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGa                      <  1:1731677/150‑1 (MQ=255)
                                                                                                                                    gAGTAGACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc                   >  1:314925/1‑150 (MQ=255)
                                                                                                                                          aCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac             >  2:1899752/1‑150 (MQ=255)
                                                                                                                                                 ccTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAAc      <  1:1899752/150‑1 (MQ=255)
                                                                                                                                                   tGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCt    >  1:618871/1‑150 (MQ=255)
                                                                                                                                                    ttAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGGAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTg   <  1:709550/149‑1 (MQ=255)
                                                                                                                                                     tAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGGCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGa  >  1:1612895/1‑150 (MQ=255)
                                                                                                                                                     tAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGa  <  1:70908/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGA  >  NC_002947/1639474‑1639772

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