Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,833,730 C→A 100% 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,7300CA98.0% 1031.6 / 9.3 308R3S (CGT→AGT) ynjEmolybdopterin synthase sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base A (145/156);  total (147/161)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 311 total aligned reads displayed.

CTCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTAT  >  CP009273/1833583‑1833876
                                                                                                                                                   |                                                                                                                                                  
ctcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgt                                                                                                                                                  <  2:2804608/150‑1 (MQ=255)
 tcCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtg                                                                                                                                                 <  2:2563977/150‑1 (MQ=255)
  ccAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAgtgt                                                                                                                                                >  2:2102806/1‑150 (MQ=255)
  ccAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTACTCAATTCGTCAGTTCAGAATGTCGCACTCAAAACTACAATGTCGGGATTTTCTATGAAAAgtgt                                                                                                                                                >  2:3294846/1‑150 (MQ=255)
    aGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGttt                                                                                                                                              <  1:628176/150‑1 (MQ=255)
     ggATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTtc                                                                                                                                             <  2:3493181/150‑1 (MQ=255)
     ggATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTtc                                                                                                                                             <  1:3257306/150‑1 (MQ=255)
      gATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTtct                                                                                                                                            <  1:2498647/150‑1 (MQ=255)
       aTGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTtctc                                                                                                                                           >  2:1039241/1‑150 (MQ=255)
       aTGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTtctc                                                                                                                                           >  1:3197475/1‑150 (MQ=255)
         gTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCaa                                                                                                                                         >  2:1100529/1‑150 (MQ=255)
          ttCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCaaa                                                                                                                                        <  2:442900/150‑1 (MQ=255)
          ttCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCaaa                                                                                                                                        <  2:191768/150‑1 (MQ=255)
           tCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAAt                                                                                                                                       <  1:797211/150‑1 (MQ=255)
            ccTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATg                                                                                                                                      <  2:1777420/150‑1 (MQ=255)
                    tAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCt                                                                                                                              <  1:4631263/150‑1 (MQ=255)
                     aGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTg                                                                                                                             <  2:4377937/150‑1 (MQ=255)
                       gtctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGc                                                                                                                           <  2:423392/149‑1 (MQ=255)
                        tctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACCAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCa                                                                                                                          <  2:4175502/150‑1 (MQ=255)
                         ctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCaa                                                                                                                         <  2:2654420/150‑1 (MQ=255)
                         ctgttctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCaa                                                                                                                         >  2:455216/1‑150 (MQ=255)
                             tctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGc                                                                                                                     >  2:4183957/1‑150 (MQ=255)
                              ctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCt                                                                                                                    >  2:3749861/1‑150 (MQ=255)
                              ctgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCt                                                                                                                    <  1:1967972/150‑1 (MQ=255)
                               tgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACACTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCtt                                                                                                                   <  1:1279121/150‑1 (MQ=255)
                               tgGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCtt                                                                                                                   <  2:4140923/150‑1 (MQ=255)
                                 gATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTa                                                                                                                 <  1:1569054/150‑1 (MQ=255)
                                  atatGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAg                                                                                                                >  2:653179/1‑150 (MQ=255)
                                   tatGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAgg                                                                                                               <  1:1590367/150‑1 (MQ=255)
                                     tGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggc                                                                                                             <  1:4542167/150‑1 (MQ=255)
                                     tGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTAAAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggc                                                                                                             >  1:3903847/1‑150 (MQ=255)
                                     tGGCGCAGTGGTCAGAAAATTACAACAAACNGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggc                                                                                                             >  1:4276155/1‑150 (MQ=255)
                                        cgcAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctgg                                                                                                          >  1:3787343/1‑150 (MQ=255)
                                        cgcAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctgg                                                                                                          >  2:267110/1‑150 (MQ=255)
                                          cAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGggctggct                                                                                                        >  1:4418561/1‑150 (MQ=255)
                                           aGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCtt                                                                                                       <  2:3854139/150‑1 (MQ=255)
                                              ggTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTgcg                                                                                                    <  1:2165929/150‑1 (MQ=255)
                                                 cAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGctt                                                                                                 >  2:2607234/1‑150 (MQ=255)
                                                  aGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttc                                                                                                >  2:1440080/1‑150 (MQ=255)
                                                    aaaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              >  1:3420185/1‑150 (MQ=255)
                                                    aaaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              <  2:1686521/150‑1 (MQ=255)
                                                    aaaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGcttctt                                                                                              >  1:2959558/1‑150 (MQ=255)
                                                     aaaTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCttcttc                                                                                             <  1:2269567/150‑1 (MQ=255)
                                                       aTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGt                                                                                           >  2:2184558/1‑150 (MQ=255)
                                                       aTTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGt                                                                                           >  2:75219/1‑150 (MQ=255)
                                                        ttACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTg                                                                                          >  2:288651/1‑150 (MQ=255)
                                                         tACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTgg                                                                                         >  1:3421226/1‑150 (MQ=255)
                                                            aacaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGcc                                                                                      <  1:2371137/150‑1 (MQ=255)
                                                             acaaACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCg                                                                                     >  1:3727372/1‑150 (MQ=255)
                                                                aaCTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTg                                                                                  >  1:1530411/1‑150 (MQ=255)
                                                                    gCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACt                                                                              >  2:3146802/1‑150 (MQ=255)
                                                                     cGATAACGCATAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACACTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTg                                                                             <  1:18082/150‑1 (MQ=255)
                                                                       aTAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACACTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGc                                                                           <  1:1695683/150‑1 (MQ=255)
                                                                       aTAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGc                                                                           <  1:2835365/150‑1 (MQ=255)
                                                                        tAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCg                                                                          >  1:1544271/1‑150 (MQ=255)
                                                                          aCGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGaa                                                                        >  1:2516295/1‑150 (MQ=255)
                                                                            gCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGc                                                                      >  1:2029454/1‑150 (MQ=255)
                                                                            gCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGc                                                                      <  2:2411359/150‑1 (MQ=255)
                                                                             cAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGcc                                                                     <  1:2575137/150‑1 (MQ=255)
                                                                             cAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGcc                                                                     >  1:101646/1‑150 (MQ=255)
                                                                              aaaGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCct                                                                    >  2:272256/1‑150 (MQ=255)
                                                                               aaGTTTTTCTAAATGCGTCAGTTCAGAATGGCGCAATCAAAACTACAATGTAGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGATGGCTTGAGCTTATTCGTGGTCCGCTGAAATGGCGAAGCctc                                                                   <  2:1155430/150‑1 (MQ=255)
                                                                                aGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCctct                                                                  >  1:4067010/1‑150 (MQ=255)
                                                                                  tttttCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTa                                                                <  2:187156/150‑1 (MQ=255)
                                                                                    tttCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTaca                                                              <  1:3780952/150‑1 (MQ=255)
                                                                                        tcAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTg                                                          >  1:3545843/1‑150 (MQ=255)
                                                                                           aTGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGAcc                                                       <  1:1434877/150‑1 (MQ=255)
                                                                                            tGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCa                                                      <  2:1086826/150‑1 (MQ=255)
                                                                                              cGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGc                                                    >  1:1149009/1‑150 (MQ=255)
                                                                                               gTCAGTTCAGAATGGCGCACTCAAAACTACGATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCt                                                   >  2:1608484/1‑150 (MQ=255)
                                                                                               gTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCt                                                   <  1:683929/150‑1 (MQ=255)
                                                                                               gTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCt                                                   >  2:3603836/1‑150 (MQ=255)
                                                                                                 cAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTc                                                 >  1:989322/1‑150 (MQ=255)
                                                                                                  aGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTca                                                >  2:3070426/1‑150 (MQ=255)
                                                                                                     tCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaaca                                             >  1:4484040/1‑150 (MQ=255)
                                                                                                        gAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaaca                                          >  1:2640654/1‑150 (MQ=255)
                                                                                                         aaTGGCGCACTCAAAACTACAATGTCGGGANTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTcaacaacaa                                         >  1:4283810/1‑150 (MQ=255)
                                                                                                          aTGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACaaa                                        <  2:3947202/150‑1 (MQ=255)
                                                                                                             gcgcACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATg                                     <  2:335153/150‑1 (MQ=255)
                                                                                                              cgcACTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGATTCAACAACAAAATgg                                    <  2:50967/150‑1 (MQ=255)
                                                                                                                  cTCAAAACTACAATGTCGGTATTTTCGATGATAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAATATGGCaaa                                <  1:3096607/150‑1 (MQ=255)
                                                                                                                  cTCAAAACTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCaaa                                >  1:3680508/1‑150 (MQ=255)
                                                                                                                       aaCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgat                           <  2:322707/150‑1 (MQ=255)
                                                                                                                        aCTACAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgata                          <  1:1270516/150‑1 (MQ=255)
                                                                                                                           aCAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatagat                       >  1:1031655/1‑150 (MQ=255)
                                                                                                                            cAATGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCgatagata                      <  1:2244073/150‑1 (MQ=255)
                                                                                                                             aaTGTCGGGTTTTTCGATGAAAAGTGTTTATCAAATGACCGCGCTGGCAATGGCTTTATGGCTGGCTTGCGATTCTTCGTGGGCCGCTGTACTGGCGAAGCCTCTTAAACTTGACCAGCTTCAACAACAATATGGCAAAGCGATAGATAc                     <  2:1311263/150‑1 (MQ=255)
                                                                                                                             aaTGTCGGGATTTTCGATGAAAAGTGTTTATCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATAc                     <  2:3894818/150‑1 (MQ=255)
                                                                                                                              aTGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACt                    <  2:3934691/150‑1 (MQ=255)
                                                                                                                               tGTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCATCTTCAACAACAAAATGGCAAAGCGATAGATACTc                   <  2:4296297/150‑1 (MQ=255)
                                                                                                                                gTCGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCg                  >  2:2338113/1‑150 (MQ=255)
                                                                                                                                  cGGGATTTTCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGcc                >  2:3183355/1‑150 (MQ=255)
                                                                                                                                        tttCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAgcg          <  1:4612872/150‑1 (MQ=255)
                                                                                                                                          tCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGcgcg        <  1:2688416/150‑1 (MQ=255)
                                                                                                                                          tCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGcgcg        >  2:4626221/1‑150 (MQ=255)
                                                                                                                                          tCGATGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGcgcg        <  2:677758/150‑1 (MQ=255)
                                                                                                                                             aTGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGttt     >  1:3576617/1‑150 (MQ=255)
                                                                                                                                              tGAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGtttt    <  2:128136/150‑1 (MQ=255)
                                                                                                                                               gAAAAGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTta   >  1:2785103/1‑150 (MQ=255)
                                                                                                                                                aaaCGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTtat  <  2:2284768/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCCAGGATGTTCCTGCTGATAGTTCTGTTCTGGATATGGCGCAGTGGTCAGAAAATTACAACAAACTGCGATAACGCAAAGTTTTTCTCAATGCGTCAGTTCAGAATGGCGCACTCAAAACTACAATGTCGGGATTTTCGATGAAACGTGTTTCTCAAATGACCGCGCTGGCAATGGCTTTAGGGCTGGCTTGCGCTTCTTCGTGGGCCGCTGAACTGGCGAAGCCTCTTACACTTGACCAGCTTCAACAACAAAATGGCAAAGCGATAGATACTCGCCCCAGCGCGTTTTAT  >  CP009273/1833583‑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: