Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,330,955 G→T 37.0% D290Y (GAT→TAT)  acnA → aconitate hydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,330,9550GT37.0% 469.4 / inf 324D290Y (GAT→TAT) acnAaconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (90/114);  new base T (68/52);  total (158/166)
Fisher's exact test for biased strand distribution p-value = 3.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

CCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTT  >  CP009273/1330808‑1331103
                                                                                                                                                   |                                                                                                                                                    
ccGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTAt                                                                                                                                                    >  1:4537881/1‑150 (MQ=255)
ccGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGAt                                                                                                                                                    <  1:1216574/150‑1 (MQ=255)
 cGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATg                                                                                                                                                   >  1:3561305/1‑150 (MQ=255)
 cGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATg                                                                                                                                                   >  1:2511059/1‑150 (MQ=255)
      tCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATGTTATGGTGATGGTCTg                                                                                                                                              >  2:4433668/1‑150 (MQ=255)
       ccATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTgg                                                                                                                                             <  2:44222/150‑1 (MQ=255)
         aTGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGAt                                                                                                                                           >  2:3493698/1‑150 (MQ=255)
         aTGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGAt                                                                                                                                           <  1:1880200/150‑1 (MQ=255)
         aTGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGAt                                                                                                                                           >  2:1696850/1‑150 (MQ=255)
           gCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCATCTGCTGCGCTAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTc                                                                                                                                         <  2:3071622/150‑1 (MQ=255)
           gCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTc                                                                                                                                         <  1:170390/150‑1 (MQ=255)
               aTCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTa                                                                                                                                     <  2:716698/150‑1 (MQ=255)
               aTCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTa                                                                                                                                     <  1:3097300/150‑1 (MQ=255)
                tCCCTGATGCAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTAc                                                                                                                                    <  2:4185868/150‑1 (MQ=255)
                  ccGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCg                                                                                                                                  >  2:1790923/1‑150 (MQ=255)
                   cGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGt                                                                                                                                 >  2:4861314/1‑150 (MQ=255)
                       tGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTTTGGTCTGGATTCACTACCGTTGGc                                                                                                                             <  2:2878605/150‑1 (MQ=255)
                         tAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCgg                                                                                                                           <  2:4838147/150‑1 (MQ=255)
                         tAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCgg                                                                                                                           <  2:314347/150‑1 (MQ=255)
                            tGGGCTTCAAACTTACCTGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATc                                                                                                                        <  2:3330118/150‑1 (MQ=255)
                            tGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATc                                                                                                                        <  1:1423437/150‑1 (MQ=255)
                                 ttCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGcc                                                                                                                   <  1:205319/150‑1 (MQ=255)
                                  tCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGcca                                                                                                                  <  1:858924/150‑1 (MQ=255)
                                   cAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGccac                                                                                                                 >  1:3304002/1‑150 (MQ=255)
                                    aaaCTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGccacc                                                                                                                >  1:2380227/1‑150 (MQ=255)
                                    aaaCTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGccacc                                                                                                                >  2:4257292/1‑150 (MQ=255)
                                     aaCTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGccacca                                                                                                               >  2:4603151/1‑150 (MQ=255)
                                      aCTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCAt                                                                                                              <  1:1836844/150‑1 (MQ=255)
                                      aCTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCAt                                                                                                              >  1:3967907/1‑150 (MQ=255)
                                       cTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCAtt                                                                                                             <  2:1929705/150‑1 (MQ=255)
                                         tACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGc                                                                                                           <  1:53946/150‑1 (MQ=255)
                                          aCCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGcc                                                                                                          <  1:47532/150‑1 (MQ=255)
                                          aCCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGcc                                                                                                          <  1:3714439/150‑1 (MQ=255)
                                            cGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCaa                                                                                                        >  2:4083613/1‑150 (MQ=255)
                                                  aTTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTc                                                                                                  <  1:2635458/150‑1 (MQ=255)
                                                   ttACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCg                                                                                                 >  2:2594384/1‑150 (MQ=255)
                                                    tACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGc                                                                                                >  1:197397/1‑150 (MQ=255)
                                                     aCGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGcc                                                                                               <  2:4212881/150‑1 (MQ=255)
                                                      cGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCa                                                                                              <  2:356181/150‑1 (MQ=255)
                                                           aGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAata                                                                                         >  2:320496/1‑150 (MQ=255)
                                                             gTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATg                                                                                       >  1:2454354/1‑150 (MQ=255)
                                                             gTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATg                                                                                       <  1:2749746/150‑1 (MQ=255)
                                                                 tACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGc                                                                                   >  2:2476737/1‑150 (MQ=255)
                                                                       cacaGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCtg                                                                             >  1:4197283/1‑150 (MQ=255)
                                                                        acaGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCtgt                                                                            <  1:1247339/150‑1 (MQ=255)
                                                                        acaGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCtgt                                                                            <  1:2089156/150‑1 (MQ=255)
                                                                        acaGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCtgt                                                                            <  2:2765252/150‑1 (MQ=255)
                                                                        acaGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCtgt                                                                            >  1:1298925/1‑150 (MQ=255)
                                                                          aGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTgg                                                                          <  2:2655531/150‑1 (MQ=255)
                                                                          aGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTgg                                                                          <  2:2999670/150‑1 (MQ=255)
                                                                             ccTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGctt                                                                       >  2:708289/1‑150 (MQ=255)
                                                                              cTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGcttc                                                                      <  1:2175153/150‑1 (MQ=255)
                                                                               tGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGcttct                                                                     >  2:3637383/1‑150 (MQ=255)
                                                                               tGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGcttct                                                                     >  2:1069107/1‑150 (MQ=255)
                                                                               tGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGcttat                                                                     <  2:4292196/150‑3 (MQ=255)
                                                                                  ttCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTcc                                                                  <  1:1532089/150‑1 (MQ=255)
                                                                                   tctcACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTccc                                                                 >  1:67334/1‑150 (MQ=255)
                                                                                       aCTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATc                                                             >  1:2262099/1‑150 (MQ=255)
                                                                                        cTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCg                                                            >  2:2109227/1‑150 (MQ=255)
                                                                                         tGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGa                                                           >  2:4593345/1‑150 (MQ=255)
                                                                                          gtgaCCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGAt                                                          <  2:4484795/147‑1 (MQ=255)
                                                                                          gTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGAt                                                          <  2:1496762/150‑1 (MQ=255)
                                                                                          gTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGAt                                                          >  2:4781300/1‑150 (MQ=255)
                                                                                          gTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGAt                                                          <  2:2830313/150‑1 (MQ=255)
                                                                                           ttACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATg                                                         >  1:1026171/1‑150 (MQ=255)
                                                                                            tACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGc                                                        >  1:897801/1‑150 (MQ=255)
                                                                                            tACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACAGTTGGCGGATCGCGCCAACATTGCCAATATGTCGCCAGAATATGTTGCCACCTGTGGCTTCTTCCCAATCGATGc                                                        <  1:1185332/150‑1 (MQ=255)
                                                                                            tACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGc                                                        >  1:4119685/1‑150 (MQ=255)
                                                                                            tACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGc                                                        <  2:3955512/150‑1 (MQ=255)
                                                                                              cccAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTg                                                      <  2:2171589/150‑1 (MQ=255)
                                                                                              cccAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTg                                                      <  1:4660905/150‑1 (MQ=255)
                                                                                               ccAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATAGGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGt                                                     >  1:818818/1‑150 (MQ=255)
                                                                                                cAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTa                                                    >  2:131764/1‑150 (MQ=255)
                                                                                                     gctgcGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCt                                               <  2:1389816/150‑1 (MQ=255)
                                                                                                     gctgcGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCt                                               <  1:34311/150‑1 (MQ=255)
                                                                                                     gctgcGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCt                                               <  1:33599/150‑1 (MQ=255)
                                                                                                           cAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTa                                         <  1:5042788/150‑1 (MQ=255)
                                                                                                               cATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATg                                     >  2:3887761/1‑150 (MQ=255)
                                                                                                                    cGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGttt                                <  2:1124508/150‑1 (MQ=255)
                                                                                                                    cGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGttt                                >  1:3138252/1‑150 (MQ=255)
                                                                                                                     gtggtgGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTa                               >  1:4439277/1‑150 (MQ=255)
                                                                                                                      tggtggGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTaa                              >  1:2205040/1‑150 (MQ=255)
                                                                                                                      tggtggGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTaa                              <  2:2334527/150‑1 (MQ=255)
                                                                                                                         tggGGAAATTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCg                           >  1:2538992/1‑150 (MQ=255)
                                                                                                                               atttcttcGAATTTTTTGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCa                     <  1:3218566/148‑1 (MQ=255)
                                                                                                                               aaTTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCa                     <  1:4981800/150‑1 (MQ=255)
                                                                                                                               aaTTCGTCGAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCa                     >  2:2198381/1‑150 (MQ=255)
                                                                                                                               aaTTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCa                     <  2:2039696/150‑1 (MQ=255)
                                                                                                                                  tcgtcgAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCg                  <  2:4711942/150‑1 (MQ=255)
                                                                                                                                   cgtcgAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGa                 >  1:3822274/1‑150 (MQ=255)
                                                                                                                                   cgtcgAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGa                 >  1:795610/1‑150 (MQ=255)
                                                                                                                                     tcgAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAg               >  1:571798/1‑150 (MQ=255)
                                                                                                                                      cgAATTTTATGGTTATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGa              <  2:4563134/150‑1 (MQ=255)
                                                                                                                                       gAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGAt             >  2:1632641/1‑150 (MQ=255)
                                                                                                                                        aaTTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATc            <  1:2873321/150‑1 (MQ=255)
                                                                                                                                            ttATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGt        >  2:5133440/1‑150 (MQ=255)
                                                                                                                                             tATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTc       <  1:1467093/150‑1 (MQ=255)
                                                                                                                                                ggTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAg    >  1:3400143/1‑150 (MQ=255)
                                                                                                                                                  tGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGtt  <  2:2396702/150‑1 (MQ=255)
                                                                                                                                                  tGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGtt  <  2:4940910/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTT  >  CP009273/1330808‑1331103

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: