Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,373,985 C→A 100% intergenic (+37/‑117) yhcB → / → degQ DUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,373,9850CA98.8% 692.4 / ‑1.4 227intergenic (+37/‑117)yhcB/degQDUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (86/136);  total (89/138)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 231 total aligned reads displayed.

TCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTG  >  CP009273/3373840‑3374129
                                                                                                                                                 |                                                                                                                                                
tCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtatc                                                                                                                                                >  1:626612/1‑148 (MQ=255)
    cTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccc                                                                                                                                            >  1:1094391/1‑148 (MQ=255)
    cTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccc                                                                                                                                            >  1:921300/1‑148 (MQ=255)
       gCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTctccctcc                                                                                                                                         <  1:2640302/148‑1 (MQ=255)
       gCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccctcc                                                                                                                                         <  1:2782268/148‑1 (MQ=255)
       gCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccctcc                                                                                                                                         <  2:1048622/148‑1 (MQ=255)
        cAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCtccctccc                                                                                                                                        <  2:2583597/148‑1 (MQ=255)
          gagTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCg                                                                                                                                      >  1:3522764/1‑148 (MQ=255)
           agTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGa                                                                                                                                     <  1:1494793/148‑1 (MQ=255)
           agTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGa                                                                                                                                     >  2:1901697/1‑148 (MQ=255)
            gTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGAc                                                                                                                                    <  2:2801805/148‑1 (MQ=255)
                 aaGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCa                                                                                                                               >  1:1862304/1‑148 (MQ=255)
                 aaGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCa                                                                                                                               >  2:3098925/1‑148 (MQ=255)
                 aaGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCa                                                                                                                               <  1:1661263/148‑1 (MQ=255)
                 aaGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCa                                                                                                                               >  2:756931/1‑148 (MQ=255)
                  aGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGAGCCCTCCTCTTATCTCCCTCCCCGACTATCAt                                                                                                                              >  2:1083060/1‑148 (MQ=255)
                    ccAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAttt                                                                                                                            >  2:2361935/1‑148 (MQ=255)
                    ccAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAttt                                                                                                                            >  1:2823672/1‑148 (MQ=255)
                      aGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTaa                                                                                                                          <  1:646215/148‑1 (MQ=255)
                       gCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAAt                                                                                                                         <  1:1893626/148‑1 (MQ=255)
                         aaCGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCt                                                                                                                       <  2:2428489/148‑1 (MQ=255)
                         aaCGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCt                                                                                                                       <  1:2373220/148‑1 (MQ=255)
                          aagATCTGGCACCGGTGCAGAGGCCTCGAGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAATCGCGACTTATTTATTTTTCGGGCGCAGCAATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTg                                                                                                                      <  2:3470960/146‑1 (MQ=255)
                           cGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTgg                                                                                                                     >  2:796548/1‑148 (MQ=255)
                             aTCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGgtg                                                                                                                   >  2:3314137/1‑148 (MQ=255)
                              tCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGgtgt                                                                                                                  <  2:3285114/148‑1 (MQ=255)
                              tCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGgtgt                                                                                                                  >  2:2497506/1‑148 (MQ=255)
                               cAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGtc                                                                                                                 <  2:3246878/148‑1 (MQ=255)
                                     ccGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTg                                                                                                           <  2:1176170/148‑1 (MQ=255)
                                      cGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGt                                                                                                          >  1:776120/1‑148 (MQ=255)
                                       ggTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGtt                                                                                                         <  2:1419803/148‑1 (MQ=255)
                                          gCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGc                                                                                                      <  2:1532548/148‑1 (MQ=255)
                                          gCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGc                                                                                                      <  2:2903798/148‑1 (MQ=255)
                                             gATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGt                                                                                                   <  1:1060745/148‑1 (MQ=255)
                                             gATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGt                                                                                                   <  1:379044/148‑1 (MQ=255)
                                              aTGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTc                                                                                                  >  2:3519643/1‑148 (MQ=255)
                                                   tCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGaaa                                                                                             >  2:2225711/1‑148 (MQ=255)
                                                    cgcgACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGAGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGACGTCTGaaaa                                                                                            >  1:1244144/1‑148 (MQ=255)
                                                        aCTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCa                                                                                        >  1:3224569/1‑148 (MQ=255)
                                                          tagtCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAAt                                                                                      <  1:2993772/145‑1 (MQ=255)
                                                          tATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAAt                                                                                      <  1:2657430/148‑1 (MQ=255)
                                                          tATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAAt                                                                                      >  1:678955/1‑148 (MQ=255)
                                                           aTTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATa                                                                                     >  1:1343995/1‑148 (MQ=255)
                                                              cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACa                                                                                  <  1:456501/148‑1 (MQ=255)
                                                              cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACa                                                                                  <  1:3507665/148‑1 (MQ=255)
                                                              cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACa                                                                                  <  1:1818402/148‑1 (MQ=255)
                                                              cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACa                                                                                  <  2:2689653/148‑1 (MQ=255)
                                                                gAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATc                                                                                >  1:3460745/1‑148 (MQ=255)
                                                                   ggCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaaa                                                                             >  2:1617320/1‑148 (MQ=255)
                                                                   ggCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaaa                                                                             <  1:1376303/148‑1 (MQ=255)
                                                                   ggCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaaa                                                                             <  1:1480964/148‑1 (MQ=255)
                                                                   ggCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaaa                                                                             >  2:3514813/1‑148 (MQ=255)
                                                                         tCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGttt                                                                       <  2:532619/148‑1 (MQ=255)
                                                                            ggCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGa                                                                    <  2:2404407/148‑1 (MQ=255)
                                                                              ccTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                  >  1:1935703/1‑148 (MQ=255)
                                                                              ccTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                  <  1:2948198/148‑1 (MQ=255)
                                                                                 gctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtctc                                                               >  2:1988470/1‑148 (MQ=255)
                                                                                   tgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCtt                                                             <  1:3446534/148‑1 (MQ=255)
                                                                                    gcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCttt                                                            <  2:1343995/148‑1 (MQ=255)
                                                                                    gcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCttt                                                            >  2:446124/1‑148 (MQ=255)
                                                                                      gTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTc                                                          <  2:1725524/148‑1 (MQ=255)
                                                                                      gTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTc                                                          <  1:118602/148‑1 (MQ=255)
                                                                                       tACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCt                                                         >  1:661205/1‑148 (MQ=255)
                                                                                       tACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCttttat                                                         >  1:3302393/1‑146 (MQ=255)
                                                                                       tACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCttttat                                                         >  2:2047231/1‑146 (MQ=255)
                                                                                       tACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCt                                                         >  2:2578262/1‑148 (MQ=255)
                                                                                        aCTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCtt                                                        <  1:2781850/148‑1 (MQ=255)
                                                                                        aCTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCAGTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTAtt                                                        <  2:3026710/148‑1 (MQ=255)
                                                                                         cTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTa                                                       <  2:910554/148‑1 (MQ=255)
                                                                                         cTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTa                                                       >  1:2165251/1‑148 (MQ=255)
                                                                                           ggCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATc                                                     <  2:1752326/148‑1 (MQ=255)
                                                                                            gcgcGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAGTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCa                                                    >  1:3135210/1‑148 (MQ=255)
                                                                                             cgcgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCAt                                                   <  1:2238250/148‑1 (MQ=255)
                                                                                                 aaGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAg                                               <  2:1792259/148‑1 (MQ=255)
                                                                                                  aGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAgg                                              >  2:1063073/1‑148 (MQ=255)
                                                                                                   gcgcGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGt                                             <  2:2983656/148‑1 (MQ=255)
                                                                                                    cgcgACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTa                                            <  2:1862304/148‑1 (MQ=255)
                                                                                                    cgcgACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTa                                            <  2:262060/148‑1 (MQ=255)
                                                                                                       gACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACga                                         >  1:2260788/1‑148 (MQ=255)
                                                                                                        aCTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACgag                                        <  1:1723729/148‑1 (MQ=255)
                                                                                                           aaTTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGc                                     <  2:2992987/148‑1 (MQ=255)
                                                                                                            atttatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCa                                    >  1:576181/1‑148 (MQ=255)
                                                                                                            atttatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCa                                    <  1:1617320/148‑1 (MQ=255)
                                                                                                              ttatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAgg                                  <  1:524232/148‑1 (MQ=255)
                                                                                                                atttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaa                                <  1:3418085/148‑1 (MQ=255)
                                                                                                                atttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaa                                <  2:626612/148‑1 (MQ=255)
                                                                                                                  ttttCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaata                              <  1:832138/148‑1 (MQ=255)
                                                                                                                        ggCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaa                        <  1:489036/148‑1 (MQ=255)
                                                                                                                        ggCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaa                        >  2:3306278/1‑148 (MQ=255)
                                                                                                                          cgcAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaaaa                      >  1:874008/1‑148 (MQ=255)
                                                                                                                             aGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaaca                   <  1:1421380/148‑1 (MQ=255)
                                                                                                                                 aTTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAcc               <  2:1266569/148‑1 (MQ=255)
                                                                                                                                  ttGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTAAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAccc              <  1:211514/148‑1 (MQ=255)
                                                                                                                                       cccTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTg         <  1:1688092/148‑1 (MQ=255)
                                                                                                                                        cctcctCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGt        <  2:1562803/148‑1 (MQ=255)
                                                                                                                                        cctcctCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGt        <  2:2681614/148‑1 (MQ=255)
                                                                                                                                        cctcctCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGt        >  2:2148313/1‑148 (MQ=255)
                                                                                                                                           cctcTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGa     >  1:1544695/1‑148 (MQ=255)
                                                                                                                                           cctcTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGa     >  2:2753529/1‑148 (MQ=255)
                                                                                                                                              cttatCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTg  >  1:910498/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTG  >  CP009273/3373840‑3374129

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: ‑