Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 350,559 G→T 29.0% intergenic (+256/‑74) prpE → / → codB propionate‑‑CoA ligase/cytosine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273350,5590GT29.0% 569.4 / 202.9 289intergenic (+256/‑74)prpE/codBpropionate‑‑CoA ligase/cytosine transporter
Reads supporting (aligned to +/- strand):  ref base G (99/106);  new base T (39/45);  total (138/151)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 295 total aligned reads displayed.

AGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCG  >  CP009273/350410‑350708
                                                                                                                                                     |                                                                                                                                                     
   cgcTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTAtt                                                                                                                                                    <  1:5193819/150‑1 (MQ=255)
    gcTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTg                                                                                                                                                   <  1:899880/150‑1 (MQ=255)
     cTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGc                                                                                                                                                  >  2:763801/1‑150 (MQ=255)
         cgTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCtttt                                                                                                                                              <  2:1072417/150‑1 (MQ=255)
          gTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACTAAAACTATTGCttttt                                                                                                                                             >  1:3287583/1‑150 (MQ=255)
           tCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTa                                                                                                                                            <  1:4145598/150‑1 (MQ=255)
              cATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCt                                                                                                                                         <  1:147782/150‑1 (MQ=255)
              cATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCt                                                                                                                                         >  1:2270187/1‑150 (MQ=255)
                tCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTc                                                                                                                                       <  2:1619719/150‑1 (MQ=255)
                  cGACTCGATGCTCAGTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCATACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAg                                                                                                                                     <  2:3715522/150‑1 (MQ=255)
                  cGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAg                                                                                                                                     <  1:4846540/150‑1 (MQ=255)
                   gACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGa                                                                                                                                    >  2:4040405/1‑150 (MQ=255)
                   gACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGa                                                                                                                                    >  1:4670427/1‑150 (MQ=255)
                     cacCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATg                                                                                                                                  >  2:4517743/1‑150 (MQ=255)
                            cTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGaa                                                                                                                           <  2:3415963/150‑1 (MQ=255)
                             tCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAAt                                                                                                                          <  2:3794717/150‑1 (MQ=255)
                              cAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATg                                                                                                                         <  1:2083747/150‑1 (MQ=255)
                               aGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgc                                                                                                                        <  1:1782057/150‑1 (MQ=255)
                                 aTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcgg                                                                                                                      >  2:3869941/1‑150 (MQ=255)
                                      tGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAt                                                                                                                 >  1:1360416/1‑150 (MQ=255)
                                      tGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAt                                                                                                                 >  1:3997818/1‑150 (MQ=255)
                                        aTGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAttt                                                                                                               <  1:1067657/150‑1 (MQ=255)
                                         tGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAtttt                                                                                                              <  1:3759552/150‑1 (MQ=255)
                                            gACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTg                                                                                                           >  2:433764/1‑150 (MQ=255)
                                            gACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCACACGAAAACTAATGCTTTATATCTTCAGATGAATAGAATGCGGCGGATTTTTTg                                                                                                           >  1:4803901/1‑150 (MQ=255)
                                                  gACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTc                                                                                                     <  2:214755/150‑1 (MQ=255)
                                                       gTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACa                                                                                                >  1:4562776/1‑150 (MQ=255)
                                                        tCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAg                                                                                               <  2:1934019/150‑1 (MQ=255)
                                                          ttATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCa                                                                                             <  1:1729970/150‑1 (MQ=255)
                                                            aTCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCaaa                                                                                           >  1:4351033/1‑150 (MQ=255)
                                                               aGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAg                                                                                        <  2:1390211/150‑1 (MQ=255)
                                                                 gCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAggg                                                                                      <  2:1263258/150‑1 (MQ=255)
                                                                     aCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGAGTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGac                                                                                  >  2:4752379/1‑149 (MQ=255)
                                                                      cccACTGTTTTTACACCGATAATTTTTCCCCCACCTTTGTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAAt                                                                                 <  1:267038/150‑1 (MQ=255)
                                                                       ccACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAAtt                                                                                <  2:170651/150‑1 (MQ=255)
                                                                        cACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAAttt                                                                               <  2:1976362/150‑1 (MQ=255)
                                                                        cACTGATTTTACACCGATAATTTTTCCCCCACCTTATTGCACTCATTCATATAAAAAATATATTTCACCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAAttt                                                                               <  1:1946229/150‑1 (MQ=255)
                                                                          cTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCg                                                                             <  1:3061669/150‑1 (MQ=255)
                                                                            gTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCgtg                                                                           >  2:2982406/1‑150 (MQ=255)
                                                                               tttACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCg                                                                        >  2:3010476/1‑150 (MQ=255)
                                                                                ttACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGc                                                                       <  2:4940816/150‑1 (MQ=255)
                                                                                 tACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCa                                                                      >  2:1576720/1‑150 (MQ=255)
                                                                                 tACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCa                                                                      <  2:1856188/150‑1 (MQ=255)
                                                                                   cacCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAg                                                                    >  1:3449162/1‑150 (MQ=255)
                                                                                     ccGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGAt                                                                  <  2:1686558/150‑1 (MQ=255)
                                                                                     ccGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCAAGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGAt                                                                  <  1:1280551/150‑1 (MQ=255)
                                                                                     ccGATAATTTTTCCCCCACCTTTTTGCACTCAGTCATATAAAAAATATATGTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGAt                                                                  <  2:631825/150‑1 (MQ=255)
                                                                                      cGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATa                                                                 <  1:676843/150‑1 (MQ=255)
                                                                                        aTAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATaac                                                               >  1:2361912/1‑150 (MQ=255)
                                                                                        aTAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATaac                                                               >  2:5231837/1‑150 (MQ=255)
                                                                                          aaTTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATaacaa                                                             >  2:3115226/1‑150 (MQ=255)
                                                                                           aTTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATaacaac                                                            <  1:252128/150‑1 (MQ=255)
                                                                                             atttCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACtt                                                          <  2:3529667/149‑1 (MQ=255)
                                                                                              tttCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACttt                                                         >  2:2801051/1‑150 (MQ=255)
                                                                                                tCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAg                                                       <  2:1033307/150‑1 (MQ=255)
                                                                                                tCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAg                                                       <  2:4632479/150‑1 (MQ=255)
                                                                                                 cccccACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGc                                                      >  2:1158934/1‑150 (MQ=255)
                                                                                                  ccccACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGcc                                                     <  1:1820317/150‑1 (MQ=255)
                                                                                                   cccACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCa                                                    >  2:5264891/1‑150 (MQ=255)
                                                                                                   cccACCTTTTTGCACTCATTCATATAAAAAATATATTTACCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTTGGTTTCAAACAGCAAAAAGGGGGAAATTCGTGTCGCAAGATAACAACTATAGCCa                                                    >  1:4239798/1‑150 (MQ=255)
                                                                                                     cACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAgg                                                  <  2:1459104/150‑1 (MQ=255)
                                                                                                     cACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAgg                                                  <  1:3071447/150‑1 (MQ=255)
                                                                                                        cTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGc                                               >  1:4944063/1‑150 (MQ=255)
                                                                                                           tttGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAg                                            >  2:3255268/1‑150 (MQ=255)
                                                                                                            ttGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGt                                           >  1:599306/1‑150 (MQ=255)
                                                                                                              gCACTCATTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTcc                                         <  2:324447/150‑1 (MQ=255)
                                                                                                               cACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTccc                                        >  1:5317651/1‑150 (MQ=255)
                                                                                                               cACTCATTCATATAAAAAATATATGTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTccc                                        <  2:3829789/150‑1 (MQ=255)
                                                                                                               cACTCAGTCATATAAAAAATATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTccc                                        <  2:3663430/150‑1 (MQ=255)
                                                                                                                       catATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCgg                                <  1:594534/150‑1 (MQ=255)
                                                                                                                            aaaaaaTATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAATGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCgg                           >  2:913098/1‑150 (MQ=255)
                                                                                                                            aaaaaaTATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTGTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCgg                           >  2:2202495/1‑150 (MQ=255)
                                                                                                                             aaaaaTATATTTCCCCACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGa                          >  2:3364437/1‑150 (MQ=255)
                                                                                                                             aaaaaTATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGa                          >  1:2026325/1‑150 (MQ=255)
                                                                                                                              aaaaTATATTTCCCCTCGAATACGAGTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATTACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGAGCGGaa                         <  2:5286321/150‑1 (MQ=255)
                                                                                                                              aaaaTATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGaa                         >  1:4763996/1‑150 (MQ=255)
                                                                                                                               aaaTATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGaaa                        >  2:4810121/1‑150 (MQ=255)
                                                                                                                                aaTATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAg                       <  1:4496452/150‑1 (MQ=255)
                                                                                                                                 atatatTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAgg                      >  2:140850/1‑150 (MQ=255)
                                                                                                                                    tatTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGt                   >  2:4897073/1‑150 (MQ=255)
                                                                                                                                    tatTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGt                   <  1:1558083/150‑1 (MQ=255)
                                                                                                                                     atTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTa                  >  2:2340264/1‑150 (MQ=255)
                                                                                                                                      tttCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTAt                 >  2:2144912/1‑150 (MQ=255)
                                                                                                                                         ccccACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTgg              >  2:549442/1‑150 (MQ=255)
                                                                                                                                          cccACGAAAACTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAACAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGc             <  1:4049150/150‑1 (MQ=255)
                                                                                                                                          cccACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGc             <  1:2471175/150‑1 (MQ=255)
                                                                                                                                           ccACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCa            >  1:3589670/1‑150 (MQ=255)
                                                                                                                                           ccACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCa            <  2:2391394/150‑1 (MQ=255)
                                                                                                                                           ccACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCa            >  1:1009255/1‑150 (MQ=255)
                                                                                                                                            cACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCAt           <  1:407345/150‑1 (MQ=255)
                                                                                                                                               gAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGa        >  2:4183370/1‑150 (MQ=255)
                                                                                                                                               gAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGa        <  2:2048430/150‑1 (MQ=255)
                                                                                                                                                aaaaCTATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGAc       <  1:2576471/150‑1 (MQ=255)
                                                                                                                                                aaaaCGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGAc       <  1:2791430/150‑1 (MQ=255)
                                                                                                                                                aaaaCGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGAc       >  1:382655/1‑150 (MQ=255)
                                                                                                                                                aaaaCGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGAc       >  1:3409981/1‑150 (MQ=255)
                                                                                                                                                 aaaCGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACg      >  1:3179324/1‑150 (MQ=255)
                                                                                                                                                   aCGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGtt    <  1:757648/150‑1 (MQ=255)
                                                                                                                                                     gATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCg  <  2:4133003/150‑1 (MQ=255)
                                                                                                                                                     gATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCg  >  1:1427364/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAAAAAGGGGGAATTTCGTGTCGCAAGATAACAACTTTAGCCAGGGGCCAGTCCCGCAGTCGGCGCGGAAAGGGGTATTGGCATTGACGTTCG  >  CP009273/350410‑350708

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