Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,881,078 G→T 71.0% G198G (GGC→GGA cysH ← phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,881,0780GT71.0% 334.6 / 281.8 296G198G (GGC→GGAcysHphosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent
Reads supporting (aligned to +/- strand):  ref base G (37/49);  new base T (126/84);  total (163/133)
Fisher's exact test for biased strand distribution p-value = 9.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 296 total aligned reads displayed.

CAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTGC  >  CP009273/2880931‑2881225
                                                                                                                                                   |                                                                                                                                                   
 aaTATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCa                                                                                                                                                  >  2:2974382/1‑150 (MQ=255)
  atatTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCAt                                                                                                                                                 <  2:4194726/150‑1 (MQ=255)
    atTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGt                                                                                                                                               >  1:3797949/1‑150 (MQ=255)
    atTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGt                                                                                                                                               >  1:1446321/1‑150 (MQ=255)
      ttACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAATTATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGttt                                                                                                                                             >  2:3525903/1‑150 (MQ=255)
        aCCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGttttt                                                                                                                                           <  2:2811767/150‑1 (MQ=255)
           cgtCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCa                                                                                                                                        <  2:2311531/148‑1 (MQ=255)
           cTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCa                                                                                                                                        >  1:1045955/1‑150 (MQ=255)
           cTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCa                                                                                                                                        >  1:2238571/1‑150 (MQ=255)
                   aaCCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGa                                                                                                                                >  2:3334453/1‑150 (MQ=255)
                    aCCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGAt                                                                                                                               <  1:3116612/150‑1 (MQ=255)
                      ccACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATaa                                                                                                                             >  1:401964/1‑150 (MQ=255)
                        acaTTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAAt                                                                                                                           <  2:3877803/150‑1 (MQ=255)
                        acaTTCCCTTTTTAAGCAAAAGAAACGCGTTTCTTCTACCGCCATGACTGGTTACCATTTACGGGTGGTTTGTGTATCGACCACCTATAAATATCCTTCATCCCATAATGTGTGTTATTTCATTACATGTTTTTGCAGGTACTGATAAAt                                                                                                                           <  1:2495482/150‑1 (MQ=255)
                          aTTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAg                                                                                                                         >  1:3601616/1‑150 (MQ=255)
                          aTTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAg                                                                                                                         <  2:4864127/150‑1 (MQ=255)
                           ttCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGt                                                                                                                        <  1:695802/150‑1 (MQ=255)
                            tCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGtt                                                                                                                       <  2:3789190/150‑1 (MQ=255)
                            tCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGtt                                                                                                                       >  2:4006680/1‑150 (MQ=255)
                             cccTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTc                                                                                                                      <  1:3657879/150‑1 (MQ=255)
                             cccTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTc                                                                                                                      >  1:1098092/1‑150 (MQ=255)
                             ccacttttAATACAAATAAACGCGTTTCTTCTTCCGCCTTGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCATTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGTTACTGATAAATAGTTc                                                                                                                      <  1:736420/146‑1 (MQ=255)
                                tttttAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGt                                                                                                                   <  2:687101/150‑1 (MQ=255)
                                 ttttAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGtt                                                                                                                  <  2:5208859/150‑1 (MQ=255)
                                   ttAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTAt                                                                                                                <  1:1691159/150‑1 (MQ=255)
                                    tAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGTTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATc                                                                                                               >  2:1435159/1‑150 (MQ=255)
                                     aaGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATcc                                                                                                              <  1:2342681/150‑1 (MQ=255)
                                        ccAAAGAAACGCGAACCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGGTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGAGAAATAGTTGGGTTATCCCAg                                                                                                           >  2:3252736/1‑150 (MQ=255)
                                         cAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGt                                                                                                          >  2:2238244/1‑150 (MQ=255)
                                          aaaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCTAGTc                                                                                                         <  2:3819070/150‑1 (MQ=255)
                                          aaaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTc                                                                                                         <  2:405674/150‑1 (MQ=255)
                                          aaaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTc                                                                                                         <  1:39694/150‑1 (MQ=255)
                                           aaGAAACGCGTTTCGTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCg                                                                                                        <  1:4990451/150‑1 (MQ=255)
                                            aGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGa                                                                                                       >  1:4118825/1‑150 (MQ=255)
                                             gAAACGCGTTTCTTCTTCCGCCATGCCGGGTACCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGTGAGATATATCAGGCCATGATTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGAt                                                                                                      >  2:172207/1‑150 (MQ=255)
                                             gAAACGCGTATCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGAAACGCCCACCGATAAATATACTTCATCCCATAATGGGTGATATTTAATTCCATTTTTTTGAAGGTAATGATAAATAGTTCGGTTATCCCAGTAGAt                                                                                                      >  1:3655380/1‑150 (MQ=255)
                                              aaaCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATa                                                                                                     >  2:543201/1‑150 (MQ=255)
                                               aaCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATaa                                                                                                    <  2:2389148/150‑1 (MQ=255)
                                                  gcgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCg                                                                                                 >  2:4961004/1‑150 (MQ=255)
                                                   cgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCgg                                                                                                >  1:888327/1‑150 (MQ=255)
                                                   cgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCgg                                                                                                >  1:2947064/1‑150 (MQ=255)
                                                   cgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCgg                                                                                                <  2:4491292/150‑1 (MQ=255)
                                                   cgTTTCATATTCCGCCATGACGGGTTCCCATTTACTGGTTGTATGGGTATCGCCCACCGATAAATAACCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTATTGCAGGTACTGATAAATAGTTCTGTAATACCAGTCTATAATCgg                                                                                                >  1:4718665/1‑150 (MQ=255)
                                                    gTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgc                                                                                               >  1:424024/1‑150 (MQ=255)
                                                    gTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgc                                                                                               <  2:5089351/150‑1 (MQ=255)
                                                         ttcttcCGCCATGCCGGGGTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgcagccc                                                                                          <  2:4631369/150‑3 (MQ=255)
                                                            ttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACttt                                                                                       <  1:1293349/150‑1 (MQ=255)
                                                              ccGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTaa                                                                                     >  1:613684/1‑150 (MQ=255)
                                                               cGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTaaa                                                                                    <  1:2147402/150‑1 (MQ=255)
                                                                gCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAAt                                                                                   <  1:331666/150‑1 (MQ=255)
                                                                  cATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATAc                                                                                 >  1:2206863/1‑150 (MQ=255)
                                                                   aTGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACg                                                                                >  1:4542246/1‑150 (MQ=255)
                                                                      ccGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCa                                                                             <  1:2108810/150‑1 (MQ=255)
                                                                      ccGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCa                                                                             <  1:4040340/150‑1 (MQ=255)
                                                                       cGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCAc                                                                            >  2:4116439/1‑150 (MQ=255)
                                                                       cGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCAc                                                                            >  2:2666877/1‑150 (MQ=255)
                                                                          gTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCt                                                                         <  1:3251016/150‑1 (MQ=255)
                                                                          gTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAAACGGCAGCACTTTAAATACGCCACGCt                                                                         >  1:99526/1‑150 (MQ=255)
                                                                           ttCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTg                                                                        <  1:2621703/150‑1 (MQ=255)
                                                                              ccATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAAt                                                                     >  2:2501/1‑150 (MQ=255)
                                                                              ccATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAAt                                                                     >  1:5160931/1‑150 (MQ=255)
                                                                               cATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAAtt                                                                    >  2:2663168/1‑150 (MQ=255)
                                                                               cATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAAtt                                                                    <  1:192210/150‑1 (MQ=255)
                                                                                  ttACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGcc                                                                 <  1:851317/150‑1 (MQ=255)
                                                                                   tACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCa                                                                <  2:1075092/150‑1 (MQ=255)
                                                                                     cGGGTTGTATGGGTATCGCCCACAGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGc                                                              <  2:4756740/150‑1 (MQ=255)
                                                                                          tGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCgg                                                         >  2:2773786/1‑150 (MQ=255)
                                                                                           gTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGt                                                        >  1:791793/1‑150 (MQ=255)
                                                                                            tATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTa                                                       >  2:1927618/1‑150 (MQ=255)
                                                                                               gggTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAAt                                                    <  2:130852/150‑1 (MQ=255)
                                                                                                   aTCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGc                                                >  2:4306026/1‑150 (MQ=255)
                                                                                                       cccACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGa                                            >  2:2439246/1‑150 (MQ=255)
                                                                                                        ccACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGAc                                           >  2:1902178/1‑150 (MQ=255)
                                                                                                        ccACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGAc                                           <  2:1243230/150‑1 (MQ=255)
                                                                                                         cACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACt                                          <  1:3880040/150‑1 (MQ=255)
                                                                                                            cGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGcc                                       <  2:892009/150‑1 (MQ=255)
                                                                                                            cGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGcc                                       <  1:3429645/150‑1 (MQ=255)
                                                                                                              aTAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCgg                                     >  1:4337960/1‑150 (MQ=255)
                                                                                                                aaaTATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAt                                   >  1:1648333/1‑150 (MQ=255)
                                                                                                                aaaTATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAt                                   >  1:4308367/1‑150 (MQ=255)
                                                                                                                 aaTATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAtt                                  >  2:1828588/1‑150 (MQ=255)
                                                                                                                 aaTATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAtt                                  >  1:375431/1‑150 (MQ=255)
                                                                                                                 aaTATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAtt                                  >  1:4051954/1‑150 (MQ=255)
                                                                                                                  atatCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTg                                 >  1:4808449/1‑150 (MQ=255)
                                                                                                                  atatCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTg                                 >  1:1833987/1‑150 (MQ=255)
                                                                                                                    atCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCCGGTACTGATAATGAGTTCGGTTATCCCAGTCGATAAGCGGCAGCACTTTAAATACGCCACGCTGACTTGCCCGCACCGGTACATTGGCTCGACTGCCGGATTGtt                               >  2:679077/1‑150 (MQ=255)
                                                                                                                    atCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGtt                               >  1:1718882/1‑150 (MQ=255)
                                                                                                                    atCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGtt                               <  2:1572033/150‑1 (MQ=255)
                                                                                                                      ccTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCa                             >  2:2628819/1‑150 (MQ=255)
                                                                                                                      ccTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCa                             >  1:2657922/1‑150 (MQ=255)
                                                                                                                           aTCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGgcg                        <  2:791793/150‑1 (MQ=255)
                                                                                                                               cATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAgg                    <  2:1153057/150‑1 (MQ=255)
                                                                                                                                 tAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGcc                  >  1:4389612/1‑150 (MQ=255)
                                                                                                                                  aaTGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCa                 >  1:1070785/1‑150 (MQ=255)
                                                                                                                                  aaTGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCa                 <  2:4437280/150‑1 (MQ=255)
                                                                                                                                    tGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGc               >  2:2167399/1‑150 (MQ=255)
                                                                                                                                       gTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCaaa            >  2:4261814/1‑150 (MQ=255)
                                                                                                                                          atatTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCa         >  1:2616397/1‑150 (MQ=255)
                                                                                                                                            atTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAgg       >  1:4126156/1‑150 (MQ=255)
                                                                                                                                             tttCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGt      >  1:96752/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTGC  >  CP009273/2880931‑2881225

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