Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 724,813 A→T 100% intergenic (+17/‑226) ybfA → / → rhsC DUF2517 family protein/Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273724,8130AT95.9% 876.7 / 20.4 267intergenic (+17/‑226)ybfA/rhsCDUF2517 family protein/Rhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base A (3/8);  new base T (100/156);  total (103/164)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 267 total aligned reads displayed.

GCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAT  >  CP009273/724667‑724959
                                                                                                                                                  |                                                                                                                                                  
gcTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCAt                                                                                                                                                 >  1:1892028/1‑150 (MQ=255)
  tGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTg                                                                                                                                               <  2:3147665/150‑1 (MQ=255)
   gCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTCCAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAATCCGGTGTGGATGGTTCAGGACTAAAAGTCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGa                                                                                                                                              <  1:571654/150‑1 (MQ=255)
   gCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGATAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGa                                                                                                                                              <  1:612175/150‑1 (MQ=255)
   gCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGa                                                                                                                                              <  1:571179/150‑1 (MQ=255)
       ttCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAAc                                                                                                                                          <  1:4110889/150‑1 (MQ=255)
       ttCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAAc                                                                                                                                          <  1:2268386/150‑1 (MQ=255)
            gCTTTTCTGGAAAGACCGCGCCCGCATTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCa                                                                                                                                     <  1:2561715/150‑1 (MQ=255)
                   tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaaa                                                                                                                              >  2:2581789/1‑150 (MQ=255)
                   tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaaa                                                                                                                              >  2:2919248/1‑150 (MQ=255)
                   tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaaa                                                                                                                              <  1:3371751/150‑1 (MQ=255)
                   tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAAACGCCAGCGGaaa                                                                                                                              >  2:3893014/1‑150 (MQ=255)
                    ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAAc                                                                                                                             <  1:1000686/150‑1 (MQ=255)
                      aaaGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGc                                                                                                                           >  2:3815539/1‑150 (MQ=255)
                         gACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTgg                                                                                                                        <  2:2382304/150‑1 (MQ=255)
                              cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                   <  2:2335748/150‑1 (MQ=255)
                              cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                   >  2:1422181/1‑150 (MQ=255)
                               gcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGttt                                                                                                                  >  2:1537807/1‑150 (MQ=255)
                                  cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                               <  1:2406447/150‑1 (MQ=255)
                                   gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTg                                                                                                              <  1:4337229/150‑1 (MQ=255)
                                    cTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGt                                                                                                             <  2:2701147/150‑1 (MQ=255)
                                      tttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTAt                                                                                                           >  2:3508131/1‑150 (MQ=255)
                                       ttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATc                                                                                                          <  1:2233288/150‑1 (MQ=255)
                                        tACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCa                                                                                                         <  1:3743255/150‑1 (MQ=255)
                                        tACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCa                                                                                                         <  1:792786/150‑1 (MQ=255)
                                            gCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGaa                                                                                                     <  1:3301278/150‑1 (MQ=255)
                                            gCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGaa                                                                                                     <  1:3118250/150‑1 (MQ=255)
                                             cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAg                                                                                                    <  1:3496007/150‑1 (MQ=255)
                                               aCCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAt                                                                                                  <  1:4265643/150‑1 (MQ=255)
                                                ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATg                                                                                                 <  1:188535/150‑1 (MQ=255)
                                                    cATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCa                                                                                             >  1:3010557/1‑150 (MQ=255)
                                                        gcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAgg                                                                                         <  1:52141/150‑1 (MQ=255)
                                                        gcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAgg                                                                                         >  2:3911697/1‑150 (MQ=255)
                                                          gTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGaa                                                                                       <  1:2391652/150‑1 (MQ=255)
                                                          gTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGaa                                                                                       <  1:2391639/150‑1 (MQ=255)
                                                             tGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAatat                                                                                    <  2:4337470/150‑1 (MQ=255)
                                                             tGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTGTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTATTTTGTATCAGGAAGATGTTCAGAGGAatat                                                                                    <  2:401481/150‑1 (MQ=255)
                                                               gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                  >  1:4492353/1‑150 (MQ=255)
                                                                 cGAAAACCAGCGAGAAACCGGTGTGGATGGATCAGGCCGAAACGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGTAGATGTTCAGAGGAATATTTAg                                                                                <  2:2599963/150‑1 (MQ=255)
                                                                  gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                               >  2:737510/1‑150 (MQ=255)
                                                                   aaaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGAc                                                                              >  2:2469211/1‑150 (MQ=255)
                                                                   aaaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGAc                                                                              >  2:4503180/1‑150 (MQ=255)
                                                                   aaaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGAc                                                                              >  1:998443/1‑150 (MQ=255)
                                                                      aCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATc                                                                           <  2:1892028/150‑1 (MQ=255)
                                                                       ccAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGACGATGTTCAGAGGAATATTTAGACATCg                                                                          >  2:3674280/1‑150 (MQ=255)
                                                                         aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGtt                                                                        >  1:1827231/1‑150 (MQ=255)
                                                                           cGATAAACCGGTGTGGATGGATCAGGCCGACAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTtat                                                                      <  2:3839837/150‑1 (MQ=255)
                                                                           cGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTtat                                                                      >  1:969866/1‑150 (MQ=255)
                                                                             aTAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAc                                                                    <  2:3066320/150‑1 (MQ=255)
                                                                                 aCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAAt                                                                <  2:844862/150‑1 (MQ=255)
                                                                                  ccGGTGTTGATTGATCAGGCGGATAAGGCAACTGGTTATTTTTATTGATTGGTGAGTTTTGTATTCATTGAAAACCTCCAGCGGTAACGCTGGCGGTTTTTTTTATCAGTAAGATGTTCAGTGGAATATTTAGACATCGTTATACCAATc                                                               <  1:4268853/150‑1 (MQ=255)
                                                                                   cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCg                                                              <  2:3877303/150‑1 (MQ=255)
                                                                                    ggTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGa                                                             >  1:742306/1‑150 (MQ=255)
                                                                                      tgtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                           <  1:915132/150‑1 (MQ=255)
                                                                                      tgtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                           <  1:1433150/150‑1 (MQ=255)
                                                                                        tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACa                                                         <  2:4138479/150‑1 (MQ=255)
                                                                                        tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACa                                                         <  2:4332052/150‑1 (MQ=255)
                                                                                        tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACa                                                         <  2:505629/150‑1 (MQ=255)
                                                                                         ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                        <  1:4466322/150‑1 (MQ=255)
                                                                                         ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                        <  2:1298891/150‑1 (MQ=255)
                                                                                            tggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatata                                                     <  1:35318/150‑1 (MQ=255)
                                                                                             ggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatatat                                                    <  1:1648895/150‑1 (MQ=255)
                                                                                              gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAtt                                                   <  2:4251673/150‑1 (MQ=255)
                                                                                                tCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCa                                                 >  2:1493844/1‑150 (MQ=255)
                                                                                                 cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                <  2:1485950/150‑1 (MQ=255)
                                                                                                 cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                <  1:336845/150‑1 (MQ=255)
                                                                                                 cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                >  1:3515498/1‑150 (MQ=255)
                                                                                                 cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                <  2:2988449/150‑1 (MQ=255)
                                                                                                       gAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAata                                          >  2:444049/1‑150 (MQ=255)
                                                                                                       gAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAata                                          >  2:1079846/1‑150 (MQ=255)
                                                                                                          aaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatat                                       <  1:1547293/150‑1 (MQ=255)
                                                                                                              cAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGTAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATaa                                   <  2:4173661/150‑1 (MQ=255)
                                                                                                               aaCTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGAGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATaaa                                  <  1:576933/150‑1 (MQ=255)
                                                                                                                 cTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAata                                <  1:1673128/150‑1 (MQ=255)
                                                                                                                 cTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAata                                >  2:3375146/1‑150 (MQ=255)
                                                                                                                 cTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAata                                >  2:47248/1‑150 (MQ=255)
                                                                                                                 cTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAata                                >  2:3375210/1‑150 (MQ=255)
                                                                                                                  tGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGACACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAatat                               <  1:3106331/150‑1 (MQ=255)
                                                                                                                   ggTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAtt                              >  1:1305684/1‑150 (MQ=255)
                                                                                                                     tGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAtttt                            >  1:4147710/1‑150 (MQ=255)
                                                                                                                        tttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCt                         >  2:1682416/1‑150 (MQ=255)
                                                                                                                        tttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCt                         >  2:3243905/1‑150 (MQ=255)
                                                                                                                        tttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCt                         <  1:3656135/150‑1 (MQ=255)
                                                                                                                          tttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTaa                       <  1:3890676/150‑1 (MQ=255)
                                                                                                                          tttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTaa                       >  2:589974/1‑150 (MQ=255)
                                                                                                                           ttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAAt                      >  1:1287884/1‑150 (MQ=255)
                                                                                                                            tATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAAtt                     <  1:1088090/150‑1 (MQ=255)
                                                                                                                            tATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAAtt                     <  1:710116/150‑1 (MQ=255)
                                                                                                                            tATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATAGAAACATATATTCATGAAATATATATAAATATTTTCCTAAtt                     >  1:25267/1‑150 (MQ=255)
                                                                                                                             attgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTg                    >  2:3883364/1‑150 (MQ=255)
                                                                                                                             attgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTg                    >  2:3518619/1‑150 (MQ=255)
                                                                                                                             attgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTg                    >  2:3518417/1‑150 (MQ=255)
                                                                                                                             attgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTg                    >  2:3517541/1‑150 (MQ=255)
                                                                                                                               tgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGtt                  <  1:4195211/150‑1 (MQ=255)
                                                                                                                                  ttgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCtt               <  1:237939/150‑1 (MQ=255)
                                                                                                                                     cTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATc            >  2:4287148/1‑150 (MQ=255)
                                                                                                                                     cTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATc            >  2:2040377/1‑150 (MQ=255)
                                                                                                                                      tGAGATTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCt           <  2:617975/150‑1 (MQ=255)
                                                                                                                                        aGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGa         >  1:4005722/1‑150 (MQ=255)
                                                                                                                                               aaTTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGatat  >  1:2035070/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAT  >  CP009273/724667‑724959

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: