Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,322,609 C→A 81.7% 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,6090CA81.7% 535.4 / 152.5 262W223L (TGG→TTG) yfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (31/17);  new base A (124/90);  total (155/107)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 266 total aligned reads displayed.

GCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCAC  >  CP009273/2322462‑2322754
                                                                                                                                                   |                                                                                                                                                 
gCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCca                                                                                                                                                  >  1:1638931/1‑149 (MQ=255)
 cTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCaac                                                                                                                                                 >  1:1425385/1‑149 (MQ=255)
     gATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCg                                                                                                                                             >  1:3639622/1‑149 (MQ=255)
         ttcttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtgt                                                                                                                                         >  2:4126766/1‑149 (MQ=255)
         ttcttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtgt                                                                                                                                         >  2:2817823/1‑149 (MQ=255)
          tcttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTa                                                                                                                                        >  1:4113128/1‑149 (MQ=255)
           cttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTAc                                                                                                                                       <  1:1008372/149‑1 (MQ=255)
             tcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCa                                                                                                                                     <  2:4379443/149‑1 (MQ=255)
                gCGGATAGCGGGCAAACAGCTCCACCTCTACACAGTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCg                                                                                                                                  <  2:1349292/149‑1 (MQ=255)
                 cGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCgg                                                                                                                                 <  1:3957007/149‑1 (MQ=255)
                  ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGt                                                                                                                                <  1:1023233/149‑1 (MQ=255)
                  ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGt                                                                                                                                >  2:3478809/1‑149 (MQ=255)
                   gATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGtt                                                                                                                               <  2:2941870/149‑1 (MQ=255)
                   gATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGtt                                                                                                                               >  1:929004/1‑149 (MQ=255)
                    aTAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTa                                                                                                                              >  1:2733112/1‑149 (MQ=255)
                    aTAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTa                                                                                                                              <  1:3404046/149‑1 (MQ=255)
                     tAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTaa                                                                                                                             <  1:3548762/149‑1 (MQ=255)
                           gCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagcc                                                                                                                       <  2:1724087/149‑1 (MQ=255)
                            cAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAAccagcca                                                                                                                      <  2:1302615/149‑1 (MQ=255)
                            cAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagcca                                                                                                                      <  1:4225417/149‑1 (MQ=255)
                            cAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagcca                                                                                                                      <  2:2751031/149‑1 (MQ=255)
                             aaaCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagccag                                                                                                                     >  1:2427161/1‑149 (MQ=255)
                              aaCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACcagccagc                                                                                                                    <  1:2340338/149‑1 (MQ=255)
                               aCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCa                                                                                                                   >  2:4521448/1‑149 (MQ=255)
                                 aGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaaa                                                                                                                 >  1:1981066/1‑149 (MQ=255)
                                  gCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaaaa                                                                                                                >  1:2918726/1‑149 (MQ=255)
                                   cTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAac                                                                                                               <  2:3207195/149‑1 (MQ=255)
                                      cACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccac                                                                                                            >  2:4011159/1‑149 (MQ=255)
                                       aCCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccacc                                                                                                           >  2:1787762/1‑149 (MQ=255)
                                        ccTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCAccc                                                                                                          >  1:247005/1‑149 (MQ=255)
                                          tctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTg                                                                                                        <  2:3373231/149‑1 (MQ=255)
                                           ctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGa                                                                                                       <  1:4065351/149‑1 (MQ=255)
                                           ctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGa                                                                                                       <  2:1035246/149‑1 (MQ=255)
                                            tACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGAc                                                                                                      <  2:4213280/149‑1 (MQ=255)
                                            tACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGAc                                                                                                      >  1:305831/1‑149 (MQ=255)
                                             acacATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACg                                                                                                     >  1:3155411/1‑149 (MQ=255)
                                             acacATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACg                                                                                                     >  1:1859708/1‑149 (MQ=255)
                                              cacaTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACgg                                                                                                    <  1:1388851/149‑1 (MQ=255)
                                                  ttGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCACAACCACCCTGACGGAttt                                                                                                >  2:3250081/1‑149 (MQ=255)
                                                  ttGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAttt                                                                                                >  1:3383402/1‑149 (MQ=255)
                                                  ttGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAttt                                                                                                <  2:2536957/149‑1 (MQ=255)
                                                      ccACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTg                                                                                            >  1:3963854/1‑149 (MQ=255)
                                                         cTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGttt                                                                                         >  2:4716203/1‑149 (MQ=255)
                                                          tGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTa                                                                                        >  2:316198/1‑149 (MQ=255)
                                                           ggCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTat                                                                                       <  1:4503043/149‑1 (MQ=255)
                                                               aCCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                   <  2:332904/149‑1 (MQ=255)
                                                                   gCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCa                                                                               <  2:3658690/149‑1 (MQ=255)
                                                                     gTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCaaa                                                                             >  1:845036/1‑149 (MQ=255)
                                                                       cGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAAcc                                                                           <  1:3122572/149‑1 (MQ=255)
                                                                             tGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACCGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTcc                                                                     <  1:704208/149‑1 (MQ=255)
                                                                             tGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTcc                                                                     >  2:4277792/1‑149 (MQ=255)
                                                                             tGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTcc                                                                     >  2:1914338/1‑149 (MQ=255)
                                                                              ggCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                    >  1:4155197/1‑149 (MQ=255)
                                                                              ggCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                    <  1:2028808/149‑1 (MQ=255)
                                                                              ggCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                    >  2:2645170/1‑149 (MQ=255)
                                                                               gCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAg                                                                   <  1:3620172/149‑1 (MQ=255)
                                                                                cAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAgg                                                                  >  1:2299934/1‑149 (MQ=255)
                                                                                 aaaaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGt                                                                 >  2:1637006/1‑149 (MQ=255)
                                                                                  aaaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtt                                                                <  2:3908857/149‑1 (MQ=255)
                                                                                  aaaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtt                                                                >  2:1933279/1‑149 (MQ=255)
                                                                                        gATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAcc                                                          >  1:330550/1‑149 (MQ=255)
                                                                                         aTTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAccc                                                         >  1:2697323/1‑149 (MQ=255)
                                                                                                  cGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGAc                                                <  2:4129886/149‑1 (MQ=255)
                                                                                                     cGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAc                                             >  2:4043981/1‑149 (MQ=255)
                                                                                                     cGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAc                                             <  2:1205468/149‑1 (MQ=255)
                                                                                                     cGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAc                                             <  2:278462/149‑1 (MQ=255)
                                                                                                       cATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCg                                           >  1:3826735/1‑149 (MQ=255)
                                                                                                       cATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCg                                           >  1:3792591/1‑149 (MQ=255)
                                                                                                         tAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                         >  1:182550/1‑149 (MQ=255)
                                                                                                         tAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                         >  1:3412419/1‑149 (MQ=255)
                                                                                                          aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGc                                        >  2:2308307/1‑149 (MQ=255)
                                                                                                          aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGc                                        >  2:3821765/1‑149 (MQ=255)
                                                                                                          aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGc                                        >  2:1948363/1‑149 (MQ=255)
                                                                                                          aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGc                                        <  2:4620244/149‑1 (MQ=255)
                                                                                                           aCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCa                                       >  1:4585309/1‑149 (MQ=255)
                                                                                                            ccggTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCAt                                      >  1:2793544/3‑149 (MQ=255)
                                                                                                             gggTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATg                                     >  2:4514087/1‑149 (MQ=255)
                                                                                                                ttAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGAt                                  <  2:1448220/149‑1 (MQ=255)
                                                                                                                ttAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGAt                                  >  1:4364800/1‑149 (MQ=255)
                                                                                                                  aGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGc                                >  1:3276997/1‑149 (MQ=255)
                                                                                                                  aGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGc                                <  2:722753/149‑1 (MQ=255)
                                                                                                                   gTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCg                               <  1:2533290/149‑1 (MQ=255)
                                                                                                                    tACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGa                              <  2:82718/149‑1 (MQ=255)
                                                                                                                    tACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGa                              >  1:1553851/1‑149 (MQ=255)
                                                                                                                     aCCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGac                             >  1:4430442/1‑149 (MQ=255)
                                                                                                                       cGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGacac                           >  2:1628985/1‑149 (MQ=255)
                                                                                                                          tttGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGa                        <  2:1166055/149‑1 (MQ=255)
                                                                                                                            tGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAtt                      >  1:1602609/1‑149 (MQ=255)
                                                                                                                               tGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTc                   >  2:2094019/1‑149 (MQ=255)
                                                                                                                                ggTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCg                  >  2:4034228/1‑149 (MQ=255)
                                                                                                                                   gCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCt               <  2:248196/149‑1 (MQ=255)
                                                                                                                                     cACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCg             >  2:296277/1‑149 (MQ=255)
                                                                                                                                       cTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTc           <  2:1829579/149‑1 (MQ=255)
                                                                                                                                          ccAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAg        >  2:941665/1‑149 (MQ=255)
                                                                                                                                          ccAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAg        >  2:2202086/1‑149 (MQ=255)
                                                                                                                                             gTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcac     >  1:4322579/1‑149 (MQ=255)
                                                                                                                                              tAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcacc    <  2:325832/149‑1 (MQ=255)
                                                                                                                                                aCCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcaccac  <  2:2434667/149‑1 (MQ=255)
                                                                                                                                                aCCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcaccac  >  1:526971/1‑149 (MQ=255)
                                                                                                                                                aCCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcaccac  <  1:3853829/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCAC  >  CP009273/2322462‑2322754

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: