Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,508,852 G→T 12.5% 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,8520GT12.5% 879.7 / 66.8 252P141T (CCC→ACC) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (111/109);  new base T (15/17);  total (126/126)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

GGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGT  >  CP009273/3508703‑3509001
                                                                                                                                                     |                                                                                                                                                     
ggTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAggg                                                                                                                                                       >  2:2293626/1‑150 (MQ=255)
 gTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGt                                                                                                                                                      <  1:418807/150‑1 (MQ=255)
     cAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGAGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGttttt                                                                                                                                                  >  1:1908894/1‑150 (MQ=255)
       gACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTg                                                                                                                                                >  2:395846/1‑150 (MQ=255)
        aCAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGt                                                                                                                                               <  2:3457145/150‑1 (MQ=255)
         cAGCGGTGCATCAGGCGGATCTCAATAGACGACGGATGCATAATCTGCGAGCGACCTGCTATAGGCGTAATACTCACAATAGAAAATGGCATTATGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTTTa                                                                                                                                              <  1:38914/150‑1 (MQ=255)
              gTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAACTAGGGTTTTTTGTAGCGGc                                                                                                                                         >  1:767030/1‑150 (MQ=255)
               tGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGcc                                                                                                                                        >  2:2951405/1‑150 (MQ=255)
                gCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCg                                                                                                                                       >  2:2826866/1‑150 (MQ=255)
                   tAAGGCGGATCGCAATAGTCGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTATGACTCACTATTGAAAAAGGCTTTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGTAc                                                                                                                                    <  1:4579895/150‑1 (MQ=255)
                   tAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAAc                                                                                                                                    <  2:47110/150‑1 (MQ=255)
                   tAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAAc                                                                                                                                    >  1:1486648/1‑150 (MQ=255)
                     aGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCAGGAAATAGGGTTTTTTGTAGCGGCCGAACgg                                                                                                                                  <  1:4134721/150‑1 (MQ=255)
                        cGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCAc                                                                                                                               <  1:531065/150‑1 (MQ=255)
                         ggATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACg                                                                                                                              >  1:2348092/1‑150 (MQ=255)
                         ggATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACg                                                                                                                              <  2:1561220/150‑1 (MQ=255)
                            tCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTa                                                                                                                           >  1:2135222/1‑150 (MQ=255)
                              gCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTaaa                                                                                                                         <  1:626214/150‑1 (MQ=255)
                                aaTAGACGACGGTTGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACt                                                                                                                       <  2:4477551/150‑1 (MQ=255)
                                   aGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAc                                                                                                                    <  1:25239/150‑1 (MQ=255)
                                    gacgacGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAcc                                                                                                                   <  2:300424/150‑1 (MQ=255)
                                    gacgacGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAcc                                                                                                                   <  2:290074/150‑1 (MQ=255)
                                     acgacgGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACcg                                                                                                                  <  2:1732179/150‑1 (MQ=255)
                                         cgGATGCATCATCTGCGCTCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACACTGCTACCGGGAAATATGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCa                                                                                                              >  1:2113372/1‑150 (MQ=255)
                                         cgGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCa                                                                                                              >  1:1341876/1‑150 (MQ=255)
                                          ggATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAg                                                                                                             >  2:3441422/1‑150 (MQ=255)
                                           gATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGa                                                                                                            <  1:1997361/150‑1 (MQ=255)
                                             tGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGAtt                                                                                                          >  1:3168191/1‑150 (MQ=255)
                                                atcatcTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGACATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTa                                                                                                       >  1:688677/1‑150 (MQ=255)
                                                atcatcTGCGCGCGCCATGCTATCGGCGTAAGACTAACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTa                                                                                                       <  1:121576/150‑1 (MQ=255)
                                                   atcTGCGCGCGCCATGCTATCTGCGTAAGACTCACAATAGAAAAATGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACg                                                                                                    <  1:4639344/150‑1 (MQ=255)
                                                   atcTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACg                                                                                                    <  1:4115312/150‑1 (MQ=255)
                                                      tGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGaca                                                                                                 <  1:461404/150‑1 (MQ=255)
                                                        cgcgcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGAcaca                                                                                               <  2:3113365/150‑1 (MQ=255)
                                                        cgcgcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGAcaca                                                                                               <  1:3039893/150‑1 (MQ=255)
                                                        cgcgcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGAcaca                                                                                               >  2:2046954/1‑150 (MQ=255)
                                                        cgcgcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGAcaca                                                                                               <  2:3166609/150‑1 (MQ=255)
                                                        cgcgcgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGAcaca                                                                                               <  1:935257/150‑1 (MQ=255)
                                                            cgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGcc                                                                                           >  1:4433377/1‑150 (MQ=255)
                                                             gcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCg                                                                                          <  1:102187/150‑1 (MQ=255)
                                                               cATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGtt                                                                                        <  1:3300189/150‑1 (MQ=255)
                                                                 tGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGt                                                                                      >  2:2008113/1‑150 (MQ=255)
                                                                  gCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTa                                                                                     >  1:1776602/1‑150 (MQ=255)
                                                                  gCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTa                                                                                     <  1:1191943/150‑1 (MQ=255)
                                                                  gCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTa                                                                                     <  1:2583747/150‑1 (MQ=255)
                                                                  gCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTa                                                                                     >  1:3794227/1‑150 (MQ=255)
                                                                  gCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTa                                                                                     >  1:3793990/1‑150 (MQ=255)
                                                                    tATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAAc                                                                                   >  1:4534103/1‑150 (MQ=255)
                                                                      tCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTATTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCg                                                                                 <  1:3368956/150‑1 (MQ=255)
                                                                        ggCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTa                                                                               <  1:3461947/150‑1 (MQ=255)
                                                                          cGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGc                                                                             <  2:2111444/150‑1 (MQ=255)
                                                                           gTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCg                                                                            <  1:2819968/150‑1 (MQ=255)
                                                                            tAAGACTCACAATAGAAAAAGGCATTCTGCTCTTTTTCAGCGAAGTGATACAACTCTGCTTACGTGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGAATTTAACGAAACAGGACGTTGAAACAGTAGCgg                                                                           >  2:2155456/1‑150 (MQ=255)
                                                                            tAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCgg                                                                           >  1:4506074/1‑150 (MQ=255)
                                                                              aGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGtt                                                                         >  2:271889/1‑150 (MQ=255)
                                                                               gACTCACAATAGACAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTc                                                                        <  2:4221876/150‑1 (MQ=255)
                                                                               gACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCGGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTc                                                                        <  2:3727751/150‑1 (MQ=255)
                                                                                aCTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCa                                                                       <  2:1272592/150‑1 (MQ=255)
                                                                                 cTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCaa                                                                      <  2:1298221/150‑1 (MQ=255)
                                                                                  tCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCaaa                                                                     >  1:2044052/1‑150 (MQ=255)
                                                                                  tCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCaaa                                                                     >  2:2665457/1‑150 (MQ=255)
                                                                                   cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAat                                                                    >  2:1415789/1‑150 (MQ=255)
                                                                                         aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaa                                                              >  2:830440/1‑150 (MQ=255)
                                                                                          gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaaa                                                             <  1:2150187/150‑1 (MQ=255)
                                                                                           aaaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAAc                                                            <  2:4156822/150‑1 (MQ=255)
                                                                                            aaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACa                                                           >  2:1323667/1‑150 (MQ=255)
                                                                                             aaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAg                                                          >  1:2731254/1‑150 (MQ=255)
                                                                                             aaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAg                                                          <  2:183930/150‑1 (MQ=255)
                                                                                               aGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTa                                                        >  2:3824399/1‑150 (MQ=255)
                                                                                                ggCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTAc                                                       >  1:4348713/1‑150 (MQ=255)
                                                                                                 gCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTAcc                                                      <  2:3609239/150‑1 (MQ=255)
                                                                                                     tCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGccc                                                  >  2:2176212/1‑150 (MQ=255)
                                                                                                     tCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGccc                                                  <  1:1560069/150‑1 (MQ=255)
                                                                                                      cTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCg                                                 >  2:1174284/1‑150 (MQ=255)
                                                                                                       tGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGa                                                >  2:2324064/1‑150 (MQ=255)
                                                                                                        gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGac                                               >  1:2900363/1‑150 (MQ=255)
                                                                                                           cTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGNTCAAATATAAAAACAGTACCGCCCGacgga                                            <  1:4279463/150‑1 (MQ=255)
                                                                                                                  aGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATc                                     <  1:1529734/150‑1 (MQ=255)
                                                                                                                     gAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCACGACGGAACGGATCCTg                                  <  1:3930690/150‑1 (MQ=255)
                                                                                                                        gTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCt                               <  2:3866938/150‑1 (MQ=255)
                                                                                                                         tGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCtt                              >  2:4390226/1‑150 (MQ=255)
                                                                                                                         tGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCtt                              >  2:4390070/1‑150 (MQ=255)
                                                                                                                            tACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttg                           <  2:325227/150‑1 (MQ=255)
                                                                                                                              cAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgtt                         >  1:2319933/1‑150 (MQ=255)
                                                                                                                               aaCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgttg                        <  2:469381/150‑1 (MQ=255)
                                                                                                                                 ctctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGaa                      >  1:1013081/1‑150 (MQ=255)
                                                                                                                                  tctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCCCCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAc                     <  2:1403308/150‑1 (MQ=255)
                                                                                                                                   ctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAct                    >  2:338730/1‑150 (MQ=255)
                                                                                                                                   ctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAct                    <  2:2519056/150‑1 (MQ=255)
                                                                                                                                      cTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                 <  1:229747/150‑1 (MQ=255)
                                                                                                                                      cTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                 >  2:3903688/1‑150 (MQ=255)
                                                                                                                                      cTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                 <  1:3255448/150‑1 (MQ=255)
                                                                                                                                         ccGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcgc              >  1:4024460/1‑150 (MQ=255)
                                                                                                                                           gggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCgg            >  1:240796/1‑150 (MQ=255)
                                                                                                                                              aaaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAAc         >  2:4074564/1‑150 (MQ=255)
                                                                                                                                              aaaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAAc         >  2:4074695/1‑150 (MQ=255)
                                                                                                                                                aTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg       >  1:3571004/1‑150 (MQ=255)
                                                                                                                                                aTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg       >  1:1869593/1‑150 (MQ=255)
                                                                                                                                                aTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg       >  1:426420/1‑150 (MQ=255)
                                                                                                                                                     gTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGt  <  2:2197895/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGT  >  CP009273/3508703‑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: