Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,769,440 C→A 46.7% A166A (GCC→GCA aroD → 3‑dehydroquinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,769,4400CA46.7% 184.3 / inf 263A166A (GCC→GCAaroD3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (71/69);  new base A (57/66);  total (128/135)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 271 total aligned reads displayed.

TTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGG  >  CP009273/1769291‑1769589
                                                                                                                                                     |                                                                                                                                                     
ttATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGcc                                                                                                                                                       <  2:610497/150‑1 (MQ=255)
      aCCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATAtt                                                                                                                                                 >  2:53269/1‑150 (MQ=255)
         ggTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCt                                                                                                                                              <  1:1472763/150‑1 (MQ=255)
             atgatCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGa                                                                                                                                          >  1:1842764/1‑150 (MQ=255)
               gatCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGAtt                                                                                                                                        <  2:2121565/150‑1 (MQ=255)
               gatCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGAtt                                                                                                                                        >  2:719515/1‑150 (MQ=255)
               gatCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGAtt                                                                                                                                        <  2:608037/150‑1 (MQ=255)
                  cAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTgcg                                                                                                                                     >  1:4351269/1‑150 (MQ=255)
                   aGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTgcgc                                                                                                                                    >  1:3875685/1‑150 (MQ=255)
                   aGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTgcgc                                                                                                                                    >  1:1740313/1‑150 (MQ=255)
                    ggTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCt                                                                                                                                   <  2:3765332/150‑1 (MQ=255)
                     gTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTg                                                                                                                                  >  1:3744196/1‑150 (MQ=255)
                        aaaGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATg                                                                                                                               <  1:1791929/150‑1 (MQ=255)
                        aaaGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATg                                                                                                                               <  1:1245675/150‑1 (MQ=255)
                        aaaGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATg                                                                                                                               <  2:4179600/150‑1 (MQ=255)
                           gAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCg                                                                                                                            >  1:4249474/1‑150 (MQ=255)
                             aaCCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCATCGACGCCGATATTCCAAAGATTGCGCTGATGCCGCa                                                                                                                          >  1:2299705/1‑150 (MQ=255)
                              aCCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCaa                                                                                                                         <  1:3734542/150‑1 (MQ=255)
                              aCCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCaa                                                                                                                         <  1:32704/150‑1 (MQ=255)
                              aCCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCACAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCaa                                                                                                                         >  1:3538398/1‑150 (MQ=255)
                                 gTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGt                                                                                                                      <  2:1596790/150‑1 (MQ=255)
                                  tCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTa                                                                                                                     <  1:3086918/150‑1 (MQ=255)
                                    gCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTAcc                                                                                                                   <  1:3831035/150‑1 (MQ=255)
                                    gCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTAcc                                                                                                                   >  1:2591986/1‑150 (MQ=255)
                                     ccTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCa                                                                                                                  <  1:863209/150‑1 (MQ=255)
                                      cTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAg                                                                                                                 <  2:2212957/150‑1 (MQ=255)
                                       tACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGc                                                                                                                >  2:2842952/1‑150 (MQ=255)
                                       tACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGc                                                                                                                <  2:2627388/150‑1 (MQ=255)
                                        aCGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCg                                                                                                               >  1:3132700/1‑150 (MQ=255)
                                        aCGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCg                                                                                                               >  2:4025236/1‑150 (MQ=255)
                                         cGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGa                                                                                                              >  2:2137858/1‑150 (MQ=255)
                                          gCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGAt                                                                                                             <  2:1083788/150‑1 (MQ=255)
                                           cccACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGAtg                                                                                                            >  2:4475332/1‑150 (MQ=255)
                                             cACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGAtgtg                                                                                                          >  1:1072710/1‑150 (MQ=255)
                                                 cgcATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGa                                                                                                      >  2:4051352/1‑150 (MQ=255)
                                                  gcATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGAc                                                                                                     <  2:114130/150‑1 (MQ=255)
                                                    atgatgTGAAAGTAGTCATGTCCAACCATGACTTCCATAAATCGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGt                                                                                                   >  1:1112157/1‑150 (MQ=255)
                                                    atgatgTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGt                                                                                                   <  1:3139724/150‑1 (MQ=255)
                                                     tgatgTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGtt                                                                                                  >  1:1583993/1‑150 (MQ=255)
                                                      gatgTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGttg                                                                                                 >  2:2900454/1‑150 (MQ=255)
                                                       atgTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGttgc                                                                                                >  1:3335032/1‑150 (MQ=255)
                                                         gtgAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGttgctt                                                                                              >  2:429380/1‑150 (MQ=255)
                                                                tagtCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGa                                                                                       <  1:2563227/150‑1 (MQ=255)
                                                                 agtCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGAc                                                                                      <  2:2463208/150‑1 (MQ=255)
                                                                 agtCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGAc                                                                                      <  2:2463462/150‑1 (MQ=255)
                                                                      tGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTgg                                                                                 >  1:3122331/1‑150 (MQ=255)
                                                                       gTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGga                                                                                <  2:4425023/150‑1 (MQ=255)
                                                                        tCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGgag                                                                               <  1:4570332/150‑1 (MQ=255)
                                                                         ccAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGgaga                                                                              >  1:3993320/1‑150 (MQ=255)
                                                                         ccAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGgaga                                                                              >  2:3048791/1‑150 (MQ=255)
                                                                            aCCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGc                                                                           <  1:2241118/150‑1 (MQ=255)
                                                                              cATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAg                                                                         >  2:3877480/1‑150 (MQ=255)
                                                                                tGACTTCCATAAAACGCCGGTAGCCGAAGATATCTTTGCCCGTCTGCGCAAAATGCAATCCTTCGTGTCAGATATTCCTAATTTTGCGCTGTTGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTTCCGCGACCCTGGATATGCAgta                                                                       <  2:2379047/150‑3 (MQ=255)
                                                                                 gACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAg                                                                      >  1:4060665/1‑150 (MQ=255)
                                                                                 gACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAg                                                                      >  2:770093/1‑150 (MQ=255)
                                                                                   cTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCa                                                                    >  1:2124467/1‑150 (MQ=255)
                                                                                    ttCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAg                                                                   >  1:4142336/1‑150 (MQ=255)
                                                                                       cATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTAt                                                                >  2:1160682/1‑150 (MQ=255)
                                                                                       cATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTAt                                                                <  1:4072704/150‑1 (MQ=255)
                                                                                        aTAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATg                                                               <  2:1269464/150‑1 (MQ=255)
                                                                                          aaaaCGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGcc                                                             >  2:4438280/1‑150 (MQ=255)
                                                                                          aaaaCGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGcc                                                             >  2:252352/1‑150 (MQ=255)
                                                                                          aaaaCGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGcc                                                             <  2:1867025/150‑1 (MQ=255)
                                                                                            aaCGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGa                                                           >  1:47651/1‑150 (MQ=255)
                                                                                                  ggAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTcc                                                     <  1:3918499/150‑1 (MQ=255)
                                                                                                    aaGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCaa                                                   <  2:3162118/150‑1 (MQ=255)
                                                                                                     aGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAAt                                                  <  2:3378946/150‑1 (MQ=255)
                                                                                                     aGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAAt                                                  <  1:16833/150‑1 (MQ=255)
                                                                                                      gCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAAtt                                                 <  1:2759720/150‑1 (MQ=255)
                                                                                                      gCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGACGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCAGCGACCATGGAGATGCAGGAGCAGTATGAAGATCGTCCAAtt                                                 >  1:1189485/1‑150 (MQ=255)
                                                                                                        cGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTAt                                               <  1:735453/150‑1 (MQ=255)
                                                                                                         gaagaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATc                                              >  2:1159948/1‑150 (MQ=255)
                                                                                                         gaagaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATc                                              <  2:568914/150‑1 (MQ=255)
                                                                                                          aagaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCa                                             <  1:4020756/150‑1 (MQ=255)
                                                                                                           agaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCAc                                            >  2:465577/1‑150 (MQ=255)
                                                                                                            gaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACg                                           <  2:1377384/150‑1 (MQ=255)
                                                                                                               aTCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATg                                        <  2:1660754/150‑1 (MQ=255)
                                                                                                                tCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGt                                       >  2:2785082/1‑150 (MQ=255)
                                                                                                                 cATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTc                                      <  2:742749/150‑1 (MQ=255)
                                                                                                                 cATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTc                                      >  1:1006717/1‑150 (MQ=255)
                                                                                                                 cATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTc                                      >  1:1006502/1‑150 (MQ=255)
                                                                                                                    tGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGTAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGAt                                   >  2:1299248/1‑150 (MQ=255)
                                                                                                                    tGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGAt                                   <  2:4399727/150‑1 (MQ=255)
                                                                                                                       ccGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGc                                <  2:3463000/150‑1 (MQ=255)
                                                                                                                       ccGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGc                                <  1:207950/150‑1 (MQ=255)
                                                                                                                        cGTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCa                               >  2:2054559/1‑150 (MQ=255)
                                                                                                                         gTCTGCGCAAAATGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCaa                              >  1:578962/1‑150 (MQ=255)
                                                                                                                           cTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCaaaa                            >  1:2690400/1‑150 (MQ=255)
                                                                                                                                cAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTgg                       >  1:2659780/1‑150 (MQ=255)
                                                                                                                                 aaaaTGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGc                      >  2:2995299/1‑150 (MQ=255)
                                                                                                                                 aaaaTGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGc                      >  1:3222656/1‑150 (MQ=255)
                                                                                                                                   aaTGCAATCCTTCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGt                    >  1:953618/1‑150 (MQ=255)
                                                                                                                                     tGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTaa                  >  2:844270/1‑150 (MQ=255)
                                                                                                                                      gCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAAt                 <  1:2776216/150‑1 (MQ=255)
                                                                                                                                        aaTCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAAttt               >  2:407911/1‑150 (MQ=255)
                                                                                                                                           ccTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTtctc            <  1:3177736/150‑1 (MQ=255)
                                                                                                                                             ttCGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGt          <  1:645463/150‑1 (MQ=255)
                                                                                                                                               cGACGCAGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGACAAAAACTGGCGTAATTTCTCGTct        <  2:1952553/150‑1 (MQ=255)
                                                                                                                                                gACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctg       <  2:1815367/150‑1 (MQ=255)
                                                                                                                                                     cGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggctgg  >  2:1460115/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAGTAGTCATGTCCAACCATGACTTCCATAAAACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGG  >  CP009273/1769291‑1769589

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: