Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,833,730 C→A 83.8% 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,7300CA83.8% 682.1 / 144.3 272R3S (CGT→AGT) ynjEmolybdopterin synthase sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (18/26);  new base A (126/102);  total (144/128)
Fisher's exact test for biased strand distribution p-value = 9.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 276 total aligned reads displayed.

TCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTAT  >  CP009273/1833584‑1833876
                                                                                                                                                  |                                                                                                                                                  
  cAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACgtgt                                                                                                                                                >  1:3957407/1‑149 (MQ=255)
  cAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACgtgt                                                                                                                                                >  1:1939089/1‑149 (MQ=255)
  cAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtgt                                                                                                                                                <  1:2359601/149‑1 (MQ=255)
   aGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGtt                                                                                                                                               >  2:1605205/1‑149 (MQ=255)
   aGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGtt                                                                                                                                               >  2:3597840/1‑149 (MQ=255)
    ggATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGttt                                                                                                                                              <  2:3659408/149‑1 (MQ=255)
     gATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTtc                                                                                                                                             <  2:3679912/149‑1 (MQ=255)
       tGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTtctc                                                                                                                                           >  1:4610492/1‑149 (MQ=255)
        gTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAACCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCa                                                                                                                                          <  1:4186989/149‑1 (MQ=255)
         ttCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCaa                                                                                                                                         >  2:565585/1‑149 (MQ=255)
          tCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCaaa                                                                                                                                        <  1:3207731/149‑1 (MQ=255)
            ctgctgATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATg                                                                                                                                      >  2:1497715/1‑149 (MQ=255)
                    aGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCt                                                                                                                              <  1:28641/149‑1 (MQ=255)
                       tctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGc                                                                                                                           >  1:2030947/1‑149 (MQ=255)
                         tgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCaa                                                                                                                         <  2:2454301/149‑1 (MQ=255)
                           ttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATg                                                                                                                       <  2:3428945/149‑1 (MQ=255)
                               ggATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCGGGCAATGGCtt                                                                                                                   <  1:1322155/149‑1 (MQ=255)
                                gATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                  >  2:4415219/1‑149 (MQ=255)
                                gATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                  <  2:1714174/149‑1 (MQ=255)
                                gATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCACACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCttt                                                                                                                  >  1:4104664/1‑149 (MQ=255)
                                 atatGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTa                                                                                                                 >  1:3824385/1‑149 (MQ=255)
                                  tatGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAg                                                                                                                >  1:31296/1‑149 (MQ=255)
                                  tatGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAg                                                                                                                <  2:1947111/149‑1 (MQ=255)
                                    tGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAggg                                                                                                              <  1:1513492/149‑1 (MQ=255)
                                     ggCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggc                                                                                                             <  2:4240797/149‑1 (MQ=255)
                                     ggCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggc                                                                                                             >  2:4653708/1‑149 (MQ=255)
                                     ggCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggc                                                                                                             <  1:3992629/149‑1 (MQ=255)
                                      gcgcAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggct                                                                                                            >  2:1285739/1‑149 (MQ=255)
                                       cgcAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctg                                                                                                           <  2:3498749/149‑1 (MQ=255)
                                          aGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctggct                                                                                                        <  2:202453/149‑1 (MQ=255)
                                          aGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctggct                                                                                                        <  1:2720500/149‑1 (MQ=255)
                                             ggTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTgc                                                                                                     <  1:4500349/149‑1 (MQ=255)
                                               tCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTgcgc                                                                                                   >  2:731603/1‑149 (MQ=255)
                                                cAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGct                                                                                                  <  1:4186925/149‑1 (MQ=255)
                                                 aGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGctt                                                                                                 >  1:3403117/1‑149 (MQ=255)
                                                   aaaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttct                                                                                               <  2:1627672/149‑1 (MQ=255)
                                                   aaaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttct                                                                                               >  1:3048069/1‑149 (MQ=255)
                                                    aaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              <  1:1392554/149‑1 (MQ=255)
                                                    aaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              >  1:4025080/1‑149 (MQ=255)
                                                     aaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCttcttc                                                                                             >  1:3226360/1‑149 (MQ=255)
                                                         acaacaAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTgg                                                                                         <  1:1040159/149‑1 (MQ=255)
                                                          caacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTggg                                                                                        <  2:3219582/149‑1 (MQ=255)
                                                          caacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTggg                                                                                        >  2:4493559/1‑149 (MQ=255)
                                                            acaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGcc                                                                                      >  1:324986/1‑149 (MQ=255)
                                                              aaaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGc                                                                                    >  1:1843501/1‑149 (MQ=255)
                                                              aaaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGc                                                                                    >  2:3623380/1‑149 (MQ=255)
                                                                   gCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAAc                                                                               <  2:4106490/149‑1 (MQ=255)
                                                                   gCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAAc                                                                               >  1:1199353/1‑149 (MQ=255)
                                                                    cGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACt                                                                              <  1:3278540/149‑1 (MQ=255)
                                                                       tAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGc                                                                           <  1:2596406/149‑1 (MQ=255)
                                                                       tAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGc                                                                           >  2:117208/1‑149 (MQ=255)
                                                                        aaCGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCg                                                                          <  2:2288086/149‑1 (MQ=255)
                                                                        aaCGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCg                                                                          <  2:2387877/149‑1 (MQ=255)
                                                                             aaaGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGcc                                                                     >  2:1893352/1‑149 (MQ=255)
                                                                             aaaGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGcc                                                                     <  2:3254976/149‑1 (MQ=255)
                                                                              aaGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCct                                                                    >  1:4446356/1‑149 (MQ=255)
                                                                               aGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCctc                                                                   >  2:1976387/1‑149 (MQ=255)
                                                                                 tttttCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCtt                                                                 >  1:3668621/1‑149 (MQ=255)
                                                                                    ttCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTaca                                                              <  2:2765683/149‑1 (MQ=255)
                                                                                     tctcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTacac                                                             >  1:4026049/1‑149 (MQ=255)
                                                                                     tctcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTacac                                                             <  1:1673600/149‑1 (MQ=255)
                                                                                       tcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACtt                                                           >  1:1634717/1‑149 (MQ=255)
                                                                                       tcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACtt                                                           <  1:2913842/149‑1 (MQ=255)
                                                                                       tcAATGCGTCAGTTCAGAATGGCGCACTCAAAACGACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACtt                                                           <  1:1007668/149‑1 (MQ=255)
                                                                                         aaTGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGa                                                         <  2:4270030/149‑1 (MQ=255)
                                                                                                cAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCtt                                                  <  1:3615790/149‑1 (MQ=255)
                                                                                                cAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCtt                                                  >  2:4226247/1‑149 (MQ=255)
                                                                                                   ttCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaa                                               <  1:673792/149‑1 (MQ=255)
                                                                                                    tCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaac                                              <  1:1438835/149‑1 (MQ=255)
                                                                                                      aGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaa                                            <  2:662968/149‑1 (MQ=255)
                                                                                                      aGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaa                                            >  1:2860257/1‑149 (MQ=255)
                                                                                                      aGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaa                                            <  2:1103210/149‑1 (MQ=255)
                                                                                                          tGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACaaa                                        >  1:4378263/1‑149 (MQ=255)
                                                                                                           ggCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACaaaa                                       <  2:3146888/149‑1 (MQ=255)
                                                                                                                 cTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCaa                                 <  2:4427894/149‑1 (MQ=255)
                                                                                                                 cTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCaa                                 <  1:4099106/149‑1 (MQ=255)
                                                                                                                   cAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAg                               <  2:2604456/149‑1 (MQ=255)
                                                                                                                   cAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAg                               <  2:4206961/149‑1 (MQ=255)
                                                                                                                    aaaaCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGc                              >  1:4294379/1‑149 (MQ=255)
                                                                                                                       aCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgat                           <  1:4229568/149‑1 (MQ=255)
                                                                                                                          aCAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgataga                        >  2:2955945/1‑149 (MQ=255)
                                                                                                                          aCAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgat                           >  1:2394330/1‑146 (MQ=255)
                                                                                                                          aCAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATAGCAAAGCgataga                        <  1:593580/149‑1 (MQ=255)
                                                                                                                            aaTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatagata                      >  2:2866034/1‑149 (MQ=255)
                                                                                                                            aaTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatagata                      >  2:2512520/1‑149 (MQ=255)
                                                                                                                            aaTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatagata                      >  2:2976462/1‑149 (MQ=255)
                                                                                                                             aTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATAc                     <  2:869039/149‑1 (MQ=255)
                                                                                                                               gTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTc                   <  1:2088695/149‑1 (MQ=255)
                                                                                                                                tCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCg                  >  1:3067670/1‑149 (MQ=255)
                                                                                                                                 cGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGc                 <  2:580769/149‑1 (MQ=255)
                                                                                                                                    gATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGcccc              >  1:3333996/1‑149 (MQ=255)
                                                                                                                                    gATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGcccc              >  1:3694761/1‑149 (MQ=255)
                                                                                                                                      ttttCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAg            >  1:1464731/1‑149 (MQ=255)
                                                                                                                                      ttttCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAg            <  2:131897/149‑1 (MQ=255)
                                                                                                                                       tttCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAgc           <  2:1453687/149‑1 (MQ=255)
                                                                                                                                           gATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGt       >  1:787103/1‑149 (MQ=255)
                                                                                                                                            aTGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtt      >  2:422992/1‑149 (MQ=255)
                                                                                                                                              gAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtttt    >  2:4511987/1‑149 (MQ=255)
                                                                                                                                               aaaaGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTta   >  2:4524181/1‑149 (MQ=255)
                                                                                                                                               aaaCGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTta   <  1:1605205/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTAT  >  CP009273/1833584‑1833876

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: