Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,530,481 G→T 14.8% R104L (CGT→CTT)  greB → transcript cleavage factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,530,4810GT14.8% 1014.8 / 122.3 338R104L (CGT→CTT) greBtranscript cleavage factor
Reads supporting (aligned to +/- strand):  ref base G (137/150);  new base T (32/19);  total (169/169)
Fisher's exact test for biased strand distribution p-value = 6.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

CTGCGTGAAATCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCG  >  CP009273/3530333‑3530630
                                                                                                                                                    |                                                                                                                                                     
cTGCTTGAATTCGACCGTCGCGTGCGATATCTCACTAACTGCCTGTAAAATCTCAAAATCGTCGATTACTCCCCTCTGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGt                                                                                                                                                      <  2:2243941/150‑1 (MQ=255)
  gCGTGAAATCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTAt                                                                                                                                                    <  2:1497574/150‑1 (MQ=255)
   cGTGAATTCGACCGTCGCTTGCGCTATCTCACTAAATGCTTGTAAATTCTCAAAATGGTCGATTACACCCCTCAGCTGTAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTAtt                                                                                                                                                   <  1:715364/150‑1 (MQ=255)
   cGTGAAATCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTTCGGAtt                                                                                                                                                   >  1:5006735/1‑150 (MQ=255)
    gTGAAATCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTg                                                                                                                                                  >  2:2622139/1‑150 (MQ=255)
     tGAAATCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGt                                                                                                                                                 >  1:627431/1‑150 (MQ=255)
      taaaTCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTATTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTc                                                                                                                                                <  1:4433560/149‑1 (MQ=255)
      gAACTCCACCGTCTCTTGCGCTATCTCACTAAATGCCTGCAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTc                                                                                                                                                <  1:2826720/150‑1 (MQ=255)
      gAAATCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCCGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTc                                                                                                                                                <  2:4656058/150‑1 (MQ=255)
      gAAATCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTc                                                                                                                                                >  1:3223058/1‑150 (MQ=255)
      gAAATCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTc                                                                                                                                                >  1:4359367/1‑150 (MQ=255)
      gAAAGCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGACGGCAACGTCTTTTTTGGCGCGTGGGTGGTGTTTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTc                                                                                                                                                <  2:3954811/150‑1 (MQ=255)
       aaaTCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGACGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCg                                                                                                                                               <  1:1244330/150‑1 (MQ=255)
        aaTCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCgg                                                                                                                                              >  2:297507/1‑150 (MQ=255)
            gACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCTATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTAc                                                                                                                                          <  1:2524109/150‑1 (MQ=255)
               cGTCGCGTGCGCTATCTCACTAAATGCCTGTAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGAt                                                                                                                                       <  2:637355/150‑1 (MQ=255)
               cGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGTAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGAt                                                                                                                                       <  2:2037615/150‑1 (MQ=255)
                gTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACGCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATg                                                                                                                                      <  2:137061/150‑1 (MQ=255)
                    cgTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAAt                                                                                                                                  >  1:3308811/1‑150 (MQ=255)
                     gTGCGCTATCTAACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAAtt                                                                                                                                 <  1:3955972/150‑1 (MQ=255)
                      tGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAAttt                                                                                                                                >  2:612795/1‑150 (MQ=255)
                        cgcTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCCGCAGGAAGGCAAAGTCTTTTTTGGCGCGGGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAAttttt                                                                                                                              <  1:1517955/150‑1 (MQ=255)
                         gcTATCTCACAAAATGCCTGGAAAATCTCAAAATCGTCGATTACTACACTCAGCAGGAAGGCAAAGTATTATTTGGCGCGAGGGTGGAGATAGAAAACGACGAAGGCGTTACTAACCGAAACCGTATTGTCGGCTACTATGAAAttatag                                                                                                                             >  2:2151480/1‑146 (MQ=255)
                          cTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTgg                                                                                                                            <  1:761852/150‑1 (MQ=255)
                               tcACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTa                                                                                                                       >  1:4095890/1‑150 (MQ=255)
                                cACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTaa                                                                                                                      >  2:3434609/1‑150 (MQ=255)
                                cACTAAATGCCTGGAAAAACTCAAAATCGTCGATTACTCCCCTCACCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGATATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTaa                                                                                                                      <  2:1617406/150‑1 (MQ=255)
                                  cTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAg                                                                                                                    >  1:850152/1‑150 (MQ=255)
                                    aaaTGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGAt                                                                                                                  >  1:4762819/1‑150 (MQ=255)
                                    aaaTGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGAt                                                                                                                  >  1:4103898/1‑150 (MQ=255)
                                     aattcatGGAAATTCACAAAATCGTCTATTACTCCCCTCAGCAGGAAGGCAAAGTATTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGAtt                                                                                                                 <  2:4732134/144‑1 (MQ=255)
                                     aaTGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGAtt                                                                                                                 >  1:2299754/1‑150 (MQ=255)
                                       tGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTAc                                                                                                               <  1:1714985/150‑1 (MQ=255)
                                          cTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACAtc                                                                                                            >  2:4172638/1‑150 (MQ=255)
                                           tGGAAAATCTCAACATCGTCGATTACTCCCCTCAGCAGGAAGGCAACGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGATATTTTTGGCCGTAAAGATTACAtct                                                                                                           <  1:3015317/150‑1 (MQ=255)
                                           tGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACAtct                                                                                                           >  1:4655253/1‑150 (MQ=255)
                                               aaaTCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTAt                                                                                                       >  1:3171534/1‑150 (MQ=255)
                                                aaTCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATc                                                                                                      >  2:410055/1‑150 (MQ=255)
                                                aaTCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATc                                                                                                      >  2:2828837/1‑150 (MQ=255)
                                                  tctcAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCTTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACAAATCTATCGa                                                                                                    >  2:2685759/1‑150 (MQ=255)
                                                   ctcACAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGAt                                                                                                   <  1:2128278/150‑1 (MQ=255)
                                                   ctcAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGAt                                                                                                   <  2:2749631/150‑1 (MQ=255)
                                                     cAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTc                                                                                                 >  2:3323004/1‑150 (MQ=255)
                                                     cAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGAGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACAAATCTATCGATTc                                                                                                 >  2:4360661/1‑150 (MQ=255)
                                                      aaaaTCGTCGATTACTCCCCTCAGCAGTAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTTAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTcc                                                                                                <  1:1009338/150‑1 (MQ=255)
                                                       aaaTCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTccc                                                                                               >  1:326581/1‑150 (MQ=255)
                                                              cgATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGc                                                                                        <  2:1895583/150‑1 (MQ=255)
                                                               gATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGcc                                                                                       >  1:494520/1‑150 (MQ=255)
                                                                aTTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGccc                                                                                      >  2:3688710/1‑150 (MQ=255)
                                                                 ttACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCcg                                                                                     >  2:1213589/1‑150 (MQ=255)
                                                                 ttACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCcg                                                                                     >  2:2114574/1‑150 (MQ=255)
                                                                 ttACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCcg                                                                                     >  1:1961817/1‑150 (MQ=255)
                                                                    cTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCgcgc                                                                                  >  2:4230890/1‑150 (MQ=255)
                                                                     tCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAAATTTTGTCCGTAAAGATTACATCTCTATAGATTCCCCGATGTCCCGCGCa                                                                                 >  2:4795226/1‑150 (MQ=255)
                                                                     tCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCa                                                                                 >  2:2310023/1‑150 (MQ=255)
                                                                      ccccTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCAt                                                                                >  2:2423107/1‑150 (MQ=255)
                                                                       cccTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCAtt                                                                               <  2:1043422/150‑1 (MQ=255)
                                                                         cTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGc                                                                             >  1:3393166/1‑150 (MQ=255)
                                                                           cagcagGAATGAAATGTCTTTTTTGGCGCTTGGGTGGAGATTTAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTg                                                                           <  1:2016086/150‑1 (MQ=255)
                                                                            agcagGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGa                                                                          >  2:1732173/1‑150 (MQ=255)
                                                                            agcagGAAGGCAAAGTATTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGa                                                                          <  1:1043337/150‑1 (MQ=255)
                                                                              ccggcaGGCACAGTCTTTTTGGGCGCGTGGGTGGAGATGGACACCGACGAGGGCGGGACTCCCCGTTTCCGGATGGTCGGCTACGATGAAATTTTTGGCCGTACAGATTACATCTCTAGCGATTCCCCGCTGGCCCGCGCATTGCTGaaa                                                                        <  2:3558042/145‑1 (MQ=255)
                                                                               agGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGaaaa                                                                       >  1:3097012/1‑150 (MQ=255)
                                                                                      caatttCTTTTTTGGCGGTTGGTTGTAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGt                                                                <  1:2285155/145‑1 (MQ=255)
                                                                                      cAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGt                                                                >  1:564421/1‑150 (MQ=255)
                                                                                       aaaGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTc                                                               >  1:142789/1‑150 (MQ=255)
                                                                                       aaaGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTc                                                               <  2:494520/150‑1 (MQ=255)
                                                                                       aaaGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTc                                                               <  1:2486562/150‑1 (MQ=255)
                                                                                           tCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCAACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGATGAAAAAATAAGTCGGCg                                                           >  2:510654/1‑150 (MQ=255)
                                                                                                tttGGAGCTTGGTTGTATTTTTAAAACGTAGATGTCGTGACTCACCGTTTCCGTATTTTCGGCTACGATGAAATTTTTGGACGTAAAGATTACATCTCTATCGATTCCCCTATGGCCCGCGCATTTCTGAAAAAAGAAGTCGGCGATCTg                                                      <  2:776971/150‑1 (MQ=255)
                                                                                                  tGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGc                                                    <  1:886960/150‑1 (MQ=255)
                                                                                                     cgcgTGGGTGTAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCggt                                                 <  2:3770307/150‑1 (MQ=255)
                                                                                                      gcgTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCggtg                                                <  1:2898073/150‑1 (MQ=255)
                                                                                                       cgTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCggtgg                                               >  1:3884810/1‑150 (MQ=255)
                                                                                                       cgTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGACCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCggtgg                                               >  2:3156619/1‑150 (MQ=255)
                                                                                                         tgggtggAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGgtggtg                                             <  2:5177420/150‑1 (MQ=255)
                                                                                                          gggtggAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGa                                            >  1:3650126/1‑150 (MQ=255)
                                                                                                             tggAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATa                                         >  1:4747529/1‑150 (MQ=255)
                                                                                                             tggAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATa                                         <  1:1788860/150‑1 (MQ=255)
                                                                                                              ggAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCGGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATAc                                        >  1:3733775/1‑150 (MQ=255)
                                                                                                                 gaTTGAAAACGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATAcccc                                     >  1:5114383/1‑150 (MQ=255)
                                                                                                                 gaTTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCACCGATGGCCCGCGCATTGATGAAAAAAGAATACGGCGATCTGGAGGTTGTGAATaaacc                                     >  2:4770261/1‑146 (MQ=255)
                                                                                                                  aTTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCt                                    >  2:3251/1‑150 (MQ=255)
                                                                                                                  aTTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCt                                    >  1:2057902/1‑150 (MQ=255)
                                                                                                                      aaaaCGACGATGGCGTGACTCACCGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCg                                >  1:3349156/1‑150 (MQ=255)
                                                                                                                        aaCGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCggg                              <  2:937144/150‑1 (MQ=255)
                                                                                                                        aaCGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCggg                              >  2:5055899/1‑150 (MQ=255)
                                                                                                                                 ggCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGc                     >  1:2617967/1‑150 (MQ=255)
                                                                                                                                  gCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCt                    <  1:3501937/150‑1 (MQ=255)
                                                                                                                                   cGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTg                   >  2:307356/1‑150 (MQ=255)
                                                                                                                                      gACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTa                <  1:183737/150‑1 (MQ=255)
                                                                                                                                      gACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTa                <  1:4172638/150‑1 (MQ=255)
                                                                                                                                        cTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATg              <  1:1143547/150‑1 (MQ=255)
                                                                                                                                          cACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGAGGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGtt            <  2:1011892/150‑1 (MQ=255)
                                                                                                                                           cccGTTTCCGTATGGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTa           <  2:4904099/149‑1 (MQ=255)
                                                                                                                                           aCCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTa           >  2:431469/1‑150 (MQ=255)
                                                                                                                                            ccGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTaa          >  1:4301518/1‑150 (MQ=255)
                                                                                                                                             cGTTTCCTTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAAt         <  1:1542954/150‑1 (MQ=255)
                                                                                                                                              gTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATg        <  1:4830257/150‑1 (MQ=255)
                                                                                                                                                   cGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATc   <  1:4063580/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCGTGAAATCGACCGTCGCGTGCGCTATCTCACTAAATGCCTGGAAAATCTCAAAATCGTCGATTACTCCCCTCAGCAGGAAGGCAAAGTCTTTTTTGGCGCGTGGGTGGAGATTGAAAACGACGATGGCGTGACTCACCGTTTCCGTATTGTCGGCTACGATGAAATTTTTGGCCGTAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCG  >  CP009273/3530333‑3530630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: