Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,187,518 G→T 13.4% D230Y (GAT→TAT)  nudC → NADH pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,187,5180GT13.4% 955.6 / 96.8 291D230Y (GAT→TAT) nudCNADH pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (158/94);  new base T (18/21);  total (176/115)
Fisher's exact test for biased strand distribution p-value = 5.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 298 total aligned reads displayed.

GAATTAAAGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGT  >  CP009273/4187369‑4187667
                                                                                                                                                     |                                                                                                                                                     
gAATTAAAGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACg                                                                                                                                                       <  2:4582810/150‑1 (MQ=255)
  aTTAAAGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGAt                                                                                                                                                     <  1:3753516/150‑1 (MQ=255)
   ttAAAGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTCTGACGAtt                                                                                                                                                    >  2:616613/1‑150 (MQ=255)
   ttAAAGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGAtt                                                                                                                                                    <  1:2334016/150‑1 (MQ=255)
    tAAAGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGAttt                                                                                                                                                   >  2:1768265/1‑150 (MQ=255)
    tAAAGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGAttt                                                                                                                                                   >  1:3499563/1‑150 (MQ=255)
     aaaGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACTATTTg                                                                                                                                                  <  2:2477184/150‑1 (MQ=255)
       aGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACTATTTGcc                                                                                                                                                >  1:4117683/1‑150 (MQ=255)
       aGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGcc                                                                                                                                                <  1:2733302/150‑1 (MQ=255)
        gTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCg                                                                                                                                               <  2:543157/150‑1 (MQ=255)
           aaaaaCTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTa                                                                                                                                            >  1:3546115/1‑150 (MQ=255)
           aaaaaCTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTa                                                                                                                                            >  1:4438242/1‑150 (MQ=255)
           aaaaaCTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTa                                                                                                                                            <  2:196726/150‑1 (MQ=255)
            aaaaCTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTAc                                                                                                                                           >  2:3609074/1‑150 (MQ=255)
              aaCTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTc                                                                                                                                         >  1:569041/1‑150 (MQ=255)
              aaCTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTACTCAGTCTTTAATGACCGCGTTTATGGCGGAATATTACAGCGGCGACATCGTGATAGACCCGAAAGAATTTCTCGAGGCGAACTGGTATCGCTATTACGATTTGCCGTTACTc                                                                                                                                         >  1:3176263/1‑150 (MQ=255)
                cTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTccc                                                                                                                                       <  2:1232571/150‑1 (MQ=255)
                cTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTccc                                                                                                                                       >  2:1869054/1‑150 (MQ=255)
                    cGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCccgccg                                                                                                                                   <  2:2797429/150‑1 (MQ=255)
                     gTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCcgccgc                                                                                                                                  >  2:745804/1‑150 (MQ=255)
                      ttACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCgccgcc                                                                                                                                 >  2:2359518/1‑150 (MQ=255)
                       tACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATAGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGccc                                                                                                                                >  1:3239207/1‑150 (MQ=255)
                         cGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCgg                                                                                                                              >  1:1422974/1‑150 (MQ=255)
                          gTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGc                                                                                                                             >  2:4289348/1‑150 (MQ=255)
                            gACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCAc                                                                                                                           >  1:4699758/1‑150 (MQ=255)
                             aCTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCAcc                                                                                                                          >  1:74280/1‑150 (MQ=255)
                              cTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCg                                                                                                                         >  1:3670152/1‑150 (MQ=255)
                               ttCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGt                                                                                                                        >  1:2618013/1‑150 (MQ=255)
                                tctcAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTa                                                                                                                       >  1:2586059/1‑150 (MQ=255)
                                  tcAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACTATTTGCCGTTACTCCCGCCGCCCGGCACCGTAgc                                                                                                                     <  1:4397756/150‑1 (MQ=255)
                                  tcAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACTATTTGCCGTTACTCCCGCCGCCCGGCACCGTAgc                                                                                                                     >  2:4233720/1‑150 (MQ=255)
                                   cAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAgcg                                                                                                                    >  1:4144208/1‑150 (MQ=255)
                                        gTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACTATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCg                                                                                                               <  1:719624/150‑1 (MQ=255)
                                         tGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGt                                                                                                              >  1:963390/1‑150 (MQ=255)
                                            ccGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCTAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTAATCCCGCCGCCCGGCACCGTAGCGAGCCGTCTg                                                                                                           >  1:3730621/1‑150 (MQ=255)
                                                ttCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATag                                                                                                       >  2:4447360/1‑150 (MQ=255)
                                                  ccTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATagaa                                                                                                     <  1:2127682/150‑1 (MQ=255)
                                                   cTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATagaag                                                                                                    >  1:3086951/1‑150 (MQ=255)
                                                   cTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATagaag                                                                                                    >  1:717859/1‑150 (MQ=255)
                                                      aGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATa                                                                                                 >  1:3865470/1‑150 (MQ=255)
                                                      aGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATa                                                                                                 <  1:847019/150‑1 (MQ=255)
                                                       gTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATAc                                                                                                <  2:4724726/150‑1 (MQ=255)
                                                         cTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACgg                                                                                              <  2:303828/150‑1 (MQ=255)
                                                           ttAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTg                                                                                            <  1:1108483/150‑1 (MQ=255)
                                                            tAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACTATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTgg                                                                                           >  1:2768127/1‑150 (MQ=255)
                                                            tAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACTATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTgg                                                                                           <  2:4068254/150‑1 (MQ=255)
                                                            tAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTgg                                                                                           <  2:3821346/150‑1 (MQ=255)
                                                            tAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTgg                                                                                           >  2:3592571/1‑150 (MQ=255)
                                                             aaTGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGc                                                                                          <  2:909685/150‑1 (MQ=255)
                                                               tGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGa                                                                                        >  2:1519511/1‑150 (MQ=255)
                                                                gACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGAt                                                                                       >  1:2824118/1‑150 (MQ=255)
                                                                 aCCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACTATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGAtg                                                                                      >  2:2365081/1‑150 (MQ=255)
                                                                     cgTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTc                                                                                  >  1:1841842/1‑150 (MQ=255)
                                                                      gTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCg                                                                                 >  2:264197/1‑150 (MQ=255)
                                                                       tttATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCgg                                                                                >  2:401411/1‑150 (MQ=255)
                                                                        ttATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCggg                                                                               <  1:3752089/150‑1 (MQ=255)
                                                                            ggCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCaga                                                                           >  2:1690937/1‑150 (MQ=255)
                                                                                 aaTATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATg                                                                      >  2:2747775/1‑150 (MQ=255)
                                                                                 aaTATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATg                                                                      >  1:126240/1‑150 (MQ=255)
                                                                                 aaTATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATg                                                                      <  1:4021615/150‑1 (MQ=255)
                                                                                 aaTATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATg                                                                      >  2:374049/1‑150 (MQ=255)
                                                                                   tatGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAg                                                                    <  1:1749274/150‑1 (MQ=255)
                                                                                   tatGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAg                                                                    <  2:3391563/150‑1 (MQ=255)
                                                                                      gACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGtga                                                                 >  1:1296511/1‑150 (MQ=255)
                                                                                           cggcgACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACTATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATa                                                            <  1:330594/150‑1 (MQ=255)
                                                                                               gACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACt                                                        >  2:3635476/1‑150 (MQ=255)
                                                                                                  aTCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGAc                                                     <  2:1259658/150‑1 (MQ=255)
                                                                                                  aTCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGAc                                                     <  1:281394/150‑1 (MQ=255)
                                                                                                       gATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCt                                                >  2:1359377/1‑150 (MQ=255)
                                                                                                       gATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGAGGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCt                                                >  2:4392842/1‑150 (MQ=255)
                                                                                                        aTCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTg                                               >  1:1046211/1‑150 (MQ=255)
                                                                                                        aTCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTg                                               >  1:1045991/1‑150 (MQ=255)
                                                                                                              ccGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCAc                                         <  2:3232339/150‑1 (MQ=255)
                                                                                                               cGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACt                                        <  1:4599776/150‑1 (MQ=255)
                                                                                                                 aaaGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTaa                                      <  1:4572259/150‑1 (MQ=255)
                                                                                                                 aaaGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTaa                                      >  1:318986/1‑150 (MQ=255)
                                                                                                                    gAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGa                                   <  1:1716522/150‑1 (MQ=255)
                                                                                                                       ttGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACa                                >  1:2671986/1‑150 (MQ=255)
                                                                                                                       ttGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACa                                <  1:1768265/150‑1 (MQ=255)
                                                                                                                          cTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGcc                             <  1:2616801/150‑1 (MQ=255)
                                                                                                                           tCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCa                            >  2:1435972/1‑150 (MQ=255)
                                                                                                                            cGAGGCGAACTGGTATCGCTATGACTATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCaa                           >  2:1553355/1‑150 (MQ=255)
                                                                                                                            cGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCaa                           >  1:1907257/1‑150 (MQ=255)
                                                                                                                               ggCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAAt                        <  2:4587125/150‑1 (MQ=255)
                                                                                                                                gCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATg                       <  2:1517163/150‑1 (MQ=255)
                                                                                                                                gCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATg                       <  1:420792/150‑1 (MQ=255)
                                                                                                                                  gAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTCAGGAACAGCCAAAATGAc                     >  1:446641/1‑150 (MQ=255)
                                                                                                                                    aCTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCg                   >  2:211245/1‑150 (MQ=255)
                                                                                                                                       ggTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAAc                >  1:1538962/1‑150 (MQ=255)
                                                                                                                                        gTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACt               <  2:196138/150‑1 (MQ=255)
                                                                                                                                         tATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGAGGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACtt              >  1:242418/1‑150 (MQ=255)
                                                                                                                                           tCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTaa            <  2:668267/150‑1 (MQ=255)
                                                                                                                                            cGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTaaa           <  2:3350933/150‑1 (MQ=255)
                                                                                                                                               tATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAAc        <  1:4308599/150‑1 (MQ=255)
                                                                                                                                                aTGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACg       >  2:3161368/1‑150 (MQ=255)
                                                                                                                                                aTGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACg       >  2:799181/1‑150 (MQ=255)
                                                                                                                                                  gACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGACCAGCCAAAATGACCGAACTTAAAAACGAt     >  2:4610685/1‑150 (MQ=255)
                                                                                                                                                    cGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCg   >  1:2966030/1‑150 (MQ=255)
                                                                                                                                                     gATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGt  >  1:1974194/1‑150 (MQ=255)
                                                                                                                                                     gATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGt  >  1:1376662/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAATTAAAGTTAAAAACTTGCGTTACGTGACTTCTCAGCCGTGGCCGTTTCCTCAGTCTTTAATGACCGCGTTTATGGCGGAATATGACAGCGGCGACATCGTGATCGACCCGAAAGAATTGCTCGAGGCGAACTGGTATCGCTATGACGATTTGCCGTTACTCCCGCCGCCCGGCACCGTAGCGCGCCGTCTGATAGAAGATACGGTGGCGATGTGTCGGGCAGAGTATGAGTGATGATACACTGACCGCCTGACGCACTAAGGAACAGCCAAAATGACCGAACTTAAAAACGATCGT  >  CP009273/4187369‑4187667

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