Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,257,232 G→T 87.2% R34S (CGC→AGC)  prs ← phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,257,2320GT87.2% 401.0 / 60.8 172R34S (CGC→AGC) prsphosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (13/9);  new base T (86/64);  total (99/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 172 total aligned reads displayed.

GCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCC  >  CP009273/1257087‑1257375
                                                                                                                                                 |                                                                                                                                               
   gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATCTGTACGCTGACTTCGCCATCGCTAAAGCGACCTa                                                                                                                                            <  2:2858436/148‑1 (MQ=255)
      gCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAg                                                                                                                                         >  2:1150473/1‑148 (MQ=255)
           acgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAgcggcg                                                                                                                                    >  2:2027238/1‑148 (MQ=255)
                  gCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGa                                                                                                                             >  1:2753072/1‑148 (MQ=255)
                   cATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAg                                                                                                                            >  1:521881/1‑148 (MQ=255)
                    aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGt                                                                                                                           <  2:1271334/148‑1 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAg                                                                                                                       >  2:963259/1‑148 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAg                                                                                                                       >  2:763270/1‑148 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCAAAAGTGGACTGGATGATGAAAATATCACCACCTCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTAACTACAGCGGCGTCGCCGAGTGAAg                                                                                                                       >  2:2423606/1‑148 (MQ=255)
                           aTAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAgtgt                                                                                                                    >  1:1592584/1‑148 (MQ=255)
                             aaCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTAc                                                                                                                  <  1:1416060/148‑1 (MQ=255)
                              aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACa                                                                                                                 >  2:3297726/1‑148 (MQ=255)
                              aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACa                                                                                                                 >  2:3297716/1‑148 (MQ=255)
                              aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACa                                                                                                                 >  1:1971555/1‑148 (MQ=255)
                                 aCTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGc                                                                                                              >  1:609374/1‑148 (MQ=255)
                                  cTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCg                                                                                                             >  1:1693295/1‑148 (MQ=255)
                                  cTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCg                                                                                                             >  2:2600449/1‑148 (MQ=255)
                                   tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCgg                                                                                                            >  2:2762908/1‑148 (MQ=255)
                                   tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCgg                                                                                                            >  2:2762959/1‑148 (MQ=255)
                                        ccATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGc                                                                                                       >  1:3275712/1‑148 (MQ=255)
                                          aTCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCaa                                                                                                     <  2:2086603/148‑1 (MQ=255)
                                           tCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAt                                                                                                    >  1:1488728/1‑148 (MQ=255)
                                           tCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATATGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAt                                                                                                    >  1:2118433/1‑148 (MQ=255)
                                              ggTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACg                                                                                                 <  2:703692/148‑1 (MQ=255)
                                              ggTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACg                                                                                                 >  2:835482/1‑148 (MQ=255)
                                              ggTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTAGGCAATACg                                                                                                 >  1:189543/1‑148 (MQ=255)
                                              ggTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGATACCTACAGCGTCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACg                                                                                                 >  2:2896114/1‑148 (MQ=255)
                                                ttGTCGTTAGTATGGGCACAAGTGGAATGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGtt                                                                                               >  2:618267/1‑148 (MQ=255)
                                                ttGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGtt                                                                                               <  2:254555/148‑1 (MQ=255)
                                                 tGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtg                                                                                              <  2:2936710/148‑1 (MQ=255)
                                                 tGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtg                                                                                              <  2:2936935/148‑1 (MQ=255)
                                                   tCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtgtg                                                                                            >  2:3309205/1‑148 (MQ=255)
                                                   tCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATATGTACGCTGACTTCGCCATCGCTAAAGCTTCCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGGTGGGACTACgtggtg                                                                                            >  2:2591946/1‑148 (MQ=255)
                                                      ttAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTa                                                                                         >  2:1271859/1‑148 (MQ=255)
                                                       tagtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAg                                                                                        >  1:988455/1‑148 (MQ=255)
                                                          tagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTc                                                                                     >  2:3198910/1‑148 (MQ=255)
                                                             gggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCgg                                                                                  >  2:1538467/1‑148 (MQ=255)
                                                               gCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCgggg                                                                                <  1:2763897/148‑1 (MQ=255)
                                                               gCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCgggg                                                                                <  2:2306990/148‑1 (MQ=255)
                                                                cacaAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGt                                                                               >  2:3039056/1‑148 (MQ=255)
                                                                  caAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTgg                                                                             >  1:948016/1‑148 (MQ=255)
                                                                    aGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCTAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCg                                                                           >  1:1176088/1‑148 (MQ=255)
                                                                      tGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGtt                                                                         <  1:1949001/148‑1 (MQ=255)
                                                                       ggACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTa                                                                        <  2:2818241/148‑1 (MQ=255)
                                                                         aCTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTAcc                                                                      <  2:609374/148‑1 (MQ=255)
                                                                         aCTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTAcc                                                                      <  1:2815541/148‑1 (MQ=255)
                                                                          cTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCa                                                                     >  2:454547/1‑148 (MQ=255)
                                                                          cTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCa                                                                     >  2:3250251/1‑148 (MQ=255)
                                                                            ggATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGc                                                                   >  2:2598157/1‑148 (MQ=255)
                                                                             gatgatGAAAATATCACCACCGCGTACATTTTCATTAAATTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCa                                                                  >  1:333333/1‑148 (MQ=255)
                                                                                gatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaa                                                               <  1:2600449/148‑1 (MQ=255)
                                                                                  tgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAg                                                             <  2:189543/148‑1 (MQ=255)
                                                                                  tgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAg                                                             >  1:2397405/1‑148 (MQ=255)
                                                                                   gaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGc                                                            >  1:2190402/1‑148 (MQ=255)
                                                                                   gaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGc                                                            >  1:282606/1‑148 (MQ=255)
                                                                                     aaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCtt                                                          <  2:2595990/148‑1 (MQ=255)
                                                                                     aaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCtt                                                          <  2:504975/148‑1 (MQ=255)
                                                                                     aaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCtt                                                          <  2:502219/148‑1 (MQ=255)
                                                                                        tatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCat                                                       >  2:1170199/1‑148 (MQ=255)
                                                                                        tatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCat                                                       >  1:3277166/1‑148 (MQ=255)
                                                                                        tatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCat                                                       <  1:1613794/148‑1 (MQ=255)
                                                                                         atCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCata                                                      <  1:454547/148‑1 (MQ=255)
                                                                                          tCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCatat                                                     <  1:3075011/148‑1 (MQ=255)
                                                                                          tCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCatat                                                     <  2:396011/148‑1 (MQ=255)
                                                                                           caccacCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATc                                                    <  1:2097370/148‑1 (MQ=255)
                                                                                           caccacCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATc                                                    <  1:432878/148‑1 (MQ=255)
                                                                                                  gcgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACg                                             <  1:3121633/148‑1 (MQ=255)
                                                                                                   cgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACga                                            >  1:2506122/1‑148 (MQ=255)
                                                                                                    gTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACgag                                           <  2:1379761/148‑1 (MQ=255)
                                                                                                     tACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACgaga                                          >  2:533739/1‑148 (MQ=255)
                                                                                                     tACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACgaga                                          >  2:2985794/1‑148 (MQ=255)
                                                                                                     tACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACgaga                                          >  2:1781030/1‑148 (MQ=255)
                                                                                                     tACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGAGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACgaga                                          <  2:1658348/148‑1 (MQ=255)
                                                                                                      aCATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaa                                         >  1:183847/1‑148 (MQ=255)
                                                                                                        aTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaaga                                       <  1:1488189/148‑1 (MQ=255)
                                                                                                        aTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaaga                                       <  1:3284274/148‑1 (MQ=255)
                                                                                                             cATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTc                                  >  2:22107/1‑148 (MQ=255)
                                                                                                             cATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTc                                  <  2:1592584/148‑1 (MQ=255)
                                                                                                              aTTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCa                                 <  1:2976774/148‑1 (MQ=255)
                                                                                                               ttAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAg                                >  2:532241/1‑148 (MQ=255)
                                                                                                                   tttGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCTGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCAt                            <  1:1615998/148‑1 (MQ=255)
                                                                                                                   tttGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCAt                            >  1:2419974/1‑148 (MQ=255)
                                                                                                                   tttGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCAt                            >  1:2645449/1‑148 (MQ=255)
                                                                                                                   tttGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCAt                            >  2:578713/1‑148 (MQ=255)
                                                                                                                        aCGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTc                       >  1:1103159/1‑148 (MQ=255)
                                                                                                                        aCGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTc                       >  1:2056116/1‑148 (MQ=255)
                                                                                                                        aCGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTc                       <  1:1810215/148‑1 (MQ=255)
                                                                                                                           cTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCAt                    <  1:980576/148‑1 (MQ=255)
                                                                                                                           cTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCAt                    <  1:2150963/148‑1 (MQ=255)
                                                                                                                            tGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCAtt                   >  2:2175868/1‑148 (MQ=255)
                                                                                                                              aCTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtgg                 <  2:739875/148‑1 (MQ=255)
                                                                                                                               cTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtggt                <  2:2657177/148‑1 (MQ=255)
                                                                                                                               cTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtggt                <  1:2692159/148‑1 (MQ=255)
                                                                                                                                ttCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtggtg               >  2:2366996/1‑148 (MQ=255)
                                                                                                                                ttCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtggtg               >  1:754766/1‑148 (MQ=255)
                                                                                                                                       tCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGaa        <  1:3198910/148‑1 (MQ=255)
                                                                                                                                       tCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGaa        <  1:22107/148‑1 (MQ=255)
                                                                                                                                         gCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATc      <  2:643130/148‑1 (MQ=255)
                                                                                                                                           tAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTg    <  1:2175302/148‑1 (MQ=255)
                                                                                                                                             aaGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGcc  >  1:682369/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCC  >  CP009273/1257087‑1257375

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: