Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 652,865 Δ1 bp 100% coding (118/210 nt) cspE → constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273652,8650G.96.5% 1135.4 / 29.4 226coding (118/210 nt)cspEconstitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base . (104/114);  total (106/120)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 230 total aligned reads displayed.

TGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTT  >  CP009273/652719‑653005
                                                                                                                                                  |                                                                                                                                            
tgtCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATgt                                                                                                                                             <  2:3094246/148‑2 (MQ=255)
tgtCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATgt                                                                                                                                             <  1:3156300/148‑2 (MQ=255)
  tCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGttt                                                                                                                                           <  1:1131050/148‑1 (MQ=255)
  tCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGttt                                                                                                                                           <  1:3252054/148‑1 (MQ=255)
   cTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGtttt                                                                                                                                          <  2:2432608/148‑1 (MQ=255)
    tATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGtttta                                                                                                                                         <  2:2509011/148‑2 (MQ=255)
       ttttCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTaaaa                                                                                                                                     >  1:2474530/1‑148 (MQ=255)
       ttttCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTaaaa                                                                                                                                     >  1:2179683/1‑148 (MQ=255)
       ttttCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTaaaa                                                                                                                                     >  1:2528680/1‑148 (MQ=255)
         ttCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAAct                                                                                                                                   <  1:2118071/148‑1 (MQ=255)
           cATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAActct                                                                                                                                 <  2:2979936/148‑1 (MQ=255)
             tGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTg                                                                                                                               <  2:751437/148‑1 (MQ=255)
             tGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTg                                                                                                                               <  2:751256/148‑1 (MQ=255)
              gTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTg                                                                                                                               >  2:1404075/1‑147 (MQ=255)
              gTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTg                                                                                                                               <  1:1404075/147‑1 (MQ=255)
              gTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGc                                                                                                                              >  1:3499164/1‑148 (MQ=255)
                 aaGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGa                                                                                                                           <  1:2275281/148‑1 (MQ=255)
                    gTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAgg                                                                                                                        >  1:2224214/1‑148 (MQ=255)
                     tAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGt                                                                                                                       <  2:1313059/148‑1 (MQ=255)
                     tAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGt                                                                                                                       >  1:891354/1‑148 (MQ=255)
                       aTTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTATCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCa                                                                                                                     >  2:2335181/1‑148 (MQ=255)
                        ttttGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAg                                                                                                                    <  2:2323186/148‑1 (MQ=255)
                        ttttGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAg                                                                                                                    >  2:1359823/1‑148 (MQ=255)
                             aTGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGt                                                                                                               >  2:566275/1‑148 (MQ=255)
                               gTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTag                                                                                                             <  1:168416/148‑1 (MQ=255)
                                      aTTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGa                                                                                                      >  2:1679662/1‑148 (MQ=255)
                                      aTTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGa                                                                                                      >  1:1356464/1‑148 (MQ=255)
                                      aTTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGa                                                                                                      >  1:1384819/1‑148 (MQ=255)
                                      aTTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGa                                                                                                      >  1:2805657/1‑148 (MQ=255)
                                      aTTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGa                                                                                                      >  2:504433/1‑148 (MQ=255)
                                       ttAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGaa                                                                                                     >  2:475417/1‑148 (MQ=255)
                                        tAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGaaa                                                                                                    >  2:1689918/1‑148 (MQ=255)
                                              tAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTa                                                                                              <  1:553567/148‑1 (MQ=255)
                                                aCGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAAc                                                                                            >  2:819000/1‑148 (MQ=255)
                                                 cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACg                                                                                           >  1:2822281/1‑148 (MQ=255)
                                                    tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTg                                                                                        >  1:167339/1‑148 (MQ=255)
                                                     aaGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGc                                                                                       <  2:3034871/148‑1 (MQ=255)
                                                      aGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGcc                                                                                      >  1:1235406/1‑148 (MQ=255)
                                                             tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGcc                                                                               >  2:2746376/1‑148 (MQ=255)
                                                             tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGcc                                                                               <  1:549640/148‑1 (MQ=255)
                                                              aaTGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGccc                                                                              <  1:1843302/148‑1 (MQ=255)
                                                                     ccAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTctgct                                                                       >  1:1729953/1‑148 (MQ=255)
                                                                         aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCaa                                                                   >  2:1394021/1‑148 (MQ=255)
                                                                         aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCaa                                                                   <  2:2791800/148‑1 (MQ=255)
                                                                         aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCaa                                                                   >  1:354311/1‑148 (MQ=255)
                                                                         aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCaa                                                                   >  1:147008/1‑148 (MQ=255)
                                                                         aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCaa                                                                   >  1:1922845/1‑148 (MQ=255)
                                                                          ggATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCaaa                                                                  <  1:431284/148‑1 (MQ=255)
                                                                            aTTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACg                                                                >  1:463954/1‑148 (MQ=255)
                                                                            aTTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACg                                                                <  1:237019/148‑1 (MQ=255)
                                                                              tCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTa                                                              >  1:1217636/1‑148 (MQ=255)
                                                                               cGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTaa                                                             <  2:2474530/148‑1 (MQ=255)
                                                                               cGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTaa                                                             <  2:3504613/148‑1 (MQ=255)
                                                                                  tttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCg                                                          <  1:1640906/148‑1 (MQ=255)
                                                                                   ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGc                                                         <  2:1535080/148‑1 (MQ=255)
                                                                                        tACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGt                                                    <  1:1060337/148‑1 (MQ=255)
                                                                                        tACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGt                                                    <  1:1590936/148‑1 (MQ=255)
                                                                                        tACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGt                                                    <  1:2901479/148‑1 (MQ=255)
                                                                                         aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTa                                                   <  1:843385/148‑1 (MQ=255)
                                                                                          cTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTaa                                                  <  2:866122/148‑1 (MQ=255)
                                                                                           tCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAg                                                 <  1:2581915/148‑1 (MQ=255)
                                                                                           tCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAg                                                 <  1:2290619/148‑1 (MQ=255)
                                                                                           tCCGGAAGACGACAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAg                                                 <  2:1786801/148‑1 (MQ=255)
                                                                                            ccGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGa                                                >  1:2462192/1‑148 (MQ=255)
                                                                                              ggAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCATACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAGTCACTAACGGTGCCAAAGGCCCTTCTGCTGCACACGTAATCGCTCTGTCAGATa                                              >  1:337102/1‑148 (MQ=255)
                                                                                              ggAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATa                                              <  2:769010/148‑1 (MQ=255)
                                                                                              ggAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATa                                              <  1:2727895/148‑1 (MQ=255)
                                                                                                aaGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACg                                            <  1:3462408/148‑1 (MQ=255)
                                                                                                  gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTc                                          >  2:1646014/1‑148 (MQ=255)
                                                                                                   aCGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCa                                         >  2:408059/1‑148 (MQ=255)
                                                                                                       cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCaa                                     >  2:2188025/1‑148 (MQ=255)
                                                                                                        agcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAg                                    <  2:2633430/148‑1 (MQ=255)
                                                                                                         gcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGa                                   >  1:345748/1‑148 (MQ=255)
                                                                                                           aaaGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAAt                                 <  2:3437642/148‑1 (MQ=255)
                                                                                                             aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTc                               <  2:530497/148‑1 (MQ=255)
                                                                                                               aCGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCaa                             <  2:2951448/148‑1 (MQ=255)
                                                                                                               aCGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCaa                             <  2:3006000/148‑1 (MQ=255)
                                                                                                                cGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCaaa                            <  2:3339397/148‑1 (MQ=255)
                                                                                                                cGTGTTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCaaa                            <  1:92607/148‑1 (MQ=255)
                                                                                                                  tgtTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAAc                          >  2:3046373/1‑148 (MQ=255)
                                                                                                                  tgtTCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAAc                          <  1:1114627/148‑1 (MQ=255)
                                                                                                                    ttCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAAccc                        >  2:3180382/1‑148 (MQ=255)
                                                                                                                    ttCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAAccc                        >  2:920657/1‑148 (MQ=255)
                                                                                                                     tCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCg                       <  1:2675714/148‑1 (MQ=255)
                                                                                                                     tCGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCg                       <  1:850069/148‑1 (MQ=255)
                                                                                                                      cGTACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGc                      >  1:1422776/1‑148 (MQ=255)
                                                                                                                        tACACTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCtt                    >  2:3282767/1‑148 (MQ=255)
                                                                                                                           acTTCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAAt                 <  2:696393/148‑1 (MQ=255)
                                                                                                                             ttCTCTGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCa               >  2:2571297/1‑148 (MQ=255)
                                                                                                                                  tGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg          <  1:939591/148‑1 (MQ=255)
                                                                                                                                  tGCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg          <  2:2187316/148‑1 (MQ=255)
                                                                                                                                   gCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGt         >  2:1599690/1‑148 (MQ=255)
                                                                                                                                   gCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGt         >  2:1599281/1‑148 (MQ=255)
                                                                                                                                   gCAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGt         <  2:1439295/148‑1 (MQ=255)
                                                                                                                                    cAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt        >  1:3297163/1‑148 (MQ=255)
                                                                                                                                    cAATCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt        <  1:2335062/148‑1 (MQ=255)
                                                                                                                                     aaTCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttt       >  2:1735965/1‑148 (MQ=255)
                                                                                                                                      aTCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttt      >  1:2936462/1‑148 (MQ=255)
                                                                                                                                       tCCAGACTAATGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt     <  2:2366836/148‑1 (MQ=255)
                                                                                                                                         cAGACTAATGTTTTAAAACTCGTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttttt   <  2:1932288/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTT  >  CP009273/652719‑653005

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: