Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,068,091 G→T 44.0% Q77K (CAG→AAG)  rutB ← ureidoacrylate amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,068,0910GT44.0% 258.3 / inf 265Q77K (CAG→AAG) rutBureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base G (57/91);  new base T (58/58);  total (116/149)
Fisher's exact test for biased strand distribution p-value = 7.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

AGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGT  >  CP009273/1067942‑1068240
                                                                                                                                                     |                                                                                                                                                     
aGGCTGCGGCACCAGTTCATCCACCCGTTGATAATCCCAGGAGCCTTTCGCCAGCACTATCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGTAAATGCGGTGAGCCGGGTCCGCCAGCCTCGACATACTg                                                                                                                                                       >  1:4331677/1‑150 (MQ=255)
  gCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACtttt                                                                                                                                                     <  1:2522146/150‑1 (MQ=255)
  gCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACtttt                                                                                                                                                     >  1:1612270/1‑150 (MQ=255)
   cTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTc                                                                                                                                                    <  2:4303628/150‑1 (MQ=255)
   cTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTc                                                                                                                                                    <  1:3245369/150‑1 (MQ=255)
   cTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTc                                                                                                                                                    <  2:4030315/150‑1 (MQ=255)
    tGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCa                                                                                                                                                   <  1:663678/150‑1 (MQ=255)
     gCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCAt                                                                                                                                                  <  2:3165860/150‑1 (MQ=255)
      cGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATc                                                                                                                                                 >  1:2295866/1‑150 (MQ=255)
       ggCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATcc                                                                                                                                                <  2:926245/150‑1 (MQ=255)
        gCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATccc                                                                                                                                               <  2:549422/150‑1 (MQ=255)
         cACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCa                                                                                                                                              >  2:2156057/1‑150 (MQ=255)
          acccggtCATCCACCAGTCGATAATCCCAGGAGCCTTTCGCCAGCAATGTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCac                                                                                                                                             <  1:2279765/144‑2 (MQ=255)
          aCCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAg                                                                                                                                             >  1:1007498/1‑150 (MQ=255)
          aCCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAg                                                                                                                                             <  1:833456/150‑1 (MQ=255)
          aCCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAg                                                                                                                                             <  2:532276/150‑1 (MQ=255)
               ttCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCAtt                                                                                                                                        >  2:3129877/1‑150 (MQ=255)
               ttCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATCGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCAtt                                                                                                                                        <  2:1615754/150‑1 (MQ=255)
                     caccaGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGaaa                                                                                                                                  <  1:1415893/150‑1 (MQ=255)
                       ccaGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAAcc                                                                                                                                <  2:640615/150‑1 (MQ=255)
                          gTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAga                                                                                                                             >  2:46948/1‑150 (MQ=255)
                             gATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAgatga                                                                                                                          <  1:3386925/150‑1 (MQ=255)
                              aTAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAgatgat                                                                                                                         <  1:4063427/150‑1 (MQ=255)
                                aaTCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCa                                                                                                                       <  2:15133/150‑1 (MQ=255)
                                aaTCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCa                                                                                                                       <  1:258191/150‑1 (MQ=255)
                                aaTCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTCACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCa                                                                                                                       <  2:791078/150‑1 (MQ=255)
                                 aTCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAg                                                                                                                      <  1:3114806/150‑1 (MQ=255)
                                   cccAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCa                                                                                                                    >  1:4263355/1‑150 (MQ=255)
                                     cAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATc                                                                                                                  >  1:4588935/1‑150 (MQ=255)
                                      aGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATcc                                                                                                                 >  1:3234933/1‑150 (MQ=255)
                                      aGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATcc                                                                                                                 >  1:3990305/1‑150 (MQ=255)
                                       ggAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCTTTTTGAAACCAGATGATCAGCATccc                                                                                                                <  2:2509634/150‑1 (MQ=255)
                                       ggAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATccc                                                                                                                <  1:2664369/150‑1 (MQ=255)
                                        gAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCGGCATCCCt                                                                                                               <  1:1921974/150‑1 (MQ=255)
                                         aGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTg                                                                                                              <  1:2236604/150‑1 (MQ=255)
                                         aGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTg                                                                                                              <  1:4396322/150‑1 (MQ=255)
                                              ttCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCt                                                                                                         >  2:3955395/1‑150 (MQ=255)
                                              ttCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCt                                                                                                         <  1:3464462/150‑1 (MQ=255)
                                                       aaTTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCg                                                                                                >  1:1522674/1‑150 (MQ=255)
                                                        aTTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCgg                                                                                               >  2:3357064/1‑150 (MQ=255)
                                                         tttCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGt                                                                                              <  1:608089/150‑1 (MQ=255)
                                                         tttCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGt                                                                                              <  1:2446430/150‑1 (MQ=255)
                                                           tCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCa                                                                                            <  2:1037153/150‑1 (MQ=255)
                                                             cccTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCAcg                                                                                          <  1:2343558/150‑1 (MQ=255)
                                                              ccTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCAcgg                                                                                         >  1:1203446/1‑150 (MQ=255)
                                                               cTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCAcggc                                                                                        <  2:2540415/150‑1 (MQ=255)
                                                                  cAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGt                                                                                     <  1:598095/150‑1 (MQ=255)
                                                                    gCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGttt                                                                                   <  2:117307/150‑1 (MQ=255)
                                                                    gCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGttt                                                                                   <  1:1082906/150‑1 (MQ=255)
                                                                      tGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGa                                                                                 >  1:1871710/1‑150 (MQ=255)
                                                                           ctgctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGtt                                                                            <  2:269454/150‑1 (MQ=255)
                                                                           ctgctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGtt                                                                            <  1:2408923/150‑1 (MQ=255)
                                                                            tgctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTg                                                                           >  1:216665/1‑150 (MQ=255)
                                                                             gctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTgg                                                                          >  2:2022538/1‑150 (MQ=255)
                                                                              ctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGc                                                                         >  1:1693672/1‑150 (MQ=255)
                                                                               ttACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCa                                                                        >  1:1056817/1‑150 (MQ=255)
                                                                               ttACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGACTGTTGGCa                                                                        >  2:690391/1‑150 (MQ=255)
                                                                               ttACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCa                                                                        >  2:4141038/1‑150 (MQ=255)
                                                                                tACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCaa                                                                       >  2:2226206/1‑150 (MQ=255)
                                                                                 aCGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAAt                                                                      <  1:1862786/150‑1 (MQ=255)
                                                                                  cGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATg                                                                     >  2:2222999/1‑150 (MQ=255)
                                                                                  cGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATg                                                                     >  2:1846896/1‑150 (MQ=255)
                                                                                   gCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGa                                                                    >  2:2919388/1‑150 (MQ=255)
                                                                                      tGGTTTTCAGGGCGTTTGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCg                                                                 <  2:2188598/150‑1 (MQ=255)
                                                                                      tGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCg                                                                 <  1:2206411/150‑1 (MQ=255)
                                                                                      tGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCg                                                                 <  1:2485459/150‑1 (MQ=255)
                                                                                       ggTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCgg                                                                >  2:2898751/1‑150 (MQ=255)
                                                                                       ggTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCgg                                                                <  1:3633784/150‑1 (MQ=255)
                                                                                        gTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCggg                                                               <  2:288252/150‑1 (MQ=255)
                                                                                            tCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAg                                                           <  2:717126/150‑1 (MQ=255)
                                                                                            tCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAg                                                           >  2:1682579/1‑150 (MQ=255)
                                                                                               gggCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGTAATGACCGGGCGAGTgg                                                        >  1:3785076/1‑150 (MQ=255)
                                                                                               gggCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTgg                                                        <  2:3720038/150‑1 (MQ=255)
                                                                                                 gCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGtt                                                      <  1:3935923/150‑1 (MQ=255)
                                                                                                       gATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACa                                                >  2:3595238/1‑150 (MQ=255)
                                                                                                        aTTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTCTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACAt                                               >  2:1281159/1‑150 (MQ=255)
                                                                                                          ttATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCa                                             <  2:4537993/150‑1 (MQ=255)
                                                                                                             tGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAAc                                          <  2:299656/150‑1 (MQ=255)
                                                                                                             tGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAAc                                          >  2:23297/1‑150 (MQ=255)
                                                                                                             tGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAAc                                          >  1:1148906/1‑150 (MQ=255)
                                                                                                             tGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAAc                                          >  2:3508216/1‑150 (MQ=255)
                                                                                                               aaaaTTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAAccc                                        >  1:4223466/1‑150 (MQ=255)
                                                                                                               aaaaTTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAAccc                                        <  1:227231/150‑1 (MQ=255)
                                                                                                                      ggTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCgaga                                 >  2:3372547/1‑150 (MQ=255)
                                                                                                                       gTGAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGAt                                <  2:1964770/150‑1 (MQ=255)
                                                                                                                         gAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCt                              >  1:2418768/1‑150 (MQ=255)
                                                                                                                         gAGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCt                              <  2:1025143/150‑1 (MQ=255)
                                                                                                                          aGCCGGGTCCGCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTa                             >  2:2773524/1‑150 (MQ=255)
                                                                                                                          aGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCACTGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTa                             <  2:39241/150‑1 (MQ=255)
                                                                                                                            ccGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAg                           <  2:4098485/150‑1 (MQ=255)
                                                                                                                                    gCCAGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAgccgcc                   >  1:3698186/1‑150 (MQ=255)
                                                                                                                                       aGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTgg                >  1:3609828/1‑150 (MQ=255)
                                                                                                                                       aGCCTCGACATACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTgg                >  2:3680509/1‑150 (MQ=255)
                                                                                                                                           tCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTg            >  1:3023025/1‑150 (MQ=255)
                                                                                                                                              aCATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCa         <  1:710241/150‑1 (MQ=255)
                                                                                                                                               cATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCAt        <  2:4058537/150‑1 (MQ=255)
                                                                                                                                                 tACTTTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATaa      <  1:3363836/150‑1 (MQ=255)
                                                                                                                                                  aCTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAg     <  2:1128727/150‑1 (MQ=255)
                                                                                                                                                   cTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGc    <  2:381897/150‑1 (MQ=255)
                                                                                                                                                     gTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGt  >  2:4545591/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGT  >  CP009273/1067942‑1068240

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: