Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,508,852 G→T 22.8% P141T (CCC→ACC)  dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,508,8520GT22.8% 634.0 / 128.0 239P141T (CCC→ACC) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (87/97);  new base T (19/36);  total (106/133)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

GTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGT  >  CP009273/3508704‑3509001
                                                                                                                                                    |                                                                                                                                                     
gTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGt                                                                                                                                                      <  1:1170290/150‑1 (MQ=255)
gTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGACACAACTCTGCTTCCGGGAAATAGGGt                                                                                                                                                      <  1:4554404/150‑1 (MQ=255)
 tCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGtt                                                                                                                                                     <  2:4481805/150‑1 (MQ=255)
 tCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGtt                                                                                                                                                     <  2:3590783/150‑1 (MQ=255)
 tCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGtt                                                                                                                                                     <  2:184049/150‑1 (MQ=255)
  cGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGtttt                                                                                                                                                    >  1:1246126/1‑150 (MQ=255)
      gACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTg                                                                                                                                                <  1:2682566/150‑1 (MQ=255)
         aGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAg                                                                                                                                             <  2:2972359/150‑1 (MQ=255)
           cGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCg                                                                                                                                           <  2:746367/150‑1 (MQ=255)
              tGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGcc                                                                                                                                        <  1:2348155/150‑1 (MQ=255)
               gCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCg                                                                                                                                       <  2:4210357/150‑1 (MQ=255)
               gCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCg                                                                                                                                       <  2:4209862/150‑1 (MQ=255)
               gCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCg                                                                                                                                       <  1:2298620/150‑1 (MQ=255)
                 aTAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGaa                                                                                                                                     <  2:77270/150‑1 (MQ=255)
                  tAAGGCGGATCGCAATAGACGACGGATGCATCAACTGCGCGCGCCATGCTATAGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTATGTGTAGCGGcagtac                                                                                                                                    >  1:893642/1‑145 (MQ=255)
                   aaGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACg                                                                                                                                   <  1:3000214/150‑1 (MQ=255)
                   aaGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACg                                                                                                                                   >  2:3114409/1‑150 (MQ=255)
                        ggATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACg                                                                                                                              >  2:2021254/1‑150 (MQ=255)
                         gATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGt                                                                                                                             >  2:2242269/1‑150 (MQ=255)
                          aTCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGtt                                                                                                                            <  1:140544/150‑1 (MQ=255)
                              cAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAAc                                                                                                                        <  2:3912057/150‑1 (MQ=255)
                               aaTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACt                                                                                                                       <  1:3035355/150‑1 (MQ=255)
                               aaTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACt                                                                                                                       <  2:3649640/150‑1 (MQ=255)
                                aTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTc                                                                                                                      <  2:4230158/150‑1 (MQ=255)
                                 tAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCa                                                                                                                     <  1:1307089/150‑1 (MQ=255)
                                  aGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAc                                                                                                                    >  1:3285574/1‑150 (MQ=255)
                                  aGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAc                                                                                                                    >  1:1929754/1‑150 (MQ=255)
                                   gacgacGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAcc                                                                                                                   <  2:1015625/150‑1 (MQ=255)
                                     cgacgGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACcgc                                                                                                                 >  1:1647473/1‑150 (MQ=255)
                                      gacgGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACcgcg                                                                                                                >  1:3902439/1‑150 (MQ=255)
                                         ggATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGACACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAg                                                                                                             <  2:441035/150‑1 (MQ=255)
                                           aTGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGAt                                                                                                           <  1:2923267/150‑1 (MQ=255)
                                            tGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGAtt                                                                                                          <  1:77058/150‑1 (MQ=255)
                                            tGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGAtt                                                                                                          <  1:4407670/150‑1 (MQ=255)
                                             gCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTg                                                                                                         <  2:129782/150‑1 (MQ=255)
                                              catcatCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGt                                                                                                        >  2:2726509/1‑150 (MQ=255)
                                              catcatCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGt                                                                                                        >  2:637390/1‑150 (MQ=255)
                                               atcatcTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTa                                                                                                       <  2:3992324/150‑1 (MQ=255)
                                                 catcTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAAc                                                                                                     >  1:645115/1‑150 (MQ=255)
                                                   tcTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGa                                                                                                   <  1:1908152/150‑1 (MQ=255)
                                                     tGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGaca                                                                                                 >  2:4596477/1‑150 (MQ=255)
                                                          gcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGc                                                                                            <  1:2251490/150‑1 (MQ=255)
                                                          gcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGc                                                                                            <  1:3917254/150‑1 (MQ=255)
                                                          gcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGc                                                                                            <  2:4566324/150‑1 (MQ=255)
                                                            gcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCg                                                                                          <  2:2360381/150‑1 (MQ=255)
                                                             ccATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGt                                                                                         >  1:4191775/1‑150 (MQ=255)
                                                              cATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGtt                                                                                        <  1:3187610/150‑1 (MQ=255)
                                                                 gCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTa                                                                                     >  2:1528322/1‑150 (MQ=255)
                                                                         cGTAAGACTCACAATAGAAAAAGGCAGTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGc                                                                             <  1:3976316/150‑1 (MQ=255)
                                                                          gTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCg                                                                            >  1:4270638/1‑150 (MQ=255)
                                                                            aaGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGt                                                                          >  2:2803186/1‑150 (MQ=255)
                                                                             aGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGtc                                                                         <  1:3947212/150‑2 (MQ=255)
                                                                               aCTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCa                                                                       <  2:911322/150‑1 (MQ=255)
                                                                                cTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCaa                                                                      >  2:3606441/1‑150 (MQ=255)
                                                                                   acaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAata                                                                   <  1:3787698/150‑1 (MQ=255)
                                                                                   acaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAata                                                                   >  1:4267947/1‑150 (MQ=255)
                                                                                    caATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAatat                                                                  >  1:3031122/1‑150 (MQ=255)
                                                                                       tAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAATCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaa                                                               >  2:1207434/1‑150 (MQ=255)
                                                                                        aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaa                                                              >  1:2785330/1‑150 (MQ=255)
                                                                                        aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaa                                                              >  1:3375902/1‑150 (MQ=255)
                                                                                        aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaa                                                              <  1:202982/150‑1 (MQ=255)
                                                                                        aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaa                                                              >  1:138355/1‑150 (MQ=255)
                                                                                         gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaaa                                                             >  2:2295928/1‑150 (MQ=255)
                                                                                           aaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACa                                                           <  2:2452822/150‑1 (MQ=255)
                                                                                           aaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACa                                                           <  2:3188612/150‑1 (MQ=255)
                                                                                           aaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACa                                                           <  2:141910/150‑1 (MQ=255)
                                                                                                  aTTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGc                                                    <  1:2161261/150‑1 (MQ=255)
                                                                                                     cTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCg                                                 >  2:2030607/1‑150 (MQ=255)
                                                                                                      tGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGa                                                >  1:45943/1‑150 (MQ=255)
                                                                                                      tGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGa                                                >  1:4132743/1‑150 (MQ=255)
                                                                                                      tGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGa                                                >  1:3276361/1‑150 (MQ=255)
                                                                                                        cgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacg                                              <  2:2565186/150‑1 (MQ=255)
                                                                                                           tttttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacggaa                                           >  1:2038889/1‑150 (MQ=255)
                                                                                                            ttttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacggaac                                          >  2:3059232/1‑150 (MQ=255)
                                                                                                                cAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGAt                                      <  1:2300182/150‑1 (MQ=255)
                                                                                                                cAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGAGTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGAt                                      >  2:2921898/1‑150 (MQ=255)
                                                                                                                 aGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATc                                     >  2:1605037/1‑150 (MQ=255)
                                                                                                                   cGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCt                                   <  1:3521230/150‑1 (MQ=255)
                                                                                                                   cGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCt                                   <  2:115929/150‑1 (MQ=255)
                                                                                                                    gAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTg                                  <  1:3413480/150‑1 (MQ=255)
                                                                                                                          aTACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCtttt                            >  2:3163197/1‑150 (MQ=255)
                                                                                                                          aTACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCtttt                            <  2:3580548/150‑1 (MQ=255)
                                                                                                                           tACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttg                           >  1:2959612/1‑150 (MQ=255)
                                                                                                                           tACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttg                           >  1:3988097/1‑150 (MQ=255)
                                                                                                                               aCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGa                       <  2:4355594/150‑1 (MQ=255)
                                                                                                                                  ctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAct                    <  2:2656407/150‑1 (MQ=255)
                                                                                                                                   tGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGACCGGATCCTGGCTTTTGTTGAActc                   >  2:1716864/1‑150 (MQ=255)
                                                                                                                                     cTTCCGGGAAATCGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                 <  2:1957127/150‑1 (MQ=255)
                                                                                                                                      ttCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTc                <  2:3784796/150‑1 (MQ=255)
                                                                                                                                       tCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcg               <  1:3489904/150‑1 (MQ=255)
                                                                                                                                         cGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcgcg             >  1:250887/1‑150 (MQ=255)
                                                                                                                                          gggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCgg            >  1:3905665/1‑150 (MQ=255)
                                                                                                                                           ggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGa           >  1:1040462/1‑150 (MQ=255)
                                                                                                                                           ggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGa           >  2:332339/1‑150 (MQ=255)
                                                                                                                                            gAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGaa          >  2:3329759/1‑150 (MQ=255)
                                                                                                                                               aTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg       <  1:533608/150‑1 (MQ=255)
                                                                                                                                                   ggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAg   <  1:3755712/150‑1 (MQ=255)
                                                                                                                                                   ggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAg   <  2:943566/150‑1 (MQ=255)
                                                                                                                                                    gTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGt  <  1:1694179/150‑1 (MQ=255)
                                                                                                                                                    gTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGt  <  2:4027911/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGT  >  CP009273/3508704‑3509001

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