Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,322,609 C→A 100% W223L (TGG→TTG)  yfaQ ← tandem DUF2300 domain protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,322,6090CA99.0% 950.5 / ‑0.6 292W223L (TGG→TTG) yfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base A (135/154);  total (138/154)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 298 total aligned reads displayed.

GCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACAT  >  CP009273/2322462‑2322756
                                                                                                                                                   |                                                                                                                                                   
  tgtgATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAg                                                                                                                                                 >  1:1399600/1‑150 (MQ=255)
  tgtgATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAg                                                                                                                                                 >  1:4884220/1‑150 (MQ=255)
   gtgATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGc                                                                                                                                                >  1:1517600/1‑150 (MQ=255)
   gtgATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGCGCCACTGCCAGTAACCAACAGc                                                                                                                                                <  2:1008666/150‑1 (MQ=255)
      aTTGTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgt                                                                                                                                             <  2:1638982/150‑1 (MQ=255)
       ttttCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtg                                                                                                                                            >  2:1159977/1‑150 (MQ=255)
        tttCTTCAGCGGATAGCGGGCGAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtgt                                                                                                                                           >  1:2316761/1‑150 (MQ=255)
           cttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTAcc                                                                                                                                        >  2:769598/1‑150 (MQ=255)
           cttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTAcc                                                                                                                                        <  2:1737706/150‑1 (MQ=255)
             tcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCAt                                                                                                                                      <  1:4776515/150‑1 (MQ=255)
             tcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCAt                                                                                                                                      <  1:3197376/150‑1 (MQ=255)
                  ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTGGTACCGTTTTGCTGGTGCCACTGCCAGTACCCAACAGCGGTGTACCAGCGGtt                                                                                                                                 >  2:2064318/1‑150 (MQ=255)
                  ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGtt                                                                                                                                 >  2:1232478/1‑150 (MQ=255)
                  ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGtt                                                                                                                                 >  1:1205945/1‑150 (MQ=255)
                   gATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTa                                                                                                                                >  1:2460166/1‑150 (MQ=255)
                    aTAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTaa                                                                                                                               >  2:3782612/1‑150 (MQ=255)
                      aGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAca                                                                                                                             <  2:3256886/150‑2 (MQ=255)
                         gggCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagc                                                                                                                          >  1:3426330/1‑150 (MQ=255)
                           gCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagcca                                                                                                                        <  2:3056313/150‑1 (MQ=255)
                           gCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagcca                                                                                                                        <  1:2743374/150‑1 (MQ=255)
                           gCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGACTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagcca                                                                                                                        <  2:2508411/150‑1 (MQ=255)
                              aaCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCa                                                                                                                     <  1:1043581/150‑1 (MQ=255)
                               aCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaa                                                                                                                    <  1:28938/150‑1 (MQ=255)
                               aCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaa                                                                                                                    <  1:677043/150‑1 (MQ=255)
                                cAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaaa                                                                                                                   >  1:1962229/1‑150 (MQ=255)
                                 aGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaaaa                                                                                                                  >  1:5074560/1‑150 (MQ=255)
                                   cTCCACCTCTCCACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAacc                                                                                                                <  2:1042616/150‑1 (MQ=255)
                                     ccACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccac                                                                                                              >  1:4047379/1‑150 (MQ=255)
                                       aCCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCAccc                                                                                                            <  1:1795116/150‑1 (MQ=255)
                                        ccTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGTTTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCATAACCACCCt                                                                                                           <  2:1635727/150‑1 (MQ=255)
                                          tctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGa                                                                                                         <  2:3403623/150‑1 (MQ=255)
                                          tctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGa                                                                                                         <  1:4287540/150‑1 (MQ=255)
                                           ctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGAc                                                                                                        <  1:251228/150‑1 (MQ=255)
                                               acaTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAt                                                                                                    >  2:3661958/1‑150 (MQ=255)
                                                 aTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCATGACGGAttt                                                                                                  >  1:223037/1‑150 (MQ=255)
                                                      ccACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGt                                                                                             >  2:441641/1‑150 (MQ=255)
                                                      ccACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGt                                                                                             <  2:14620/150‑1 (MQ=255)
                                                          tGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTat                                                                                         >  1:4492782/1‑150 (MQ=255)
                                                              aaCCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcata                                                                                     <  2:4844306/150‑2 (MQ=255)
                                                               aCCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCg                                                                                    <  2:2286744/150‑1 (MQ=255)
                                                                ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCgg                                                                                   <  2:3338776/150‑1 (MQ=255)
                                                                ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAGTTCCTGTTTATCATCgg                                                                                   <  2:563732/150‑1 (MQ=255)
                                                                  tGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCa                                                                                 >  1:2192865/1‑150 (MQ=255)
                                                                  tGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCa                                                                                 <  2:4076773/150‑1 (MQ=255)
                                                                  tGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCa                                                                                 <  1:5006620/150‑1 (MQ=255)
                                                                   gCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCaa                                                                                >  2:2316759/1‑150 (MQ=255)
                                                                     gTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAAc                                                                              <  1:985787/150‑1 (MQ=255)
                                                                     gTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAAc                                                                              <  1:4315398/150‑1 (MQ=255)
                                                                        ggCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAg                                                                           >  2:1861303/1‑150 (MQ=255)
                                                                           aCTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTc                                                                        <  2:3846367/150‑1 (MQ=255)
                                                                           aCTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTc                                                                        >  2:4712869/1‑150 (MQ=255)
                                                                            cTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTcc                                                                       >  2:4049800/1‑150 (MQ=255)
                                                                             tGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACATCGGTGTACCATCGTTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                      <  2:1616599/150‑1 (MQ=255)
                                                                             tGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                      <  1:3635438/150‑1 (MQ=255)
                                                                               gCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAgg                                                                    <  1:4464143/150‑1 (MQ=255)
                                                                                cAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACGGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGt                                                                   <  2:649306/150‑1 (MQ=255)
                                                                                 aaaaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtt                                                                  <  1:4098514/150‑1 (MQ=255)
                                                                                       cGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAcc                                                            >  2:4276883/1‑150 (MQ=255)
                                                                                           tAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTa                                                        >  1:2561454/1‑150 (MQ=255)
                                                                                             aGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTAcc                                                      <  2:4770007/150‑1 (MQ=255)
                                                                                              gCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTAccc                                                     <  1:3916546/150‑1 (MQ=255)
                                                                                               cACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCg                                                    >  1:171575/1‑150 (MQ=255)
                                                                                                  cGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACt                                                 >  2:3641009/1‑150 (MQ=255)
                                                                                                  cGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATGTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACt                                                 >  2:2722171/1‑150 (MQ=255)
                                                                                                   gTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTa                                                >  1:1894887/1‑150 (MQ=255)
                                                                                                    tCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAc                                               >  1:2569596/1‑150 (MQ=255)
                                                                                                    tCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAc                                               >  2:4088254/1‑150 (MQ=255)
                                                                                                    tCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAc                                               <  2:2809298/150‑1 (MQ=255)
                                                                                                     cGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAcc                                              <  1:4422579/150‑1 (MQ=255)
                                                                                                       cATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGa                                            <  2:2132106/150‑1 (MQ=255)
                                                                                                       cATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGa                                            >  2:4146664/1‑150 (MQ=255)
                                                                                                         tAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACACTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGc                                          >  2:3906384/1‑150 (MQ=255)
                                                                                                           aCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCAt                                        <  1:818915/150‑1 (MQ=255)
                                                                                                              ggTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGa                                     <  1:1609092/150‑1 (MQ=255)
                                                                                                                ttAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATg                                   <  1:1132678/150‑1 (MQ=255)
                                                                                                                   gTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGa                                <  1:4611097/150‑1 (MQ=255)
                                                                                                                    tACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGac                               <  2:1835464/150‑1 (MQ=255)
                                                                                                                     aCCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGaca                              >  2:4235589/1‑150 (MQ=255)
                                                                                                                      ccGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGacac                             >  2:1657239/1‑150 (MQ=255)
                                                                                                                      ccGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGacac                             <  1:3506469/150‑1 (MQ=255)
                                                                                                                      ccGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGacac                             <  1:92070/150‑1 (MQ=255)
                                                                                                                       cGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACAcc                            >  1:116889/1‑150 (MQ=255)
                                                                                                                       cGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACAcc                            >  1:4705092/1‑150 (MQ=255)
                                                                                                                        gTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCg                           <  1:587315/150‑1 (MQ=255)
                                                                                                                          tttGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAt                         >  1:4981569/1‑150 (MQ=255)
                                                                                                                              cTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTc                     <  2:4327844/150‑1 (MQ=255)
                                                                                                                               tGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCg                    <  2:857150/150‑1 (MQ=255)
                                                                                                                                    ccACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCg               <  1:3107226/150‑1 (MQ=255)
                                                                                                                                    ccACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCg               >  2:1994916/1‑150 (MQ=255)
                                                                                                                                    ccACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCg               <  2:3541786/150‑1 (MQ=255)
                                                                                                                                     cACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGt              <  1:3739559/150‑1 (MQ=255)
                                                                                                                                       cTGACAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTcc            <  1:3035809/150‑1 (MQ=255)
                                                                                                                                         gCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAg          <  2:4454427/150‑1 (MQ=255)
                                                                                                                                         gCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAg          >  2:4508836/1‑150 (MQ=255)
                                                                                                                                         gCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAg          <  1:1861303/150‑1 (MQ=255)
                                                                                                                                          ccAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGc         >  2:4940366/1‑150 (MQ=255)
                                                                                                                                           cAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGca        >  1:2122347/1‑150 (MQ=255)
                                                                                                                                            aGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcac       >  2:4776145/1‑150 (MQ=255)
                                                                                                                                            aGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcac       >  1:1238192/1‑150 (MQ=255)
                                                                                                                                               aaCCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcaccac    >  2:1517814/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACAT  >  CP009273/2322462‑2322756

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: