Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,833,730 C→A 86.3% R3S (CGT→AGT)  ynjE → molybdopterin synthase sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,833,7300CA86.3% 775.5 / 126.8 299R3S (CGT→AGT) ynjEmolybdopterin synthase sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (25/16);  new base A (128/130);  total (153/146)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

Only 100 of 307 total aligned reads displayed.

CTCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTATAAC  >  CP009273/1833583‑1833879
                                                                                                                                                   |                                                                                                                                                     
ctcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgt                                                                                                                                                     >  1:3030232/1‑150 (MQ=255)
ctcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgt                                                                                                                                                     <  1:1988444/150‑1 (MQ=255)
  ccAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtgt                                                                                                                                                   >  1:4432062/1‑150 (MQ=255)
  ccAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtgt                                                                                                                                                   >  2:2602662/1‑150 (MQ=255)
  ccAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtgt                                                                                                                                                   <  2:1188053/150‑1 (MQ=255)
  ccAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtgt                                                                                                                                                   <  1:447570/150‑1 (MQ=255)
   cAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGtt                                                                                                                                                  <  1:2323756/150‑1 (MQ=255)
    acgaagTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGACAGAAAATTACAACCAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGttt                                                                                                                                                 <  2:1825037/145‑1 (MQ=255)
     ggATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTtc                                                                                                                                                >  2:1950096/1‑150 (MQ=255)
       aTGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTtctc                                                                                                                                              >  1:2161531/1‑150 (MQ=255)
        tGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCa                                                                                                                                             >  2:2901313/1‑150 (MQ=255)
            ccTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATg                                                                                                                                         <  1:497920/150‑1 (MQ=255)
              tgctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGAc                                                                                                                                       <  2:15508/150‑1 (MQ=255)
                ctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcg                                                                                                                                     <  1:342736/150‑1 (MQ=255)
                 tgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcgc                                                                                                                                    >  1:2205879/1‑150 (MQ=255)
                  gATAGTTCTGTTCTGTATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcgcg                                                                                                                                   <  1:362563/150‑1 (MQ=255)
                  gATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcgcg                                                                                                                                   >  1:3686376/1‑150 (MQ=255)
                  gATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACcgcg                                                                                                                                   >  1:1376994/1‑150 (MQ=255)
                   aTAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCgcgc                                                                                                                                  >  1:1520811/1‑150 (MQ=255)
                    tAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCgcgcc                                                                                                                                 <  2:587248/150‑2 (MQ=255)
                      gttctgttctGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTgg                                                                                                                               >  1:2430900/1‑150 (MQ=255)
                      gttctgttctGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTgg                                                                                                                               <  2:3741241/150‑1 (MQ=255)
                         ctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCaa                                                                                                                            <  2:3019887/150‑1 (MQ=255)
                         ctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCaa                                                                                                                            >  2:706699/1‑150 (MQ=255)
                          tgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAAt                                                                                                                           <  2:1008720/150‑1 (MQ=255)
                           gttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATg                                                                                                                          <  1:1186545/150‑1 (MQ=255)
                             tctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGc                                                                                                                        <  2:1861049/150‑1 (MQ=255)
                               tgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCtt                                                                                                                      <  2:1394401/150‑1 (MQ=255)
                                ggATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                     >  2:220175/1‑150 (MQ=255)
                                 gATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTa                                                                                                                    >  2:1043306/1‑150 (MQ=255)
                                 gATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTa                                                                                                                    >  1:608154/1‑150 (MQ=255)
                                     tGGCGCAGTTGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCACTGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggc                                                                                                                <  1:3843297/150‑1 (MQ=255)
                                          cAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctggct                                                                                                           <  1:2820191/150‑1 (MQ=255)
                                           aGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCACTGGCTTTAGGGCTGGCtt                                                                                                          >  1:3793955/1‑150 (MQ=255)
                                            gTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTg                                                                                                         <  2:1988769/150‑1 (MQ=255)
                                             tGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTgc                                                                                                        >  2:3155458/1‑150 (MQ=255)
                                             tGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTgc                                                                                                        >  1:1566844/1‑150 (MQ=255)
                                             tGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTgc                                                                                                        >  2:886918/1‑150 (MQ=255)
                                                  aGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttc                                                                                                   <  1:4359535/150‑1 (MQ=255)
                                                   gAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttct                                                                                                  >  1:1618014/1‑150 (MQ=255)
                                                       aTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGt                                                                                              >  1:1446355/1‑150 (MQ=255)
                                                        ttACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTg                                                                                             >  1:2230806/1‑150 (MQ=255)
                                                         tagaacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTgg                                                                                            <  2:3540172/147‑1 (MQ=255)
                                                         tACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTgg                                                                                            <  2:1254601/150‑1 (MQ=255)
                                                          acaacaAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTggg                                                                                           <  1:1118795/150‑1 (MQ=255)
                                                            aacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGcc                                                                                         >  1:3384195/1‑150 (MQ=255)
                                                                    gCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACt                                                                                 <  2:2636913/150‑1 (MQ=255)
                                                                      gATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTgg                                                                               <  1:2893828/150‑1 (MQ=255)
                                                                      gATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTgg                                                                               <  2:3576390/150‑1 (MQ=255)
                                                                      gATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTgg                                                                               <  1:1498715/150‑1 (MQ=255)
                                                                           cGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGACGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAg                                                                          >  1:2179113/1‑150 (MQ=255)
                                                                            gCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGc                                                                         >  1:3542151/1‑150 (MQ=255)
                                                                              aaaGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCct                                                                       >  2:3163639/1‑150 (MQ=255)
                                                                              aaaGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCct                                                                       <  2:1106533/150‑1 (MQ=255)
                                                                                  tttttCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTa                                                                   <  1:434498/150‑1 (MQ=255)
                                                                                   ttttCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTac                                                                  <  2:2498451/150‑1 (MQ=255)
                                                                                    tttCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTaca                                                                 <  1:406310/150‑1 (MQ=255)
                                                                                      tctcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACt                                                               <  1:1164430/150‑1 (MQ=255)
                                                                                       ctcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACtt                                                              >  1:4389137/1‑150 (MQ=255)
                                                                                         cAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGa                                                            >  2:3596330/1‑150 (MQ=255)
                                                                                          aaTGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGAc                                                           >  2:2474812/1‑150 (MQ=255)
                                                                                           aTGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGAcc                                                          <  1:1299070/150‑1 (MQ=255)
                                                                                               gTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCt                                                      <  1:2523488/150‑1 (MQ=255)
                                                                                               gTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCt                                                      <  1:4249014/150‑1 (MQ=255)
                                                                                                    ttCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaac                                                 <  2:380675/150‑1 (MQ=255)
                                                                                                       aGAATGGCGCACTCAACACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaac                                              <  2:30488/150‑1 (MQ=255)
                                                                                                       aGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaac                                              >  2:312474/1‑150 (MQ=255)
                                                                                                       aGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaac                                              >  1:2595635/1‑150 (MQ=255)
                                                                                                         aaTGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaacaa                                            >  1:2977438/1‑150 (MQ=255)
                                                                                                            ggCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAAt                                         <  2:3952063/150‑1 (MQ=255)
                                                                                                            ggCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAAt                                         <  2:170727/150‑1 (MQ=255)
                                                                                                              cgcACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATgg                                       >  1:86584/1‑150 (MQ=255)
                                                                                                              cgcACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATgg                                       <  1:988909/150‑1 (MQ=255)
                                                                                                                    cAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGc                                 >  1:2006147/1‑150 (MQ=255)
                                                                                                                     aaaaCTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCg                                <  2:327245/150‑1 (MQ=255)
                                                                                                                     aaaaCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCg                                <  2:560905/150‑1 (MQ=255)
                                                                                                                      aaaCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCga                               >  1:4380145/1‑150 (MQ=255)
                                                                                                                      aaaCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCga                               >  2:1634768/1‑150 (MQ=255)
                                                                                                                       aaCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgat                              <  2:2819018/150‑1 (MQ=255)
                                                                                                                         cTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatag                            <  1:2127108/150‑1 (MQ=255)
                                                                                                                          tACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgataga                           <  2:3899404/150‑1 (MQ=255)
                                                                                                                           aCAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatagat                          <  1:1015059/150‑1 (MQ=255)
                                                                                                                            cAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatagata                         >  2:2230079/1‑150 (MQ=255)
                                                                                                                             aaTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATAc                        <  1:3155458/150‑1 (MQ=255)
                                                                                                                             aaTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATAc                        <  1:3815289/150‑1 (MQ=255)
                                                                                                                                  cGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGcc                   <  2:2359637/150‑1 (MQ=255)
                                                                                                                                  cGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGcc                   <  2:4467718/150‑1 (MQ=255)
                                                                                                                                   gggATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGccc                  <  1:39069/150‑1 (MQ=255)
                                                                                                                                    ggATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGcccc                 <  2:4132633/150‑1 (MQ=255)
                                                                                                                                     gATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCa                <  2:2977438/150‑1 (MQ=255)
                                                                                                                                          tCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGcgcg           >  2:3640644/1‑150 (MQ=255)
                                                                                                                                          tCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGcgcg           >  2:4516995/1‑150 (MQ=255)
                                                                                                                                            gATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtt         >  1:2699307/1‑150 (MQ=255)
                                                                                                                                            gATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtt         >  1:2699241/1‑150 (MQ=255)
                                                                                                                                              tGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtttt       <  1:3163689/150‑1 (MQ=255)
                                                                                                                                              tGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtttt       >  1:794976/1‑150 (MQ=255)
                                                                                                                                                aaaaGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTtat     <  2:4319837/150‑1 (MQ=255)
                                                                                                                                                aaaaGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTtat     >  1:3278958/1‑150 (MQ=255)
                                                                                                                                                aaaaGTGTTTCTCAAAGGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAACATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTtat     <  2:3895031/150‑1 (MQ=255)
                                                                                                                                                   cGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTATAAc  <  2:1789432/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTATAAC  >  CP009273/1833583‑1833879

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: