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.5% 682.7 / ‑4.3 209R34S (CGC→AGC) prsphosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (107/101);  total (108/101)
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.94e-01
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 211 total aligned reads displayed.

GCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCG  >  CP009273/1257087‑1257376
                                                                                                                                                 |                                                                                                                                                
gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTAc                                                                                                                                                <  1:64162/148‑1 (MQ=255)
gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTAc                                                                                                                                                <  1:163930/148‑1 (MQ=255)
 cGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTAcc                                                                                                                                               >  2:306722/1‑148 (MQ=255)
   gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTa                                                                                                                                             >  2:856410/1‑148 (MQ=255)
   gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTa                                                                                                                                             <  2:2407755/148‑1 (MQ=255)
   gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTa                                                                                                                                             <  2:646724/148‑1 (MQ=255)
   gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTa                                                                                                                                             >  1:260045/1‑148 (MQ=255)
    aaGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGctacctac                                                                                                                                            >  2:1307194/1‑148 (MQ=255)
    aaGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGctacctac                                                                                                                                            >  2:796718/1‑148 (MQ=255)
     aGCACGACGCAGGGCATCAACCATAACGACTACTTCCATCAGGTTTTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACTCTGACTTCGCCATCGCTAAAGCTACCTACa                                                                                                                                           <  2:2069234/148‑1 (MQ=255)
       cACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAgc                                                                                                                                         <  2:3282868/148‑1 (MQ=255)
           acgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAgcggcg                                                                                                                                     <  2:1407782/148‑1 (MQ=255)
            cgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGt                                                                                                                                    <  2:102970/148‑1 (MQ=255)
                    aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGt                                                                                                                            >  1:2431382/1‑148 (MQ=255)
                    aTCAACCAGACCGACTAATTCCATCAGGTGGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGt                                                                                                                            <  2:1049784/148‑1 (MQ=255)
                       aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGaa                                                                                                                         <  2:2994537/148‑1 (MQ=255)
                       aaCCAGAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGaa                                                                                                                         <  2:2994499/148‑1 (MQ=255)
                        aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAg                                                                                                                        >  1:1393859/1‑148 (MQ=255)
                         ccATAACGAATAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAgt                                                                                                                       <  1:58723/148‑1 (MQ=255)
                            tAACGACTAAGTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTa                                                                                                                    <  1:2129690/148‑1 (MQ=255)
                             aaCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTAc                                                                                                                   >  2:2375757/1‑148 (MQ=255)
                               cGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAg                                                                                                                 <  1:2742970/148‑1 (MQ=255)
                               cGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAg                                                                                                                 >  1:2699013/1‑148 (MQ=255)
                                gACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAgg                                                                                                                <  1:1580036/148‑1 (MQ=255)
                                  cTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCg                                                                                                              >  1:3195569/1‑148 (MQ=255)
                                   tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCgg                                                                                                             <  1:1116994/148‑1 (MQ=255)
                                     aTTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGtt                                                                                                           >  2:1791559/1‑148 (MQ=255)
                                      ttCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTg                                                                                                          >  1:3203539/1‑148 (MQ=255)
                                       tCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTgg                                                                                                         <  1:2550278/148‑1 (MQ=255)
                                       tCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTgg                                                                                                         >  2:2902970/1‑148 (MQ=255)
                                       tCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACGTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTgg                                                                                                         <  2:1407087/148‑1 (MQ=255)
                                         cATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCa                                                                                                       >  1:489393/1‑148 (MQ=255)
                                           tCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATNTGTACGCGGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAt                                                                                                     >  1:1785993/1‑148 (MQ=255)
                                            cAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATa                                                                                                    <  1:3081513/148‑1 (MQ=255)
                                            cAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATa                                                                                                    <  1:2253068/148‑1 (MQ=255)
                                            cAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATa                                                                                                    >  2:3417164/1‑148 (MQ=255)
                                               gTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGt                                                                                                 <  2:2321586/148‑1 (MQ=255)
                                                 tGTCTTTATTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCTCTAAAGCTACCTACAGCTGCGTCGCCGAGTGAAGTTTACAGGCGGTTGGCAATACGttt                                                                                               <  1:2197421/148‑2 (MQ=255)
                                                 tGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtg                                                                                               <  1:3240090/148‑1 (MQ=255)
                                                 tGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtg                                                                                               >  1:2105978/1‑148 (MQ=255)
                                                    cGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGc                                                                                            <  2:1194367/148‑1 (MQ=255)
                                                      ttAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTa                                                                                          <  1:2843274/148‑1 (MQ=255)
                                                       tagtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAg                                                                                         >  2:2791016/1‑148 (MQ=255)
                                                       tagtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAg                                                                                         >  2:3364368/1‑148 (MQ=255)
                                                          tagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTc                                                                                      <  2:1060722/148‑1 (MQ=255)
                                                          tagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTc                                                                                      >  2:547542/1‑148 (MQ=255)
                                                          tagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATGTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTc                                                                                      <  2:3300380/148‑1 (MQ=255)
                                                             ggtcacaAGTGGACTGGAAGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCgg                                                                                   >  2:3340407/4‑148 (MQ=255)
                                                                   aaGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGc                                                                             >  1:376582/1‑148 (MQ=255)
                                                                     gTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGt                                                                           <  2:1428251/148‑1 (MQ=255)
                                                                     gTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGt                                                                           >  2:2048515/1‑148 (MQ=255)
                                                                     gTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGt                                                                           <  2:1944144/148‑1 (MQ=255)
                                                                     gTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGt                                                                           >  1:5239/1‑148 (MQ=255)
                                                                        gACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTAc                                                                        <  1:1719313/148‑1 (MQ=255)
                                                                           tGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAg                                                                     <  1:2954787/148‑1 (MQ=255)
                                                                              atgatgAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaa                                                                  >  2:1670927/1‑148 (MQ=255)
                                                                              atgatgAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaa                                                                  >  2:1821733/1‑148 (MQ=255)
                                                                               tgatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaa                                                                 >  1:3027979/1‑148 (MQ=255)
                                                                                gatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaa                                                                <  1:904289/148‑1 (MQ=255)
                                                                                gatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaa                                                                <  1:3417164/148‑1 (MQ=255)
                                                                                  tgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAg                                                              >  2:1712882/1‑148 (MQ=255)
                                                                                   gaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGc                                                             >  1:1483325/1‑148 (MQ=255)
                                                                                    aaaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCt                                                            <  1:984690/148‑1 (MQ=255)
                                                                                    aaaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCt                                                            <  1:1302743/148‑1 (MQ=255)
                                                                                    aaaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCt                                                            <  1:2169121/148‑1 (MQ=255)
                                                                                      aaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                          <  2:1736395/148‑1 (MQ=255)
                                                                                           caccacCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATc                                                     <  1:874806/148‑1 (MQ=255)
                                                                                              caccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAgg                                                  <  1:1873071/148‑1 (MQ=255)
                                                                                                ccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCa                                                <  1:699691/148‑1 (MQ=255)
                                                                                                  gcgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTTCAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCTGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACg                                              <  1:1051018/148‑1 (MQ=255)
                                                                                                  gcgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACg                                              <  2:1918596/148‑1 (MQ=255)
                                                                                                   cgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACga                                             >  1:1808492/1‑148 (MQ=255)
                                                                                                      aCATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaa                                          >  2:3461792/1‑148 (MQ=255)
                                                                                                      aCATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaa                                          >  2:3459882/1‑148 (MQ=255)
                                                                                                      aCATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaa                                          >  1:2979422/1‑148 (MQ=255)
                                                                                                       cATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaag                                         <  1:1834742/148‑1 (MQ=255)
                                                                                                       cATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaag                                         <  1:3308932/148‑1 (MQ=255)
                                                                                                        aTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaaga                                        <  2:3191009/148‑1 (MQ=255)
                                                                                                          tttCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAAc                                      >  1:682446/1‑148 (MQ=255)
                                                                                                               ttAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAg                                 >  2:3163571/1‑148 (MQ=255)
                                                                                                                tAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAgg                                >  2:2316149/1‑148 (MQ=255)
                                                                                                                 aaTTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGc                               <  1:2506317/148‑1 (MQ=255)
                                                                                                                   tttGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCAt                             >  2:601701/1‑148 (MQ=255)
                                                                                                                      gTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGAATGCg                          <  2:2616524/148‑1 (MQ=255)
                                                                                                                          gCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCa                      <  2:1467669/148‑1 (MQ=255)
                                                                                                                            tGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCAtt                    <  1:952181/148‑1 (MQ=255)
                                                                                                                              aCTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtgg                  >  1:3356742/1‑148 (MQ=255)
                                                                                                                              aCTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtgg                  <  2:2165019/148‑1 (MQ=255)
                                                                                                                              aCTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtgg                  >  2:1363158/1‑148 (MQ=255)
                                                                                                                               cTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtggt                 <  2:1266638/148‑1 (MQ=255)
                                                                                                                                ttCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtggtg                >  1:1669074/1‑148 (MQ=255)
                                                                                                                                     cATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAg           >  1:880383/1‑148 (MQ=255)
                                                                                                                                      aTCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGa          >  2:3016314/1‑148 (MQ=255)
                                                                                                                                      aTCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGa          <  2:2881129/148‑1 (MQ=255)
                                                                                                                                      aTCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGa          <  2:2275585/148‑1 (MQ=255)
                                                                                                                                       tCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGaa         <  1:1732313/148‑1 (MQ=255)
                                                                                                                                         gCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATc       <  1:104102/148‑1 (MQ=255)
                                                                                                                                            aaaGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGc    >  2:412082/1‑148 (MQ=255)
                                                                                                                                             aaGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGcc   >  1:1189172/1‑148 (MQ=255)
                                                                                                                                              aGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCg  >  1:1604435/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: