Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,521,598 C→A 18.2% Q591K (CAG→AAG)  yrfF → inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,521,5980CA18.2% 922.6 / 130.9 302Q591K (CAG→AAG) yrfFinner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (120/127);  new base A (30/25);  total (150/152)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 307 total aligned reads displayed.

CGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCAT  >  CP009273/3521449‑3521747
                                                                                                                                                     |                                                                                                                                                     
cGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATc                                                                                                                                                       <  1:1939957/150‑1 (MQ=255)
 ggTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCa                                                                                                                                                      >  1:2168980/1‑150 (MQ=255)
   tGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGc                                                                                                                                                    <  2:2405937/150‑1 (MQ=255)
     aGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGccc                                                                                                                                                  >  1:25300/1‑150 (MQ=255)
       tGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTg                                                                                                                                                <  2:1105253/150‑1 (MQ=255)
          ggAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGCTTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGcc                                                                                                                                             >  1:1162415/1‑150 (MQ=255)
           gAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCt                                                                                                                                            <  2:3862964/150‑1 (MQ=255)
           gAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCt                                                                                                                                            >  2:3320965/1‑150 (MQ=255)
            aaTCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCtt                                                                                                                                           >  1:731992/1‑150 (MQ=255)
             aTCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTc                                                                                                                                          <  1:388459/150‑1 (MQ=255)
               cGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCgg                                                                                                                                        >  2:3499367/1‑150 (MQ=255)
               cGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCgg                                                                                                                                        <  1:472295/150‑1 (MQ=255)
                gCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCGGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGc                                                                                                                                       >  1:3090633/1‑150 (MQ=255)
                 cTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCa                                                                                                                                      <  1:3516481/150‑1 (MQ=255)
                   ggATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATc                                                                                                                                    >  2:2915957/1‑150 (MQ=255)
                   ggATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATc                                                                                                                                    >  2:597763/1‑150 (MQ=255)
                       aaCCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACtt                                                                                                                                <  2:1782111/150‑1 (MQ=255)
                         ccTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTa                                                                                                                              <  1:2330244/150‑1 (MQ=255)
                         ccTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTa                                                                                                                              <  1:3718573/150‑1 (MQ=255)
                          cTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTAc                                                                                                                             <  2:4183388/150‑1 (MQ=255)
                               ggCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTa                                                                                                                        <  1:3995374/150‑1 (MQ=255)
                                  aaCCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTAccc                                                                                                                     >  2:1354830/1‑150 (MQ=255)
                                     cTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACccgcc                                                                                                                  <  1:925028/150‑1 (MQ=255)
                                     cTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACccgcc                                                                                                                  >  2:912035/1‑150 (MQ=255)
                                     cTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACccgcc                                                                                                                  <  1:1356634/150‑1 (MQ=255)
                                      tCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACccgccg                                                                                                                 <  2:2723484/150‑1 (MQ=255)
                                             cgcCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAAc                                                                                                          >  1:1736267/1‑150 (MQ=255)
                                                cGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGt                                                                                                       <  2:3402523/150‑1 (MQ=255)
                                                  tttATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTgg                                                                                                     <  1:4394055/150‑1 (MQ=255)
                                                  tttATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTgg                                                                                                     <  1:4202418/150‑1 (MQ=255)
                                                    tATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGaa                                                                                                   <  2:3944499/150‑1 (MQ=255)
                                                     aTCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAAc                                                                                                  <  2:3118037/150‑1 (MQ=255)
                                                       cACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGc                                                                                                >  2:4619679/1‑150 (MQ=255)
                                                        aCGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCt                                                                                               <  2:789784/150‑1 (MQ=255)
                                                        aCGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCt                                                                                               <  2:207431/150‑1 (MQ=255)
                                                         cGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCtt                                                                                              <  1:443769/150‑1 (MQ=255)
                                                         cGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCtt                                                                                              <  2:3582765/150‑1 (MQ=255)
                                                           cATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTc                                                                                            <  1:330977/150‑1 (MQ=255)
                                                               aaaCGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGa                                                                                        >  1:977099/1‑150 (MQ=255)
                                                                aaCGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGaa                                                                                       >  2:2494453/1‑150 (MQ=255)
                                                                   ggCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACt                                                                                    <  1:2917741/150‑1 (MQ=255)
                                                                    gcgcGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTg                                                                                   <  1:2362637/150‑1 (MQ=255)
                                                                    gcgcGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTg                                                                                   <  1:4164828/150‑1 (MQ=255)
                                                                         aGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGcaca                                                                              >  2:2613386/1‑150 (MQ=255)
                                                                           cGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACaaa                                                                            <  1:1748430/150‑1 (MQ=255)
                                                                               acaATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCt                                                                        <  1:4294009/150‑1 (MQ=255)
                                                                                caATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTg                                                                       >  1:4504370/1‑150 (MQ=255)
                                                                                caATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTg                                                                       <  2:2344013/150‑1 (MQ=255)
                                                                                 aaTCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGa                                                                      >  1:675454/1‑150 (MQ=255)
                                                                                 aaTCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGa                                                                      >  1:3561563/1‑150 (MQ=255)
                                                                                  aTCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGAt                                                                     <  1:223433/150‑1 (MQ=255)
                                                                                   tCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATGTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATg                                                                    >  2:840389/1‑150 (MQ=255)
                                                                                    cACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGc                                                                   <  1:588366/150‑1 (MQ=255)
                                                                                    cACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGc                                                                   <  1:4643231/150‑1 (MQ=255)
                                                                                       tAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATa                                                                <  1:1878106/150‑1 (MQ=255)
                                                                                        aaaCAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATAc                                                               <  1:3827239/150‑1 (MQ=255)
                                                                                             gTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGt                                                          >  2:1591647/1‑150 (MQ=255)
                                                                                               ccGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGttt                                                        <  2:1534191/150‑1 (MQ=255)
                                                                                                cGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTa                                                       >  2:1421918/1‑150 (MQ=255)
                                                                                                cGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTa                                                       <  2:4540572/150‑1 (MQ=255)
                                                                                                 ggCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTaa                                                      <  1:2888639/150‑1 (MQ=255)
                                                                                                 ggCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTaa                                                      >  2:515793/1‑150 (MQ=255)
                                                                                                 ggCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTaa                                                      >  1:781737/1‑150 (MQ=255)
                                                                                                 ggCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTaa                                                      >  1:3420624/1‑150 (MQ=255)
                                                                                                  gCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAAc                                                     >  2:4287156/1‑150 (MQ=255)
                                                                                                   cccccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACg                                                    >  1:4056921/1‑150 (MQ=255)
                                                                                                   cccccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACg                                                    <  2:2168980/150‑1 (MQ=255)
                                                                                                   cccccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCNTTCCAGAAANTGGCACAANTGCTGATGCATACGCCGNTTAACg                                                    >  1:2652198/1‑150 (MQ=255)
                                                                                                     cccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGcc                                                  >  2:1239652/1‑150 (MQ=255)
                                                                                                     cccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGcc                                                  >  1:4115025/1‑150 (MQ=255)
                                                                                                      ccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCg                                                 <  2:1329291/150‑1 (MQ=255)
                                                                                                       cggcggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGa                                                <  2:2544756/150‑1 (MQ=255)
                                                                                                       cggcggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGa                                                <  2:4364368/150‑1 (MQ=255)
                                                                                                        ggcggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAg                                               <  1:740027/150‑1 (MQ=255)
                                                                                                        ggcggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAg                                               <  1:528243/150‑1 (MQ=255)
                                                                                                          cggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAggg                                             >  2:653379/1‑150 (MQ=255)
                                                                                                             aTTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTAt                                          <  2:2255135/150‑1 (MQ=255)
                                                                                                             aTTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTAt                                          <  1:389758/150‑1 (MQ=255)
                                                                                                              ttCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATc                                         <  2:1023085/150‑1 (MQ=255)
                                                                                                                 cTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTc                                      <  1:4614957/150‑1 (MQ=255)
                                                                                                                 cTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTc                                      >  1:1029381/1‑150 (MQ=255)
                                                                                                                     ttGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACaa                                  <  1:2749029/150‑1 (MQ=255)
                                                                                                                       gTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACaaaa                                <  1:3259336/150‑1 (MQ=255)
                                                                                                                        tCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACaaaaa                               >  2:2025001/1‑150 (MQ=255)
                                                                                                                         cAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAAt                              >  2:1059182/1‑150 (MQ=255)
                                                                                                                         cAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAAt                              <  2:1508900/150‑1 (MQ=255)
                                                                                                                          aGCGACGAAGGAAGCGATTTTGTTGATATGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACATATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATc                             <  1:4072584/150‑1 (MQ=255)
                                                                                                                            cgacgaAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCtt                           >  1:230788/1‑150 (MQ=255)
                                                                                                                             gacgaAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTc                          <  2:724701/150‑1 (MQ=255)
                                                                                                                             gacgaAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTc                          >  1:1512488/1‑150 (MQ=255)
                                                                                                                                   ggCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGAc                    >  2:3604046/1‑150 (MQ=255)
                                                                                                                                       gCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCa                <  1:2742469/150‑1 (MQ=255)
                                                                                                                                          aTTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATg             >  2:3163884/1‑150 (MQ=255)
                                                                                                                                          aTTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATg             <  2:2846146/150‑1 (MQ=255)
                                                                                                                                           ttttGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATgg            <  1:1323143/150‑1 (MQ=255)
                                                                                                                                            tttGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGt           <  2:4085815/150‑1 (MQ=255)
                                                                                                                                             ttgttgATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTa          >  1:1643636/1‑150 (MQ=255)
                                                                                                                                              tgttgATCAGCCCTGGCCTTAGGCATCACTGTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAATCTGGCTCTAATGCTTATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAATCATCTTCAATGACGCCAATGTTAc         <  2:1907182/150‑1 (MQ=255)
                                                                                                                                                   aTCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACgcagc    >  1:3158989/1‑150 (MQ=255)
                                                                                                                                                    tCAGCCCTGGCCTTCTGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACgcagca   >  1:299985/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCAT  >  CP009273/3521449‑3521747

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