Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,508,852 G→T 5.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,8520GT 5.8% 868.6 / 20.5 231P141T (CCC→ACC) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (115/102);  new base T (10/4);  total (125/106)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

GGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAG  >  CP009273/3508703‑3509000
                                                                                                                                                     |                                                                                                                                                    
 gTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGt                                                                                                                                                     <  2:2579320/150‑1 (MQ=255)
     cAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGttttt                                                                                                                                                 <  2:612321/150‑1 (MQ=255)
      aGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGtttttt                                                                                                                                                >  2:3237057/1‑150 (MQ=255)
       gACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTg                                                                                                                                               <  1:4415167/150‑1 (MQ=255)
         cAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTa                                                                                                                                             >  1:1815744/1‑150 (MQ=255)
           gCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGc                                                                                                                                           <  1:3666810/150‑1 (MQ=255)
                 cATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGa                                                                                                                                     >  2:3699343/1‑150 (MQ=255)
                        cGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCAc                                                                                                                              >  1:1125239/1‑150 (MQ=255)
                         ggATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACg                                                                                                                             <  2:2984559/150‑1 (MQ=255)
                          gATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGt                                                                                                                            >  1:1426235/1‑150 (MQ=255)
                            tCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTa                                                                                                                          <  1:315282/150‑1 (MQ=255)
                              tcAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTaaa                                                                                                                        <  2:3877609/149‑1 (MQ=255)
                              gCAATAGACGACGTATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTaaa                                                                                                                        <  2:1799522/150‑1 (MQ=255)
                               cAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAAc                                                                                                                       >  2:2613956/1‑150 (MQ=255)
                                aaTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACt                                                                                                                      >  1:2146953/1‑150 (MQ=255)
                                aaTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACt                                                                                                                      >  2:3653705/1‑150 (MQ=255)
                                aaTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACt                                                                                                                      >  1:1792016/1‑150 (MQ=255)
                                  tAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCa                                                                                                                    >  2:598779/1‑150 (MQ=255)
                                   aGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTTGGACTCACTATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAc                                                                                                                   <  1:55939/150‑1 (MQ=255)
                                    gacgacGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAcc                                                                                                                  >  1:1142355/1‑150 (MQ=255)
                                    gacgacGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAcc                                                                                                                  <  1:3972019/150‑1 (MQ=255)
                                       gacgGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACcgcg                                                                                                               <  2:1541923/150‑1 (MQ=255)
                                       gacgGATGCATCAACTGCGCGCGCCATGCGATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTATTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACcgcg                                                                                                               <  2:455256/150‑1 (MQ=255)
                                          ggATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAg                                                                                                            >  2:4471905/1‑150 (MQ=255)
                                             tGCATCATCTGCGCGCGCCATTCTATCGGAGTAAGACTCACAATAGACAAAGTCATTCTGCGCTTTTTCATCGAATTTATACAACTCTGCTTCCTGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGAtt                                                                                                         <  2:3951331/150‑1 (MQ=255)
                                                       gcgcgcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGacac                                                                                               <  1:3501758/150‑1 (MQ=255)
                                                        cgcgcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTa                                                                                                                                             <  2:4379051/103‑1 (MQ=255)
                                                          cgcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAgg                                                                                            >  2:3984023/1‑150 (MQ=255)
                                                           gcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGc                                                                                           <  1:3731771/150‑1 (MQ=255)
                                                            cgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGcc                                                                                          >  1:855991/1‑150 (MQ=255)
                                                            cgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGcc                                                                                          <  2:3511860/150‑1 (MQ=255)
                                                                 tGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGt                                                                                     <  1:3266126/150‑1 (MQ=255)
                                                                 tGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGt                                                                                     >  1:4156983/1‑150 (MQ=255)
                                                                    tATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAAc                                                                                  <  2:865355/150‑1 (MQ=255)
                                                                    tATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAAGAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAAc                                                                                  >  1:4042365/1‑150 (MQ=255)
                                                                          cGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGACGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGc                                                                            <  2:1492119/150‑1 (MQ=255)
                                                                          cGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCAGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGc                                                                            <  2:1507645/150‑1 (MQ=255)
                                                                           gTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCg                                                                           >  2:611518/1‑150 (MQ=255)
                                                                           gTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCg                                                                           <  2:2398662/150‑1 (MQ=255)
                                                                           gTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCg                                                                           <  2:719953/150‑1 (MQ=255)
                                                                            tAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCgg                                                                          >  2:668712/1‑150 (MQ=255)
                                                                            tAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCgg                                                                          >  2:4258934/1‑150 (MQ=255)
                                                                              aGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGtt                                                                        >  1:1875215/1‑150 (MQ=255)
                                                                              aGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGtt                                                                        >  2:1062411/1‑150 (MQ=255)
                                                                               gACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTc                                                                       >  2:1381700/1‑150 (MQ=255)
                                                                                 cTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCaa                                                                     >  2:922195/1‑150 (MQ=255)
                                                                                  tCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCaaa                                                                    >  1:729769/1‑150 (MQ=255)
                                                                                   cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGACCGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAat                                                                   >  1:1804832/1‑150 (MQ=255)
                                                                                   cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAat                                                                   <  1:565415/150‑1 (MQ=255)
                                                                                   cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAat                                                                   <  1:3653705/150‑1 (MQ=255)
                                                                                    acaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAata                                                                  <  2:4125185/150‑1 (MQ=255)
                                                                                     caATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAatat                                                                 >  2:1331893/1‑150 (MQ=255)
                                                                                         aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaa                                                             >  2:2753219/1‑150 (MQ=255)
                                                                                          gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaaa                                                            <  1:4167656/150‑1 (MQ=255)
                                                                                          gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaaa                                                            >  1:252987/1‑150 (MQ=255)
                                                                                               aGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTa                                                       <  2:1526136/150‑1 (MQ=255)
                                                                                                ggCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTAc                                                      >  1:4513061/1‑150 (MQ=255)
                                                                                                  cATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCg                                                    >  2:4332597/1‑150 (MQ=255)
                                                                                                  cATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCg                                                    >  2:4330558/1‑150 (MQ=255)
                                                                                                  cATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCg                                                    >  2:3785178/1‑150 (MQ=255)
                                                                                                  cATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCg                                                    >  2:305765/1‑150 (MQ=255)
                                                                                                   aTTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGc                                                   >  1:3237765/1‑150 (MQ=255)
                                                                                                    ttCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGcc                                                  >  1:1483281/1‑150 (MQ=255)
                                                                                                      cTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTGTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCTCATATAAAACCAGTACCGCCCg                                                >  1:3924371/1‑150 (MQ=255)
                                                                                                        gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGac                                              <  1:1779800/150‑1 (MQ=255)
                                                                                                            tttttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacggaa                                          >  1:2538136/1‑150 (MQ=255)
                                                                                                            tttttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacggaa                                          >  1:4233014/1‑150 (MQ=255)
                                                                                                            tttttCAGCGAAGTGATACAACTCTGCGTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacggaa                                          >  1:224645/1‑150 (MQ=255)
                                                                                                              tttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacggaacg                                        <  1:3217971/150‑1 (MQ=255)
                                                                                                                  aGCGAAGTGATACACCTCTGCTTCCGGGATATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATc                                    <  2:3619122/150‑1 (MQ=255)
                                                                                                                  aGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATc                                    >  2:2890020/1‑150 (MQ=255)
                                                                                                                   gCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATcc                                   <  1:2357132/150‑1 (MQ=255)
                                                                                                                    cGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCt                                  >  1:3554232/1‑150 (MQ=255)
                                                                                                                    cGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCt                                  >  1:595974/1‑150 (MQ=255)
                                                                                                                    cGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCt                                  >  2:539595/1‑150 (MQ=255)
                                                                                                                      aaGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTgg                                <  2:2765839/150‑1 (MQ=255)
                                                                                                                         tGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCtt                             <  2:1472261/150‑1 (MQ=255)
                                                                                                                          gATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCttt                            >  1:287582/1‑150 (MQ=255)
                                                                                                                           aTACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCtttt                           <  1:2944726/150‑1 (MQ=255)
                                                                                                                             aCAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgt                         >  2:2005717/1‑150 (MQ=255)
                                                                                                                                  tctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAc                    <  1:1583352/150‑1 (MQ=255)
                                                                                                                                   ctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGAAGGAACGGATCCTGGCTTTTGTTGAAct                   >  2:3674321/1‑150 (MQ=255)
                                                                                                                                    tGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAActc                  <  2:2502720/150‑1 (MQ=255)
                                                                                                                                     gCTTCCGTGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACTTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAACTATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAActct                 <  2:2993718/150‑1 (MQ=255)
                                                                                                                                     gCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAActct                 <  1:2416923/150‑1 (MQ=255)
                                                                                                                                      cTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                <  2:4058717/150‑1 (MQ=255)
                                                                                                                                      cTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                <  2:2530780/150‑1 (MQ=255)
                                                                                                                                      cTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                >  2:1669272/1‑150 (MQ=255)
                                                                                                                                       ttCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTc               >  1:3243744/1‑150 (MQ=255)
                                                                                                                                       ttCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTc               <  2:1125239/150‑1 (MQ=255)
                                                                                                                                          cGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcgcg            >  1:1334926/1‑150 (MQ=255)
                                                                                                                                            ggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGa          >  2:4065102/1‑150 (MQ=255)
                                                                                                                                            ggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGa          >  1:1459924/1‑150 (MQ=255)
                                                                                                                                              aaaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAAc        >  1:378601/1‑150 (MQ=255)
                                                                                                                                              aaaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAAc        <  2:1733665/150‑1 (MQ=255)
                                                                                                                                               aaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACt       >  1:1171060/1‑150 (MQ=255)
                                                                                                                                                 tAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGa     <  2:2131617/150‑1 (MQ=255)
                                                                                                                                                   gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATa   <  1:2404234/150‑1 (MQ=255)
                                                                                                                                                    ggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAg  >  1:429837/1‑150 (MQ=255)
                                                                                                                                                    ggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAg  >  1:429998/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAG  >  CP009273/3508703‑3509000

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: