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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,257,2320GT99.0% 709.2 / ‑1.7 222R34S (CGC→AGC) prsphosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (126/93);  total (127/95)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 225 total aligned reads displayed.

GCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCG  >  CP009273/1257087‑1257376
                                                                                                                                                 |                                                                                                                                                
gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTAc                                                                                                                                                <  2:891024/148‑1 (MQ=255)
gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTAc                                                                                                                                                <  2:2591927/148‑1 (MQ=255)
   gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTa                                                                                                                                             >  2:599174/1‑148 (MQ=255)
   gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTa                                                                                                                                             >  1:3109472/1‑148 (MQ=255)
   gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTa                                                                                                                                             >  2:1558109/1‑148 (MQ=255)
    aaGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGctacctac                                                                                                                                            >  1:3135756/1‑148 (MQ=255)
    aaGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGctacctac                                                                                                                                            >  2:2707126/1‑148 (MQ=255)
        acgacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAgcg                                                                                                                                        >  1:117590/1‑148 (MQ=255)
         cgacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAgcgg                                                                                                                                       >  1:1303932/1‑148 (MQ=255)
                 ggCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCg                                                                                                                               >  1:2228474/1‑148 (MQ=255)
                      cAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGa                                                                                                                          >  1:150521/1‑148 (MQ=255)
                      cAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATATTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGa                                                                                                                          >  2:1210463/1‑148 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGacg                                                                                                                        >  2:3738047/1‑146 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAg                                                                                                                        >  2:1847024/1‑148 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAg                                                                                                                        <  2:2128193/148‑1 (MQ=255)
                          cATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAgtg                                                                                                                      >  1:2772404/1‑148 (MQ=255)
                           aTAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAgtgt                                                                                                                     <  2:560551/148‑1 (MQ=255)
                             aaCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTAc                                                                                                                   >  2:1246194/1‑148 (MQ=255)
                              aCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACa                                                                                                                  <  1:1210646/148‑1 (MQ=255)
                               cGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAg                                                                                                                 <  2:3332052/148‑1 (MQ=255)
                               cGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAg                                                                                                                 >  2:590339/1‑148 (MQ=255)
                               cGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAg                                                                                                                 >  1:1242709/1‑148 (MQ=255)
                               cGACTAAGTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAg                                                                                                                 <  2:3331969/148‑1 (MQ=255)
                                 aCTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGc                                                                                                               >  2:3479884/1‑148 (MQ=255)
                                   tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCgg                                                                                                             <  1:107028/148‑1 (MQ=255)
                                   tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCgg                                                                                                             <  1:106992/148‑1 (MQ=255)
                                       tCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTgg                                                                                                         <  2:2907527/148‑1 (MQ=255)
                                       tCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTgg                                                                                                         >  2:2366728/1‑148 (MQ=255)
                                       tCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTgg                                                                                                         >  1:298412/1‑148 (MQ=255)
                                       tCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTgg                                                                                                         >  2:2366683/1‑148 (MQ=255)
                                       tCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTgg                                                                                                         >  2:1358623/1‑148 (MQ=255)
                                        ccATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGc                                                                                                        <  1:2566368/148‑1 (MQ=255)
                                         cATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCa                                                                                                       <  1:448453/148‑1 (MQ=255)
                                         cATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCa                                                                                                       <  1:447971/148‑1 (MQ=255)
                                           tCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAt                                                                                                     <  1:2463652/148‑1 (MQ=255)
                                           tCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAt                                                                                                     >  1:3511397/1‑148 (MQ=255)
                                            cAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATa                                                                                                    <  1:3324744/148‑1 (MQ=255)
                                                   tCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtgtg                                                                                             >  2:3450941/1‑148 (MQ=255)
                                                   tCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtgtg                                                                                             <  2:3109165/148‑1 (MQ=255)
                                                   tCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtgtg                                                                                             >  2:1752347/1‑148 (MQ=255)
                                                     gTTTGTAGGGGCACAAGTGGACTGGATGATGAATATATCACCTCCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACATTCAGCGGCGTCGCCGAGTGTAGTGTACAGGCGGTTGGCAATACGTTGTGCt                                                                                           <  2:459122/148‑1 (MQ=255)
                                                      ttAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTa                                                                                          >  1:1785790/1‑148 (MQ=255)
                                                      ttAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTa                                                                                          >  1:1786751/1‑148 (MQ=255)
                                                         gtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGtt                                                                                       <  2:1721262/148‑1 (MQ=255)
                                                          tagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTc                                                                                      <  1:965407/148‑1 (MQ=255)
                                                          tagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTc                                                                                      >  2:2474804/1‑148 (MQ=255)
                                                            ggggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCg                                                                                    >  1:1271983/1‑148 (MQ=255)
                                                               gCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCgggg                                                                                 <  2:1958949/148‑1 (MQ=255)
                                                                 acaAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGGTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTg                                                                               >  1:2830017/1‑148 (MQ=255)
                                                                    aGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCg                                                                            <  1:724444/148‑1 (MQ=255)
                                                                    aGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCg                                                                            >  1:721174/1‑148 (MQ=255)
                                                                     gTGGACTGGATGATGAAAATAACACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGt                                                                           >  1:2477299/1‑148 (MQ=255)
                                                                          cTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCa                                                                      >  2:432767/1‑148 (MQ=255)
                                                                                gatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaa                                                                >  2:813631/1‑148 (MQ=255)
                                                                                 atgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaaa                                                               >  2:2467320/1‑148 (MQ=255)
                                                                                 atgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaaa                                                               >  2:2216825/1‑148 (MQ=255)
                                                                                 atgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaaa                                                               >  2:2619901/1‑148 (MQ=255)
                                                                                 atgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaaa                                                               <  2:1998129/148‑1 (MQ=255)
                                                                                  tgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAg                                                              >  2:1013433/1‑148 (MQ=255)
                                                                                   gaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGc                                                             <  1:370887/148‑1 (MQ=255)
                                                                                    aaaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCt                                                            >  1:2531993/1‑148 (MQ=255)
                                                                                     aaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCtt                                                           <  2:2379979/148‑1 (MQ=255)
                                                                                      aaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGTGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                          <  2:3562216/148‑1 (MQ=255)
                                                                                      aaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                          >  2:1784055/1‑148 (MQ=255)
                                                                                      aaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                          <  2:588261/148‑1 (MQ=255)
                                                                                      aaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                          <  1:3319069/148‑1 (MQ=255)
                                                                                      aaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                          >  2:842251/1‑148 (MQ=255)
                                                                                      aaTATCACCACCGCGTACATTTTCATTAATTTGTAAGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                          >  1:41429/1‑148 (MQ=255)
                                                                                       atatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCa                                                         <  1:3370229/148‑1 (MQ=255)
                                                                                       atatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCAACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGGTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCa                                                         >  1:2264938/1‑148 (MQ=255)
                                                                                        tatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACTTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCat                                                        >  2:1539588/1‑148 (MQ=255)
                                                                                        tatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCat                                                        <  2:1123526/148‑1 (MQ=255)
                                                                                          tCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCatat                                                      >  2:2202987/1‑148 (MQ=255)
                                                                                          tCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCatat                                                      >  1:3148477/1‑148 (MQ=255)
                                                                                               accGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGc                                                 >  2:2727632/1‑148 (MQ=255)
                                                                                                  gcgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACg                                              >  2:3632134/1‑148 (MQ=255)
                                                                                                     tACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACgaga                                           >  1:2824846/1‑148 (MQ=255)
                                                                                                      aCATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCAGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaa                                          >  2:1793322/1‑148 (MQ=255)
                                                                                                        aTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaaga                                        <  2:117590/148‑1 (MQ=255)
                                                                                                         ttttCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAgaagaa                                       >  1:1103148/1‑148 (MQ=255)
                                                                                                         ttttCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAgaagaa                                       >  2:2999733/1‑148 (MQ=255)
                                                                                                          tttCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAAc                                      >  2:3215348/1‑148 (MQ=255)
                                                                                                            tCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCt                                    <  1:1845669/148‑1 (MQ=255)
                                                                                                              aTTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCa                                  <  2:2835125/148‑1 (MQ=255)
                                                                                                                tAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAgg                                <  1:2692177/148‑1 (MQ=255)
                                                                                                                tAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAgg                                <  1:3089690/148‑1 (MQ=255)
                                                                                                                  aTTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCa                              >  2:1414886/1‑148 (MQ=255)
                                                                                                                   tttGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCAt                             <  1:1359020/148‑1 (MQ=255)
                                                                                                                   tttGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCAt                             <  1:2366683/148‑1 (MQ=255)
                                                                                                                           cTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCAt                     >  1:1333994/1‑148 (MQ=255)
                                                                                                                             gACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtg                   <  2:2772404/148‑1 (MQ=255)
                                                                                                                                    ccATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGaaa            >  1:2419396/1‑148 (MQ=255)
                                                                                                                                    ccATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGaaa            >  1:2420160/1‑148 (MQ=255)
                                                                                                                                     cATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAg           <  1:2215451/148‑1 (MQ=255)
                                                                                                                                       tCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGaa         >  1:1478021/1‑148 (MQ=255)
                                                                                                                                       tCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGaa         >  1:1697640/1‑148 (MQ=255)
                                                                                                                                        cGCTACAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAAt        <  2:2485381/148‑1 (MQ=255)
                                                                                                                                            aaaGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGc    >  2:828503/1‑148 (MQ=255)
                                                                                                                                            aaaGCTACCTACAGCGGCGTCGACGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCATGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGc    >  2:912640/1‑148 (MQ=255)
                                                                                                                                              aGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCg  >  1:47704/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCG  >  CP009273/1257087‑1257376

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: ‑