Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,306,969 T→A 19.0% N23Y (AAC→TAC)  phnM ← ribophosphonate triphosphate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,306,9690TA19.0% 858.7 / 191.2 315N23Y (AAC→TAC) phnMribophosphonate triphosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base T (116/139);  new base A (32/28);  total (148/167)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 320 total aligned reads displayed.

GGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATGGCTGCG  >  CP009273/4306821‑4307116
                                                                                                                                                    |                                                                                                                                                   
ggCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGtt                                                                                                                                                     >  1:3883259/1‑149 (MQ=255)
  ccAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCt                                                                                                                                                   <  2:262012/149‑1 (MQ=255)
   cAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTg                                                                                                                                                  >  2:510861/1‑149 (MQ=255)
    aGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGc                                                                                                                                                 <  2:69938/149‑1 (MQ=255)
    aGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGc                                                                                                                                                 <  1:2453684/149‑1 (MQ=255)
         aCTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTc                                                                                                                                            >  2:3125247/1‑149 (MQ=255)
           tttCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCa                                                                                                                                          >  1:607224/1‑149 (MQ=255)
                   gcgGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCg                                                                                                                                  >  2:4652763/1‑149 (MQ=255)
                   gcgGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCg                                                                                                                                  >  1:3485957/1‑149 (MQ=255)
                         tGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTacc                                                                                                                            >  1:1675745/1‑149 (MQ=255)
                          gaagaaTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTacca                                                                                                                           >  2:615739/1‑149 (MQ=255)
                          gaagaaTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTacca                                                                                                                           >  2:3424586/1‑149 (MQ=255)
                             gaaTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCt                                                                                                                        >  1:4704622/1‑149 (MQ=255)
                             gaaTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCt                                                                                                                        >  2:591035/1‑149 (MQ=255)
                               aTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCg                                                                                                                      >  2:2870314/1‑149 (MQ=255)
                                 ttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGtt                                                                                                                    >  2:4409482/1‑149 (MQ=255)
                                 ttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGtt                                                                                                                    <  1:2293322/149‑1 (MQ=255)
                                 ttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGtt                                                                                                                    <  2:357027/149‑1 (MQ=255)
                                  tATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGttt                                                                                                                   <  1:4049482/149‑1 (MQ=255)
                                  tATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGttt                                                                                                                   <  1:916033/149‑1 (MQ=255)
                                     ccAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCTTTTTcc                                                                                                                <  1:1882569/149‑1 (MQ=255)
                                     ccAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTcc                                                                                                                >  2:3278535/1‑149 (MQ=255)
                                     ccAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTcc                                                                                                                <  1:1156215/149‑1 (MQ=255)
                                        gATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGc                                                                                                             >  2:1655809/1‑149 (MQ=255)
                                         aTTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCa                                                                                                            >  1:2501479/1‑149 (MQ=255)
                                         aTTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCa                                                                                                            >  1:2966564/1‑149 (MQ=255)
                                         aTTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCa                                                                                                            >  1:1419459/1‑149 (MQ=255)
                                             tCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAg                                                                                                        >  2:3903134/1‑149 (MQ=255)
                                              cGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGc                                                                                                       >  2:4025964/1‑149 (MQ=255)
                                              cGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGc                                                                                                       >  2:4420273/1‑149 (MQ=255)
                                                gTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCtt                                                                                                     <  1:3296291/149‑1 (MQ=255)
                                                 tATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTa                                                                                                    >  1:3594747/1‑149 (MQ=255)
                                                   tGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAAc                                                                                                  <  2:322354/149‑1 (MQ=255)
                                                   tGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAAc                                                                                                  >  2:4512613/1‑149 (MQ=255)
                                                     cAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGt                                                                                                <  2:506848/149‑1 (MQ=255)
                                                     cAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGt                                                                                                <  1:679328/149‑1 (MQ=255)
                                                       gCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTa                                                                                              <  2:3193659/149‑1 (MQ=255)
                                                       gCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTa                                                                                              <  1:3347132/149‑1 (MQ=255)
                                                        cTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTAt                                                                                             <  2:894793/149‑1 (MQ=255)
                                                          cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTg                                                                                           >  2:22368/1‑149 (MQ=255)
                                                           aaTCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGa                                                                                          >  1:5082216/1‑149 (MQ=255)
                                                           aaTCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGa                                                                                          <  1:3988750/149‑1 (MQ=255)
                                                              cAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATaa                                                                                       <  1:3110774/149‑1 (MQ=255)
                                                              cAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATaa                                                                                       <  1:1188447/149‑1 (MQ=255)
                                                                gCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATc                                                                                     <  1:1817281/149‑1 (MQ=255)
                                                                 ccccGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCa                                                                                    >  2:3885392/1‑149 (MQ=255)
                                                                 ccccGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCa                                                                                    >  1:5100616/1‑149 (MQ=255)
                                                                   ccGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATg                                                                                  <  2:49030/149‑1 (MQ=255)
                                                                    cGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGa                                                                                 <  1:3036436/149‑1 (MQ=255)
                                                                    cGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGa                                                                                 <  2:1862464/149‑1 (MQ=255)
                                                                     ggCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGaa                                                                                <  2:4734505/149‑1 (MQ=255)
                                                                     ggCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGaa                                                                                <  2:2016140/149‑1 (MQ=255)
                                                                       cagcagCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAaga                                                                              <  1:5002238/149‑1 (MQ=255)
                                                                        agcagcCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAagag                                                                             <  1:3541374/149‑1 (MQ=255)
                                                                         gcagcCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAgg                                                                            >  2:4446933/1‑149 (MQ=255)
                                                                          cagccagcCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGc                                                                           <  2:1093493/149‑1 (MQ=255)
                                                                          cagccagcCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGc                                                                           <  2:1406882/149‑1 (MQ=255)
                                                                          cagccagcCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGc                                                                           <  2:4180634/149‑1 (MQ=255)
                                                                            gccagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTc                                                                         >  1:3160848/1‑149 (MQ=255)
                                                                             ccagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTcc                                                                        <  2:1352712/149‑1 (MQ=255)
                                                                             ccagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTcc                                                                        >  1:493477/1‑149 (MQ=255)
                                                                              cagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTccc                                                                       <  2:4771387/149‑1 (MQ=255)
                                                                               agccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCa                                                                      <  2:4094013/149‑1 (MQ=255)
                                                                                   gccTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTACCGTTATTGATAATCATGAAGAGGCTCCCATTgg                                                                  >  2:2475207/1‑149 (MQ=255)
                                                                                   gccTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTgg                                                                  <  1:1234437/149‑1 (MQ=255)
                                                                                     cTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGt                                                                >  2:72808/1‑149 (MQ=255)
                                                                                        cGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCa                                                             <  1:588236/149‑1 (MQ=255)
                                                                                        cGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCa                                                             <  2:4625431/149‑1 (MQ=255)
                                                                                            gTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCg                                                         >  1:2007669/1‑149 (MQ=255)
                                                                                            gTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTCTTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCg                                                         >  2:4272866/1‑149 (MQ=255)
                                                                                              ccATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGt                                                       <  1:4749036/149‑1 (MQ=255)
                                                                                               cATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTc                                                      <  2:2400835/149‑1 (MQ=255)
                                                                                               cATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTc                                                      <  1:4846292/149‑1 (MQ=255)
                                                                                                aTCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCg                                                     <  1:4592208/149‑1 (MQ=255)
                                                                                                  cGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGc                                                   <  2:1567712/149‑1 (MQ=255)
                                                                                                  cGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGc                                                   <  2:3040273/149‑1 (MQ=255)
                                                                                                  cGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGc                                                   <  1:4906944/149‑1 (MQ=255)
                                                                                                    ccTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGa                                                 <  2:4819536/149‑1 (MQ=255)
                                                                                                         ggCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCa                                            >  2:2955099/1‑149 (MQ=255)
                                                                                                          gCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCAt                                           <  1:273650/149‑1 (MQ=255)
                                                                                                           cAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCAtt                                          <  1:1747266/149‑1 (MQ=255)
                                                                                                             ggcggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTAc                                        <  2:2897867/149‑1 (MQ=255)
                                                                                                             ggcggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTAc                                        >  2:923736/1‑149 (MQ=255)
                                                                                                                ggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTa                                     <  2:1233323/149‑1 (MQ=255)
                                                                                                                    ctGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGc                                 >  1:2516296/1‑149 (MQ=255)
                                                                                                                    ctGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGc                                 >  1:4196258/1‑149 (MQ=255)
                                                                                                                      ggCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCtt                               <  2:4328082/149‑1 (MQ=255)
                                                                                                                      ggCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCtt                               <  1:4465758/149‑1 (MQ=255)
                                                                                                                       gCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTc                              >  2:68005/1‑149 (MQ=255)
                                                                                                                                aaaGGCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGATTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGaa                     <  2:2830309/149‑1 (MQ=255)
                                                                                                                                 aaGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAg                    >  2:3608157/1‑149 (MQ=255)
                                                                                                                                  aGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGa                   >  1:1055000/1‑149 (MQ=255)
                                                                                                                                   ggCGCGGATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGAt                  >  1:4070207/1‑149 (MQ=255)
                                                                                                                                    gcgcGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATg                 <  2:2321561/149‑1 (MQ=255)
                                                                                                                                        ggATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTa             <  1:545389/149‑1 (MQ=255)
                                                                                                                                         gATTTCGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTAc            >  1:378734/1‑149 (MQ=255)
                                                                                                                                              cGCCGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATgg       <  2:3756533/149‑1 (MQ=255)
                                                                                                                                               gCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATGGc      >  2:4220893/1‑149 (MQ=255)
                                                                                                                                                 cGTACTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATCATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATGGCTg    <  1:3868630/149‑1 (MQ=255)
                                                                                                                                                  gTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATGGCTgc   <  2:2740790/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTCCCATTGGGTGCAGGCGGTCGGCGACGTCATTACGTACAGCTTCGTCATGGAAGATGCCTACGATGGCTGCG  >  CP009273/4306821‑4307116

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