Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,879,461 A→T 31.9% K342N (AAA→AAT yidC → membrane protein insertase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,879,4610AT31.9% 421.3 / 307.1 295K342N (AAA→AATyidCmembrane protein insertase
Reads supporting (aligned to +/- strand):  ref base A (89/112);  new base T (34/60);  total (123/172)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.06e-01

Only 100 of 300 total aligned reads displayed.

GTACTGGTTCAGCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAGATTCAGGC  >  CP009273/3879312‑3879610
                                                                                                                                                     |                                                                                                                                                     
gTACTGGTTCAGCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCaaa                                                                                                                                                       >  1:1495655/1‑150 (MQ=255)
  aCTGGTTCAGCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAAct                                                                                                                                                     <  2:2496246/150‑1 (MQ=255)
    tGGTACAGCCTGGTCAGACTGGCGCGATGAACAGCAACCTGTGGGTAGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCAGACCGTAGAAAACCGAAGGTTGTGGATCATATCTCAGCCGATGTTAAAActgc                                                                                                                                                   >  1:1458981/1‑150 (MQ=255)
     ggTTCAGCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATctgct                                                                                                                                                  >  1:4610437/1‑150 (MQ=255)
      gTTCAGCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAATATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAActgctg                                                                                                                                                 <  1:2899873/150‑1 (MQ=255)
      gTTCAGCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATctgctg                                                                                                                                                 <  2:3391223/150‑1 (MQ=255)
        tCAGCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGaa                                                                                                                                               <  1:746165/150‑1 (MQ=255)
         cAGCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGaaa                                                                                                                                              >  2:1402131/1‑150 (MQ=255)
           gCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATg                                                                                                                                            <  1:2675413/150‑1 (MQ=255)
            ccTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATgg                                                                                                                                           <  2:3202922/150‑1 (MQ=255)
               ggTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTTGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATc                                                                                                                                        >  2:4115803/1‑150 (MQ=255)
               ggTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATc                                                                                                                                        <  2:4487274/150‑1 (MQ=255)
               ggTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATc                                                                                                                                        <  1:3438971/150‑1 (MQ=255)
                  cAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCAt                                                                                                                                     >  1:189462/1‑150 (MQ=255)
                   aGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATa                                                                                                                                    <  2:2141747/150‑1 (MQ=255)
                   aGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATa                                                                                                                                    <  2:2288080/150‑1 (MQ=255)
                    gACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAg                                                                                                                                   <  1:484899/150‑1 (MQ=255)
                      cTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCAGAGCt                                                                                                                                 <  2:4132338/150‑1 (MQ=255)
                            cgATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTgg                                                                                                                           <  2:3285474/150‑1 (MQ=255)
                              aTGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGt                                                                                                                         <  2:3475637/150‑1 (MQ=255)
                              aTGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGt                                                                                                                         <  1:2312834/150‑1 (MQ=255)
                                gAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTaa                                                                                                                       >  1:3898574/1‑150 (MQ=255)
                                gAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTaa                                                                                                                       <  1:545642/150‑1 (MQ=255)
                                 aaCAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAAc                                                                                                                      >  2:1285105/1‑150 (MQ=255)
                                    aGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgg                                                                                                                   <  1:4159676/150‑1 (MQ=255)
                                    aGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgg                                                                                                                   <  1:1723209/150‑1 (MQ=255)
                                      cACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgggg                                                                                                                 >  2:313152/1‑150 (MQ=255)
                                      cACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgggg                                                                                                                 <  1:1411172/150‑1 (MQ=255)
                                         ccTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCtt                                                                                                              <  2:793035/150‑1 (MQ=255)
                                          cTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTtc                                                                                                             <  2:2714596/150‑1 (MQ=255)
                                            gtgGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTtctc                                                                                                           <  2:26352/150‑1 (MQ=255)
                                            gtgGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTtctc                                                                                                           >  1:1680024/1‑150 (MQ=255)
                                             tgGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGTTTGTGGTTCATCTCTCAGCCGCTTTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTcc                                                                                                          <  1:2540596/150‑1 (MQ=255)
                                                gTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCAtt                                                                                                       >  2:3115360/1‑150 (MQ=255)
                                                  tGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTat                                                                                                     >  1:4308223/1‑150 (MQ=255)
                                                   ggCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTatc                                                                                                    >  1:3727663/1‑150 (MQ=255)
                                                   ggCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTatc                                                                                                    >  1:1892299/1‑150 (MQ=255)
                                                    gCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTatca                                                                                                   >  1:105728/1‑150 (MQ=255)
                                                     cccGGAAATCCAGGACAATCTGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGTAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTatcat                                                                                                  <  2:4576310/150‑1 (MQ=255)
                                                     cccGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTatcat                                                                                                  >  1:1955501/1‑150 (MQ=255)
                                                       cGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTatcatca                                                                                                <  1:1137620/150‑1 (MQ=255)
                                                        ggAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTATCAGCCGCTGTTCAAACTGCTGATATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTatcatcat                                                                                               <  1:3618861/150‑1 (MQ=255)
                                                          aaaTCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTAtcatcatca                                                                                             <  2:789248/150‑1 (MQ=255)
                                                             tCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCt                                                                                          <  1:4503413/150‑1 (MQ=255)
                                                               cAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCttt                                                                                        <  2:556528/150‑1 (MQ=255)
                                                                aGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTa                                                                                       >  1:3924755/1‑150 (MQ=255)
                                                                aGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATAATCATCACCTTTa                                                                                       <  2:1361535/150‑1 (MQ=255)
                                                                  gACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTAtc                                                                                     <  1:828166/150‑1 (MQ=255)
                                                                   aCAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTAtcg                                                                                    <  1:2806931/150‑1 (MQ=255)
                                                                      aaaTGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTAtcgttc                                                                                 <  2:3231438/150‑1 (MQ=255)
                                                                      aaaTGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTAtcgttc                                                                                 <  1:3945155/150‑1 (MQ=255)
                                                                                  ttGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGt                                                                     <  2:1817976/150‑1 (MQ=255)
                                                                                      tCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCAGAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTAccc                                                                 <  1:3158606/150‑1 (MQ=255)
                                                                                      tCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTAccc                                                                 >  1:1333733/1‑150 (MQ=255)
                                                                                         gCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCCTAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCt                                                              >  2:2789738/1‑150 (MQ=255)
                                                                                         gCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCt                                                              <  2:532577/150‑1 (MQ=255)
                                                                                            ccTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGAc                                                           <  2:3771367/150‑1 (MQ=255)
                                                                                             cTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTATCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGAcc                                                          <  1:2252549/150‑1 (MQ=255)
                                                                                                gATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCaaa                                                       <  2:1333990/150‑1 (MQ=255)
                                                                                                gATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCaaa                                                       <  2:1117441/150‑1 (MQ=255)
                                                                                                 aTCTGACCGTTGATTATGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAg                                                      >  2:4545264/1‑150 (MQ=255)
                                                                                                   cTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAgcg                                                    >  1:2500523/1‑150 (MQ=255)
                                                                                                   cTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAgcg                                                    <  2:1576282/150‑1 (MQ=255)
                                                                                                    tGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAgcgc                                                   >  2:876448/1‑150 (MQ=255)
                                                                                                      aCCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAg                                                 <  1:562390/150‑1 (MQ=255)
                                                                                                      aCCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAg                                                 >  2:14592/1‑150 (MQ=255)
                                                                                                        cGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTa                                               >  1:607524/1‑150 (MQ=255)
                                                                                                          ttGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTaca                                             <  2:1967005/150‑1 (MQ=255)
                                                                                                           tGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTacac                                            >  1:952891/1‑150 (MQ=255)
                                                                                                           tGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTacac                                            <  2:3710955/150‑1 (MQ=255)
                                                                                                           tGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTacac                                            >  2:2758689/1‑150 (MQ=255)
                                                                                                             aTTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCt                                          >  1:1226566/1‑150 (MQ=255)
                                                                                                             aTTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCt                                          >  2:1190192/1‑150 (MQ=255)
                                                                                                               tACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTcc                                        <  1:3270241/150‑1 (MQ=255)
                                                                                                                  ggttggttGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATg                                     >  1:1857574/1‑150 (MQ=255)
                                                                                                                     tggttgTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCg                                  <  1:3019108/150‑1 (MQ=255)
                                                                                                                        ttgTGGTTCATCTCTCAGCCGCTGTTCAATCTGCTTAAATGGATCCATAGCTTTGTGGGTATCTGGGGCTTCTCCATTATCATCATCACCTTTTTATTTCGTGGCATAATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAg                               <  1:1351578/150‑1 (MQ=255)
                                                                                                                          gtgGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGAt                             <  1:4362836/150‑1 (MQ=255)
                                                                                                                            ggTTCATCTCTCAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGc                           >  1:4516509/1‑150 (MQ=255)
                                                                                                                            ggTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGc                           >  2:3232839/1‑150 (MQ=255)
                                                                                                                             gTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCg                          <  1:2788019/150‑1 (MQ=255)
                                                                                                                              ttCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGt                         >  2:885784/1‑150 (MQ=255)
                                                                                                                                   ctctcAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGtt                    <  2:3269229/150‑1 (MQ=255)
                                                                                                                                    tctcAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTg                   >  1:206874/1‑150 (MQ=255)
                                                                                                                                     ctcAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGc                  >  2:2982868/1‑150 (MQ=255)
                                                                                                                                     ctcAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGc                  >  1:3644240/1‑150 (MQ=255)
                                                                                                                                      tcAGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCTTCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCa                 >  2:4325874/1‑150 (MQ=255)
                                                                                                                                      tcAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCa                 >  2:918061/1‑150 (MQ=255)
                                                                                                                                        aGCCGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGc               <  2:4314365/150‑1 (MQ=255)
                                                                                                                                          ccGCTGTTCAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCg             <  2:4551649/150‑1 (MQ=255)
                                                                                                                                           cGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGa            >  1:2741919/1‑150 (MQ=255)
                                                                                                                                            gCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGaa           >  1:1825307/1‑150 (MQ=255)
                                                                                                                                             cTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAg          >  1:1757702/1‑150 (MQ=255)
                                                                                                                                             cTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAg          >  1:4214186/1‑150 (MQ=255)
                                                                                                                                               gTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAGAt        <  1:2355534/150‑1 (MQ=255)
                                                                                                                                                 tCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAGATTc      >  1:1629564/1‑150 (MQ=255)
                                                                                                                                                 tCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAGATTc      <  2:1420886/150‑1 (MQ=255)
                                                                                                                                                  cAATCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAGATTCa     <  2:2787186/150‑1 (MQ=255)
                                                                                                                                                  cAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAGATTCa     >  2:759678/1‑150 (MQ=255)
                                                                                                                                                   aaaCTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCATCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACACGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAGATTCAg    >  2:721069/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTACTGGTTCAGCCTGGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCCATTATCATCATCACCTTTATCGTTCGTGGCATCATGTACCCGCTGACCAAAGCGCAGTACACCTCCATGGCGAAGATGCGTATGTTGCAGCCGAAGATTCAGGC  >  CP009273/3879312‑3879610

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: