Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,688,591 A→C 5.7% S2R (AGT→AGG bcsA ← cellulose synthase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,688,5910AC 5.7% 789.0 / 33.6 257S2R (AGT→AGGbcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (116/126);  new base C (6/9);  total (122/135)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

TGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCAGCCATCGCTTCA  >  CP009273/3688442‑3688740
                                                                                                                                                     |                                                                                                                                                     
tGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATa                                                                                                                                                       <  1:380146/150‑1 (MQ=255)
  cTCCAGCGGAATAAAAGTCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACt                                                                                                                                                     >  1:2466447/1‑150 (MQ=255)
  cTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACt                                                                                                                                                     >  1:1950662/1‑150 (MQ=255)
  cTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACt                                                                                                                                                     >  1:1951418/1‑150 (MQ=255)
    ccAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCa                                                                                                                                                   >  1:4013520/1‑150 (MQ=255)
    ccAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCa                                                                                                                                                   >  1:134033/1‑150 (MQ=255)
    ccAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCa                                                                                                                                                   >  1:767601/1‑150 (MQ=255)
     cAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCat                                                                                                                                                  <  2:1079576/150‑1 (MQ=255)
       gCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCatga                                                                                                                                                <  2:1741196/150‑1 (MQ=255)
        cGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCatgat                                                                                                                                               <  1:3463509/150‑1 (MQ=255)
        cGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCatgat                                                                                                                                               >  2:2074202/1‑150 (MQ=255)
        cGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCatgat                                                                                                                                               <  1:451374/150‑1 (MQ=255)
        cGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCatgat                                                                                                                                               >  1:3725898/1‑150 (MQ=255)
         ggAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCatgatg                                                                                                                                              <  2:892669/150‑1 (MQ=255)
            aTAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCAc                                                                                                                                           <  2:1980616/150‑1 (MQ=255)
             tAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACt                                                                                                                                          <  1:4079031/150‑1 (MQ=255)
             tAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACt                                                                                                                                          <  1:3085580/150‑1 (MQ=255)
             tAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACt                                                                                                                                          >  2:4219258/1‑150 (MQ=255)
             tAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACt                                                                                                                                          <  1:1263909/150‑1 (MQ=255)
               aaaaTCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTcc                                                                                                                                        <  2:3229710/150‑1 (MQ=255)
                  aTCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGa                                                                                                                                     >  2:3497307/1‑150 (MQ=255)
                        gccaggATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCg                                                                                                                               <  1:2912703/150‑1 (MQ=255)
                        gccaggATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCg                                                                                                                               >  1:2795759/1‑150 (MQ=255)
                        gccaggATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCg                                                                                                                               <  1:146582/150‑1 (MQ=255)
                         ccaggATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGt                                                                                                                              >  1:3669449/1‑150 (MQ=255)
                           aggATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGttt                                                                                                                            <  1:2049133/150‑1 (MQ=255)
                            ggATCATCCAGAAACAGTCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGtttt                                                                                                                           >  2:4141920/1‑150 (MQ=255)
                            ggATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGtttt                                                                                                                           >  2:4141879/1‑150 (MQ=255)
                              atcatcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCa                                                                                                                         <  1:2264188/150‑1 (MQ=255)
                               tcatcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAg                                                                                                                        >  1:752749/1‑150 (MQ=255)
                                  tcCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGccc                                                                                                                     <  1:3366647/150‑1 (MQ=255)
                                     aGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGa                                                                                                                  >  1:1541003/1‑150 (MQ=255)
                                         aCAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGagtag                                                                                                              <  2:3790364/150‑1 (MQ=255)
                                            gCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATCCTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTc                                                                                                           <  1:2385575/150‑1 (MQ=255)
                                            gCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTc                                                                                                           >  2:923003/1‑150 (MQ=255)
                                            gCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAAAGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTc                                                                                                           <  2:4091295/150‑1 (MQ=255)
                                            gCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAAAGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTc                                                                                                           <  2:365365/150‑1 (MQ=255)
                                               gAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAAc                                                                                                        <  2:3719378/150‑1 (MQ=255)
                                                aGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACa                                                                                                       <  1:1413672/150‑1 (MQ=255)
                                                aGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACa                                                                                                       <  1:1414616/150‑1 (MQ=255)
                                                   gTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGc                                                                                                    <  1:3797153/150‑1 (MQ=255)
                                                     cgcgCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCAc                                                                                                  >  1:262769/1‑150 (MQ=255)
                                                     cgcgCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCAc                                                                                                  <  2:3268890/150‑1 (MQ=255)
                                                      gcgcTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCAcc                                                                                                 <  1:3001784/150‑1 (MQ=255)
                                                      gcgcTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCAcc                                                                                                 >  1:3984300/1‑150 (MQ=255)
                                                         cTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGt                                                                                              >  2:3316867/1‑150 (MQ=255)
                                                                ccGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAg                                                                                       >  2:4413041/1‑150 (MQ=255)
                                                                  gACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCg                                                                                     >  2:4104259/1‑150 (MQ=255)
                                                                   aCGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGt                                                                                    >  2:1170890/1‑150 (MQ=255)
                                                                    cGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGGTAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTc                                                                                   <  1:2186328/150‑1 (MQ=255)
                                                                    cGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTc                                                                                   <  1:2186340/150‑1 (MQ=255)
                                                                    cGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTc                                                                                   <  1:1627304/150‑1 (MQ=255)
                                                                    cGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTc                                                                                   <  2:298453/150‑1 (MQ=255)
                                                                    cGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTc                                                                                   <  2:1446306/150‑1 (MQ=255)
                                                                     gCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCa                                                                                  >  1:4285848/1‑150 (MQ=255)
                                                                     gCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCa                                                                                  >  1:2040202/1‑150 (MQ=255)
                                                                         cGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTAt                                                                              <  1:1075619/150‑1 (MQ=255)
                                                                         cGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTAt                                                                              >  2:1581905/1‑150 (MQ=255)
                                                                          gTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTAtc                                                                             <  2:2964195/150‑1 (MQ=255)
                                                                          gTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTAtc                                                                             <  1:3572729/150‑1 (MQ=255)
                                                                          gTGACGACGAGAATCGCGATCACGCCCGAGACGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTAtc                                                                             <  2:2964304/150‑1 (MQ=255)
                                                                           tGACGACGATAATCGCGATAACGCCCGATACGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTAtct                                                                            <  2:278574/150‑1 (MQ=255)
                                                                              cgacgaTAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCtt                                                                         <  2:229175/150‑1 (MQ=255)
                                                                                acgaTAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCa                                                                       <  2:549910/150‑1 (MQ=255)
                                                                                   aTAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAgcc                                                                    >  1:714468/1‑150 (MQ=255)
                                                                                   aTAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAgcc                                                                    >  1:2982935/1‑150 (MQ=255)
                                                                                      aTCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAgccgcc                                                                 <  2:3990754/150‑1 (MQ=255)
                                                                                      aTCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAgccgcc                                                                 <  2:3968503/150‑1 (MQ=255)
                                                                                      aTCGCGAGAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAgccgcc                                                                 <  2:2148578/150‑1 (MQ=255)
                                                                                       tCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCa                                                                <  1:3168251/150‑1 (MQ=255)
                                                                                        cgcgATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAg                                                               >  1:2098111/1‑150 (MQ=255)
                                                                                          cgATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAgcg                                                             <  2:2720433/150‑1 (MQ=255)
                                                                                          cgATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAgcg                                                             <  2:2720898/150‑1 (MQ=255)
                                                                                           gATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAgcgc                                                            >  2:2609794/1‑150 (MQ=255)
                                                                                           gATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAgcgc                                                            >  2:314648/1‑150 (MQ=255)
                                                                                                cGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCAc                                                       >  1:3112610/1‑150 (MQ=255)
                                                                                                 gCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATCCTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACt                                                      <  1:3294900/150‑1 (MQ=255)
                                                                                                 gCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACt                                                      >  2:3334992/1‑150 (MQ=255)
                                                                                                  cccGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCGTCAGCCGCCAGCGCATCACTg                                                     <  1:979506/150‑1 (MQ=255)
                                                                                                      aTAAGCCGCGCGTTGACCGTCGGGATAAGCAACCACCGGGTCAGGATCCTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGAATCACTGCGat                                                 <  1:2590800/150‑1 (MQ=255)
                                                                                                                cgTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTAc                                       >  2:1839905/1‑150 (MQ=255)
                                                                                                                  ttGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTg                                     <  1:1078659/150‑1 (MQ=255)
                                                                                                                  ttGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTg                                     <  2:4426447/150‑1 (MQ=255)
                                                                                                                    gACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGt                                   <  1:4238486/150‑1 (MQ=255)
                                                                                                                    gACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGt                                   <  1:1760147/150‑1 (MQ=255)
                                                                                                                      ccGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTg                                 <  1:2596667/150‑1 (MQ=255)
                                                                                                                       cggcggGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGc                                <  1:4413041/150‑1 (MQ=255)
                                                                                                                        ggcggGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCt                               <  1:1181069/150‑1 (MQ=255)
                                                                                                                          cggGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTa                             >  2:601140/1‑150 (MQ=255)
                                                                                                                          cggGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTa                             >  2:863943/1‑150 (MQ=255)
                                                                                                                           gggATAAGCAACCACCGGGTCAGGATCCTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAg                            >  2:2406156/1‑150 (MQ=255)
                                                                                                                           gggATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAg                            >  2:2985983/1‑150 (MQ=255)
                                                                                                                           gggATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAg                            <  2:752749/150‑1 (MQ=255)
                                                                                                                                 aGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCa                      >  1:4025844/1‑150 (MQ=255)
                                                                                                                                      ccaccGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCAt                 >  2:4064946/1‑150 (MQ=255)
                                                                                                                                         ccGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCa              <  1:1351861/150‑1 (MQ=255)
                                                                                                                                          cGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCAg             >  1:2738417/1‑150 (MQ=255)
                                                                                                                                             gTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCAGCCa          >  2:3023000/1‑150 (MQ=255)
                                                                                                                                               cAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCAGCCATc        >  2:1323276/1‑150 (MQ=255)
                                                                                                                                                 ggATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCAGCCATCGc      >  1:1126107/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCTCCAGCGGAATAAAAATCCAGGCCAGGATCATCCAGAAACAGCCGAGCGTCGCGCTGAAAGCCGACGCACCGTGACGACGATAATCGCGATAACGCCCGATAAGCCGCGCGTTGACCGGCGGGATAAGCAACCACCGGGTCAGGATACTCATGATGCACTCCCGACTGGCGTTTTCAGCCCGGAGTAGTTCAACAGGCACCAGTTCGCCAGCGTCAGTATCTCTTCAGCCGCCAGCGCATCACTGCGATATTCACCTACTGGTTGCTTAGCCGCCAGGCATTCAGCCATCGCTTCA  >  CP009273/3688442‑3688740

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: