Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,405,662 C→A 10.4% I217I (ATC→ATA yhdJ → DNA adenine methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,405,6620CA10.4% 780.3 / 55.9 231I217I (ATC→ATAyhdJDNA adenine methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (112/95);  new base A (13/11);  total (125/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 235 total aligned reads displayed.

CAGCCTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCATTACTCTGCGGAAGAACT  >  CP009273/3405513‑3405811
                                                                                                                                                     |                                                                                                                                                     
caGCCTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATc                                                                                                                                                       >  1:2791548/1‑150 (MQ=255)
      tACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTtctc                                                                                                                                                 <  2:697507/150‑1 (MQ=255)
        cAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGa                                                                                                                                               >  1:828074/1‑150 (MQ=255)
         aaTCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAGGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGAc                                                                                                                                              >  1:2354970/1‑150 (MQ=255)
         aaTCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGAc                                                                                                                                              <  2:3147271/150‑1 (MQ=255)
         aaTCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGAc                                                                                                                                              >  1:3610073/1‑150 (MQ=255)
         aaTCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATAGTTCTCGAc                                                                                                                                              >  1:828823/1‑150 (MQ=255)
          atcatcAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGACAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTTTTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGAcc                                                                                                                                             <  1:3555434/150‑1 (MQ=255)
           tcatcaAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGAccc                                                                                                                                            <  1:589559/150‑1 (MQ=255)
                aaaaaGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTg                                                                                                                                       >  2:382031/1‑150 (MQ=255)
                   aaGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTg                                                                                                                                    >  1:3061542/1‑150 (MQ=255)
                   aaGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTg                                                                                                                                    >  2:538973/1‑150 (MQ=255)
                     gTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGt                                                                                                                                  >  1:1446073/1‑150 (MQ=255)
                     gTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGt                                                                                                                                  >  1:3033088/1‑150 (MQ=255)
                      tACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTa                                                                                                                                 >  1:3217179/1‑150 (MQ=255)
                         cGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATAGTTCTCGACCCGTTTGCTGGTAGCt                                                                                                                              >  1:4656570/1‑150 (MQ=255)
                            gTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTta                                                                                                                           >  2:2659298/1‑150 (MQ=255)
                               aCGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTtacta                                                                                                                        >  2:3409554/1‑150 (MQ=255)
                                    tGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGt                                                                                                                   >  2:4837886/1‑150 (MQ=255)
                                      ggATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGc                                                                                                                 >  2:116447/1‑150 (MQ=255)
                                         ttttCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGt                                                                                                              >  2:1987159/1‑150 (MQ=255)
                                         ttttCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGt                                                                                                              >  2:2535961/1‑150 (MQ=255)
                                          tttCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTa                                                                                                             <  1:5021421/150‑1 (MQ=255)
                                          tttCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTa                                                                                                             <  1:4313563/150‑1 (MQ=255)
                                             ccGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACTCAAAAACCGGAAGCCTTACTGAAACCCATTCTACTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTCCTACCGGTGCCGTAGcc                                                                                                          >  1:3971852/1‑150 (MQ=255)
                                             ccGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGcc                                                                                                          <  1:1105835/150‑1 (MQ=255)
                                             ccGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGcc                                                                                                          >  2:2326295/1‑150 (MQ=255)
                                             ccGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGcc                                                                                                          >  1:1028657/1‑150 (MQ=255)
                                             ccGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGcc                                                                                                          >  1:1219795/1‑150 (MQ=255)
                                              cgcgcgTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCa                                                                                                         <  1:685456/150‑1 (MQ=255)
                                                cgcgTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATc                                                                                                       >  1:2214732/1‑150 (MQ=255)
                                                cgcgTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATc                                                                                                       <  2:3425671/150‑1 (MQ=255)
                                                cgcgTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATc                                                                                                       >  2:3539563/1‑150 (MQ=255)
                                                cgcgTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATc                                                                                                       >  1:4799444/1‑150 (MQ=255)
                                                 gcgtgcgtTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATca                                                                                                      <  2:692095/150‑2 (MQ=255)
                                                        ttATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAg                                                                                                                                <  1:514995/117‑1 (MQ=255)
                                                          aTTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGttt                                                                                                                                        <  2:1668462/107‑1 (MQ=255)
                                                          aTTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGAc                                                                                             <  2:1358861/150‑1 (MQ=255)
                                                                 ggatgAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAAtt                                                                                      >  2:3034546/1‑150 (MQ=255)
                                                                 ggatgAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAAtt                                                                                      <  1:638939/150‑1 (MQ=255)
                                                                 ggatgAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAAtt                                                                                      >  2:5152695/1‑150 (MQ=255)
                                                                  gatgAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTc                                                                                     >  2:4333742/1‑150 (MQ=255)
                                                                    tgAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCAt                                                                                   <  2:2634075/150‑1 (MQ=255)
                                                                      aaTATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTg                                                                                 <  1:2330272/150‑1 (MQ=255)
                                                                      aaTATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTg                                                                                 <  1:3221542/150‑1 (MQ=255)
                                                                      aaTATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTAGTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTg                                                                                 <  1:3905740/150‑1 (MQ=255)
                                                                          tGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTAt                                                                             <  2:797299/150‑1 (MQ=255)
                                                                              aaCCACCCGACGCAAAAACCGGAAGCCTTACAGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTgag                                                                         >  1:2679025/1‑150 (MQ=255)
                                                                                ccaccCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGAt                                                                       >  2:4140086/1‑150 (MQ=255)
                                                                                 caccCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTTAGATc                                                                      >  2:4738062/1‑150 (MQ=255)
                                                                                  accCGACGCAAAAACCGGATGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCa                                                                     <  2:2354970/150‑1 (MQ=255)
                                                                                  accCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCa                                                                     <  1:1764671/150‑1 (MQ=255)
                                                                                    ccGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTCTTGAGATCAAc                                                                   >  1:43733/1‑150 (MQ=255)
                                                                                    ccGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAAc                                                                   <  1:2569477/150‑1 (MQ=255)
                                                                                     cGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACa                                                                  >  2:3886369/1‑150 (MQ=255)
                                                                                      gACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGACAATTCATTGGTATTGAGATCAACAg                                                                 >  2:5107683/1‑150 (MQ=255)
                                                                                      gACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAg                                                                 >  2:2712982/1‑150 (MQ=255)
                                                                                      gACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAg                                                                 >  1:4202660/1‑150 (MQ=255)
                                                                                      gACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAg                                                                 >  1:4712066/1‑150 (MQ=255)
                                                                                          cAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAg                                                             <  1:4580078/150‑1 (MQ=255)
                                                                                           aaaaaCCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGt                                                            >  1:1637115/1‑150 (MQ=255)
                                                                                           aaaaaCCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGt                                                            >  1:3424180/1‑150 (MQ=255)
                                                                                            aaaaCCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTa                                                           <  2:3824152/150‑1 (MQ=255)
                                                                                            aaaaCCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTa                                                           >  2:1531166/1‑150 (MQ=255)
                                                                                            aaaaCCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTa                                                           <  1:1306591/150‑1 (MQ=255)
                                                                                             aaaCCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTAc                                                          <  2:1359818/150‑1 (MQ=255)
                                                                                               aCCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACAt                                                        <  2:3767495/150‑1 (MQ=255)
                                                                                               aCCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACAt                                                        <  2:2262359/150‑1 (MQ=255)
                                                                                                  ggAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCaa                                                     >  1:4089407/1‑150 (MQ=255)
                                                                                                  ggAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCaa                                                     >  1:927970/1‑150 (MQ=255)
                                                                                                    aaGACTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCaaaa                                                   <  1:460412/150‑1 (MQ=255)
                                                                                                     aGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAAt                                                  <  1:3678431/150‑1 (MQ=255)
                                                                                                     aGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATAGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAAt                                                  >  2:4505097/1‑150 (MQ=255)
                                                                                                     aGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATAGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAAt                                                  >  2:1698210/1‑150 (MQ=255)
                                                                                                      gCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATg                                                 >  2:3977773/1‑150 (MQ=255)
                                                                                                        cTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATggg                                               >  2:2420287/1‑150 (MQ=255)
                                                                                                        cTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATggg                                               >  2:111527/1‑150 (MQ=255)
                                                                                                           aCTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCt                                            <  2:2990061/150‑1 (MQ=255)
                                                                                                           aCTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCt                                            <  2:83932/150‑1 (MQ=255)
                                                                                                               aaaCGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGa                                        <  2:3600508/150‑1 (MQ=255)
                                                                                                                    cATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCt                                   >  1:1935132/1‑150 (MQ=255)
                                                                                                                    cATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCt                                   <  1:965784/150‑1 (MQ=255)
                                                                                                                      ttattCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTgg                                 >  2:300927/1‑150 (MQ=255)
                                                                                                                      ttattCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTgg                                 >  1:2424745/1‑150 (MQ=255)
                                                                                                                       tattCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGa                                >  2:3712438/1‑150 (MQ=255)
                                                                                                                        attCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGAt                               <  2:2220469/150‑1 (MQ=255)
                                                                                                                          tctcGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATgt                             >  2:1852445/1‑150 (MQ=255)
                                                                                                                           ctcGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATgtc                            >  2:4161982/1‑150 (MQ=255)
                                                                                                                                 tctTCCAACCCTGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATgtcgcgtcg                      <  2:889054/150‑1 (MQ=255)
                                                                                                                                  ctTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATgtcgcgtcgc                     >  1:5056841/1‑150 (MQ=255)
                                                                                                                                   ttCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCa                    >  2:3745555/1‑150 (MQ=255)
                                                                                                                                    tCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCAt                   <  1:954366/150‑1 (MQ=255)
                                                                                                                                     ccAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCAtt                  >  2:4497081/1‑150 (MQ=255)
                                                                                                                                         cccAGGCGATATCGGTCTCGACCCGTTTGCTGGTATCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCATTActc              <  1:3767475/150‑1 (MQ=255)
                                                                                                                                               cGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCATTACTCTGCGga        >  1:4750026/1‑150 (MQ=255)
                                                                                                                                                gATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCATTACTCTGCGgaa       <  1:2081753/150‑1 (MQ=255)
                                                                                                                                                  tatCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCATTACTCTGCGgaaga     >  2:1578210/1‑150 (MQ=255)
                                                                                                                                                  tatCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCATTACTCTGCGgaaga     <  2:1438495/150‑1 (MQ=255)
                                                                                                                                                    tCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCATTACTCTGCGGAAGAAc   <  1:3409554/150‑1 (MQ=255)
                                                                                                                                                     cGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCATTACTCTGCGGAAGAACt  <  2:458702/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCCTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACGGCTGGATGTCGCGTCGCATTACTCTGCGGAAGAACT  >  CP009273/3405513‑3405811

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: