Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 435,336 G→T 100% Q100K (CAG→AAG)  dxs ← 1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273435,3360GT98.2% 928.7 / 7.1 275Q100K (CAG→AAG) dxs1‑deoxyxylulose‑5‑phosphate synthase, thiamine triphosphate‑binding, FAD‑requiring
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base T (147/123);  total (149/126)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 275 total aligned reads displayed.

CACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGCCCAGCCCGG  >  CP009273/435189‑435479
                                                                                                                                                   |                                                                                                                                               
  caGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACgg                                                                                                                                             >  1:4518135/1‑150 (MQ=255)
    gACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGAt                                                                                                                                           <  2:3014426/150‑1 (MQ=255)
      cGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATgg                                                                                                                                         <  1:2111440/150‑1 (MQ=255)
      cGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATgg                                                                                                                                         <  1:1833578/150‑1 (MQ=255)
      cGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATgg                                                                                                                                         >  2:1531871/1‑150 (MQ=255)
       ggTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGt                                                                                                                                        <  1:3196416/150‑1 (MQ=255)
        gTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTg                                                                                                                                       <  2:3431095/150‑1 (MQ=255)
         tGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGc                                                                                                                                      >  1:2929227/1‑150 (MQ=255)
              cgATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACTCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGAtt                                                                                                                                 <  2:1789729/150‑1 (MQ=255)
               gATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGAttt                                                                                                                                <  2:2044218/150‑1 (MQ=255)
                      gCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTcgcg                                                                                                                         <  1:3517975/150‑1 (MQ=255)
                      gCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTcgcg                                                                                                                         >  2:1425133/1‑150 (MQ=255)
                       ccTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCgcgg                                                                                                                        >  1:4054017/1‑150 (MQ=255)
                         ttcttTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCgcggcg                                                                                                                      >  2:2084744/1‑150 (MQ=255)
                          tcttTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGt                                                                                                                     >  1:3430963/1‑150 (MQ=255)
                            tttttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTcc                                                                                                                   <  1:3163313/150‑1 (MQ=255)
                            tttttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTcc                                                                                                                   <  1:4235437/150‑1 (MQ=255)
                            tttttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTcc                                                                                                                   <  2:3674587/150‑1 (MQ=255)
                              tttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCgg                                                                                                                 <  1:563812/150‑1 (MQ=255)
                              tttCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCgg                                                                                                                 <  1:4010690/150‑1 (MQ=255)
                                tCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTc                                                                                                               >  1:3291944/1‑150 (MQ=255)
                                 cGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCa                                                                                                              <  2:3010943/150‑1 (MQ=255)
                                   gcagcaACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATCTTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCaaa                                                                                                            <  2:406167/150‑1 (MQ=255)
                                     agcaACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAAt                                                                                                          >  2:3482509/1‑150 (MQ=255)
                                     agcaACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAAt                                                                                                          <  1:1119966/150‑1 (MQ=255)
                                     agcaACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAAt                                                                                                          >  2:3337619/1‑150 (MQ=255)
                                       caACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAAttt                                                                                                        >  1:1012145/1‑150 (MQ=255)
                                        aaCCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAAtttt                                                                                                       >  2:400407/1‑150 (MQ=255)
                                         aCCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTa                                                                                                      <  2:3380907/150‑1 (MQ=255)
                                          ccGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTAt                                                                                                     <  1:3686908/150‑1 (MQ=255)
                                          ccGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTAt                                                                                                     <  1:94/150‑1 (MQ=255)
                                            gCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGACCGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGc                                                                                                   >  2:178644/1‑150 (MQ=255)
                                            gCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGc                                                                                                   <  1:3868531/150‑1 (MQ=255)
                                              aaTACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCgg                                                                                                 >  1:664765/1‑150 (MQ=255)
                                               aTACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGa                                                                                                >  2:4444271/1‑150 (MQ=255)
                                               aTACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGa                                                                                                >  2:49410/1‑150 (MQ=255)
                                               aTACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGa                                                                                                <  2:2535503/150‑1 (MQ=255)
                                                tACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGAt                                                                                               <  1:361537/150‑1 (MQ=255)
                                                 aCCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATa                                                                                              >  2:1236371/1‑150 (MQ=255)
                                                   cAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAg                                                                                            <  2:703740/150‑1 (MQ=255)
                                                     aTTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGcc                                                                                          >  2:3288177/1‑150 (MQ=255)
                                                         cGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGAt                                                                                      >  2:3228949/1‑150 (MQ=255)
                                                         cGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGAt                                                                                      <  2:3516811/150‑1 (MQ=255)
                                                         cGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGAt                                                                                      <  1:2692939/150‑1 (MQ=255)
                                                             aCTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGccc                                                                                  <  2:4571129/150‑1 (MQ=255)
                                                             aCTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGccc                                                                                  <  2:598332/150‑1 (MQ=255)
                                                                 aTGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCAGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCcaca                                                                              >  2:3787648/1‑150 (MQ=255)
                                                                  tGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACAt                                                                             <  1:4522432/150‑1 (MQ=255)
                                                                  tGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACAt                                                                             <  2:4467508/150‑1 (MQ=255)
                                                                       gTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCaa                                                                        >  1:4078685/1‑150 (MQ=255)
                                                                         tgatgaATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAAt                                                                      >  2:2582868/1‑150 (MQ=255)
                                                                            tgaATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCaa                                                                   >  1:1049119/1‑150 (MQ=255)
                                                                             gaATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAAt                                                                  <  1:3359888/150‑1 (MQ=255)
                                                                             gaATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAAt                                                                  <  1:3445941/150‑1 (MQ=255)
                                                                              aaTGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAAtt                                                                 <  1:2347421/150‑1 (MQ=255)
                                                                              aaTGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAAtt                                                                 <  1:4493806/150‑1 (MQ=255)
                                                                               aTGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTg                                                                >  1:827506/1‑150 (MQ=255)
                                                                                 gCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCAGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGt                                                              <  2:2869469/150‑1 (MQ=255)
                                                                                    cGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGACCGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCaa                                                           >  2:1504465/1‑150 (MQ=255)
                                                                                    cGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCaa                                                           >  2:608067/1‑150 (MQ=255)
                                                                                      aCGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAAc                                                         >  1:2765623/1‑150 (MQ=255)
                                                                                       cGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACg                                                        >  2:2507568/1‑150 (MQ=255)
                                                                                        gCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACgg                                                       >  2:2083776/1‑150 (MQ=255)
                                                                                        gCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACgg                                                       >  2:199665/1‑150 (MQ=255)
                                                                                          ttAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACgggg                                                     <  1:4452308/150‑1 (MQ=255)
                                                                                          ttAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACgggg                                                     <  1:2195214/150‑1 (MQ=255)
                                                                                           tAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGgt                                                    >  2:306508/1‑150 (MQ=255)
                                                                                             aTACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGgtgt                                                  <  1:3927185/150‑1 (MQ=255)
                                                                                                   cATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGa                                            >  2:4497480/1‑150 (MQ=255)
                                                                                                     tatTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACa                                          >  1:1409267/1‑150 (MQ=255)
                                                                                                      atTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACAt                                         >  1:1768516/1‑150 (MQ=255)
                                                                                                      atTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACAt                                         >  1:2716371/1‑150 (MQ=255)
                                                                                                         cGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGt                                      >  2:2170756/1‑150 (MQ=255)
                                                                                                         cGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGt                                      >  2:3002104/1‑150 (MQ=255)
                                                                                                         cGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGt                                      <  2:2883962/150‑1 (MQ=255)
                                                                                                         cGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGt                                      >  2:4135648/1‑150 (MQ=255)
                                                                                                          gCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTg                                     >  1:3970749/1‑150 (MQ=255)
                                                                                                          gCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTg                                     <  1:4205794/150‑1 (MQ=255)
                                                                                                             ttCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAg                                  >  2:2932272/1‑150 (MQ=255)
                                                                                                             ttCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAg                                  <  2:3575651/150‑1 (MQ=255)
                                                                                                               cgccgcGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAgcg                                >  1:2299881/1‑150 (MQ=255)
                                                                                                               cgccgcGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAgcg                                <  2:1085168/150‑1 (MQ=255)
                                                                                                                gccgcGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAgcgc                               >  2:3234171/1‑150 (MQ=255)
                                                                                                                 ccgcGCCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGcc                              <  2:1005399/150‑1 (MQ=255)
                                                                                                                   gcgcCACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCAc                            >  1:4620926/1‑150 (MQ=255)
                                                                                                                      ccACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGt                         <  1:1383441/150‑1 (MQ=255)
                                                                                                                       cACGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTc                        >  1:1463971/1‑150 (MQ=255)
                                                                                                                       aaCGGGAACGTGTGCAGAACGCAATTCATAAGCATTGAGCCGATTTTGTCGCGGCGTCCGGTCAAAAATTTATGCGGATAAGCCTGATGCCACACATCAAAAAACAAATTTACAAACTAGTTGATGTATAAATAGTGCATCGACACgtac                        >  2:643359/2‑147 (MQ=255)
                                                                                                                         cGGGAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAg                      >  1:3187306/1‑150 (MQ=255)
                                                                                                                          gggAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGt                     >  1:4342855/1‑150 (MQ=255)
                                                                                                                          gggAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGt                     >  1:1484814/1‑150 (MQ=255)
                                                                                                                           ggAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGtt                    >  2:3872757/1‑150 (MQ=255)
                                                                                                                            gtaCGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGTTTTTGTCGCGGCGTCCGGTCATATTTTTATGCGGATTAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTc                   <  2:2641119/148‑1 (MQ=255)
                                                                                                                            gAACGGGTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTc                   <  2:2342975/150‑1 (MQ=255)
                                                                                                                                 ggTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCg              >  1:1256218/1‑150 (MQ=255)
                                                                                                                                 ggTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCg              >  1:3976777/1‑150 (MQ=255)
                                                                                                                                  gTGCAGACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGt             <  1:3482505/150‑1 (MQ=255)
                                                                                                                                       gACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGCCCa        >  2:1304315/1‑150 (MQ=255)
                                                                                                                                       gACCGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGACAAAATATTATGCGGATAAGCCAGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGCCCa        >  1:1868869/1‑150 (MQ=255)
                                                                                                                                          cGCCTTTCTTACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGCCCAGcc     >  2:3072694/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CACAGACGGTGCGGCGATTTTTGCCTTCTTTTTCGGCAGCAACCGCAATACCAATTCCGGCACTGATGGAGGTTGATGAATGCCCGACGCTTAATACGTCATATTCGCTTTCGCCGCGCCACGGGAACGGGTGCAGACCGCCTTTCTGACGGATGGTGCCGATTTTGTCGCGGCGTCCGGTCAAAATTTTATGCGGATAAGCCTGATGCCCCACATCCCAAATCAATTGGTCAAACGGGGTGTTGTAGACATAGTGCAGCGCCACGGTCAGTTCGACCGTGCCCAGCCCGG  >  CP009273/435189‑435479

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