Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 411,556 C→A 24.6% D285Y (GAT→TAT)  sbcD ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273411,5560CA24.6% 695.6 / 170.3 296D285Y (GAT→TAT) sbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (124/99);  new base A (36/37);  total (160/136)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 300 total aligned reads displayed.

TTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATGGGGGAGCCGCAATAGCGA  >  CP009273/411407‑411704
                                                                                                                                                     |                                                                                                                                                    
ttCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATc                                                                                                                                                      <  1:3698213/150‑1 (MQ=255)
ttCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATc                                                                                                                                                      <  1:2924978/150‑1 (MQ=255)
    aCAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCt                                                                                                                                                  >  2:119368/1‑150 (MQ=255)
    aCAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCt                                                                                                                                                  >  1:505841/1‑150 (MQ=255)
    aCAGGCAAAGATTCGGTTAATGCCTGGATTTTGCGATGAATAACATGCAGATACACATCAGTAGTGATTTAGATATCCAGCCAGACAGGTGTCTCCTGCGATACATCGCTCCACTGAACCAGCTGAGAGGTAATCGACGCCAGATCGCCt                                                                                                                                                  >  1:5180367/1‑150 (MQ=255)
      aGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCttt                                                                                                                                                <  1:315138/150‑1 (MQ=255)
        gCAATGATTCGGTTAATGCCTGGATTTTGCGTTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCa                                                                                                                                              >  1:1697658/1‑150 (MQ=255)
        gCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCa                                                                                                                                              >  1:2771467/1‑150 (MQ=255)
        gCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCa                                                                                                                                              <  2:3677034/150‑1 (MQ=255)
           aTGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCa                                                                                                                                           >  1:3708834/1‑150 (MQ=255)
             gATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACt                                                                                                                                         <  2:1854793/150‑1 (MQ=255)
               ttCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGc                                                                                                                                       >  1:2604937/1‑150 (MQ=255)
                 cGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCa                                                                                                                                     <  1:3413455/150‑1 (MQ=255)
                   gTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCCTg                                                                                                                                   >  2:3379493/1‑150 (MQ=255)
                   gTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATg                                                                                                                                   >  1:464209/1‑150 (MQ=255)
                   gTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATg                                                                                                                                   >  1:5263457/1‑150 (MQ=255)
                        tGCCTGTATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTg                                                                                                                              <  1:4372512/150‑1 (MQ=255)
                          ccTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCg                                                                                                                            >  1:1032735/1‑150 (MQ=255)
                          ccTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCg                                                                                                                            >  2:2218334/1‑150 (MQ=255)
                          ccTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCg                                                                                                                            >  2:453832/1‑150 (MQ=255)
                              gATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTAc                                                                                                                        <  2:1628946/150‑1 (MQ=255)
                                 tttGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCgg                                                                                                                     <  1:3544560/150‑1 (MQ=255)
                                       cTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGtt                                                                                                               >  2:2733927/1‑150 (MQ=255)
                                        tGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTc                                                                                                              >  2:2632416/1‑150 (MQ=255)
                                        tGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTc                                                                                                              <  1:1970558/150‑1 (MQ=255)
                                         gAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCa                                                                                                             <  2:375946/150‑1 (MQ=255)
                                          aaTATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAg                                                                                                            <  1:979986/150‑1 (MQ=255)
                                           atatCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAgg                                                                                                           >  1:1873191/1‑150 (MQ=255)
                                           atatCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAgg                                                                                                           >  1:3456490/1‑150 (MQ=255)
                                            tatCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGt                                                                                                          <  1:4581001/150‑1 (MQ=255)
                                             atCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGtt                                                                                                         >  1:2703720/1‑150 (MQ=255)
                                                  gCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCa                                                                                                    >  1:2461079/1‑150 (MQ=255)
                                                  gCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCa                                                                                                    >  1:5025144/1‑150 (MQ=255)
                                                  gCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCa                                                                                                    <  2:3978761/150‑1 (MQ=255)
                                                  gCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCa                                                                                                    <  2:1472256/150‑1 (MQ=255)
                                                  gCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCa                                                                                                    <  2:5042883/150‑1 (MQ=255)
                                                    aGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACg                                                                                                  >  1:3323246/1‑150 (MQ=255)
                                                    aGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACg                                                                                                  <  2:3556650/150‑1 (MQ=255)
                                                      aTACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGct                                                                                                >  2:5159215/1‑150 (MQ=255)
                                                        aCTCATCAGTAGTGATTTCGATCTCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGctct                                                                                              <  2:4007325/150‑1 (MQ=255)
                                                        aCTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGctct                                                                                              >  1:2855513/1‑150 (MQ=255)
                                                         cTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGctctc                                                                                             >  2:3471191/1‑150 (MQ=255)
                                                          tcatcaGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGctctct                                                                                            <  1:3835217/150‑1 (MQ=255)
                                                           catcaGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTa                                                                                           >  2:4730913/1‑150 (MQ=255)
                                                            atcaGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGAACGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTaa                                                                                          >  1:408590/1‑150 (MQ=255)
                                                             tcaGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAAt                                                                                         >  2:674899/1‑150 (MQ=255)
                                                                gtagtGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTg                                                                                      >  1:2149714/1‑150 (MQ=255)
                                                                   gtGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCg                                                                                   >  1:4091996/1‑150 (MQ=255)
                                                                    tGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCACCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGt                                                                                  >  1:4429893/1‑150 (MQ=255)
                                                                      aTTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGttt                                                                                <  1:4030359/150‑1 (MQ=255)
                                                                      aTTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGttt                                                                                >  1:23971/1‑150 (MQ=255)
                                                                       tttCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTg                                                                               <  1:902335/150‑1 (MQ=255)
                                                                        ttCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGa                                                                              <  2:4583367/150‑1 (MQ=255)
                                                                         tCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGaa                                                                             >  1:3088627/1‑150 (MQ=255)
                                                                         tCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGaa                                                                             <  1:1446544/150‑1 (MQ=255)
                                                                          cGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGaaa                                                                            <  2:4837212/150‑1 (MQ=255)
                                                                           gATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGaaaa                                                                           >  2:1987007/1‑150 (MQ=255)
                                                                           gATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGaaaa                                                                           >  2:3014609/1‑150 (MQ=255)
                                                                           gATATCCAGCCAGACAGGTGGCACCTGCGATACATCGCGCCACTGATCCAGCTGTGCGGTAATCGACGCCAGAACGCCTTTCAGCACTGACATGGGTTGCGTAACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGACGTTTGaaaa                                                                           >  2:3123000/1‑150 (MQ=255)
                                                                            atatCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGCGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAAt                                                                          <  1:4129657/150‑1 (MQ=255)
                                                                             tatCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATg                                                                         <  1:599830/150‑1 (MQ=255)
                                                                               tCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTc                                                                       <  2:2909732/150‑1 (MQ=255)
                                                                                   gccagACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCa                                                                   <  2:2361322/150‑1 (MQ=255)
                                                                                    ccagACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAg                                                                  >  2:718808/1‑150 (MQ=255)
                                                                                      agACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCAGGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGAt                                                                >  2:4715401/1‑150 (MQ=255)
                                                                                             tGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACata                                                         <  2:3902282/150‑1 (MQ=255)
                                                                                             tGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACata                                                         >  1:2555569/1‑150 (MQ=255)
                                                                                               gCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACatagt                                                       <  2:2113044/150‑3 (MQ=255)
                                                                                               gCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATAtt                                                       >  1:4728838/1‑150 (MQ=255)
                                                                                               gCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATAtt                                                       >  1:4689489/1‑150 (MQ=255)
                                                                                               gCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATAtt                                                       <  1:3933770/150‑1 (MQ=255)
                                                                                                  ccTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTAc                                                    >  2:1972247/1‑150 (MQ=255)
                                                                                                  ccTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGGTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTAc                                                    >  2:5272838/1‑150 (MQ=255)
                                                                                                    tGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTActc                                                  <  1:1596962/150‑1 (MQ=255)
                                                                                                     gCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTActct                                                 <  1:119368/150‑1 (MQ=255)
                                                                                                        aTACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTAc                                              <  2:3445397/150‑1 (MQ=255)
                                                                                                           cATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGc                                           >  2:5068965/1‑150 (MQ=255)
                                                                                                           cATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGc                                           >  2:798082/1‑150 (MQ=255)
                                                                                                              cgcgCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCAtt                                        >  2:884841/1‑150 (MQ=255)
                                                                                                                cgcCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATtca                                      >  2:3775553/1‑150 (MQ=255)
                                                                                                                  ccACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATtcatc                                    <  2:4425227/150‑1 (MQ=255)
                                                                                                                   cACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATtcatca                                   <  2:3291228/150‑1 (MQ=255)
                                                                                                                   cACTGTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATtcatca                                   <  2:3795878/150‑1 (MQ=255)
                                                                                                                    aCTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCaa                                  <  1:4948895/150‑1 (MQ=255)
                                                                                                                    aCTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCaa                                  >  1:4991534/1‑150 (MQ=255)
                                                                                                                       gTTCCAGCTGTGCGGTAATCGACGCCAGATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGGTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAAc                               >  2:4620171/1‑150 (MQ=255)
                                                                                                                           cAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAg                           <  2:3003672/150‑1 (MQ=255)
                                                                                                                             gCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTg                         <  2:1758203/150‑1 (MQ=255)
                                                                                                                             gCTGTGCGGTAATCGACGCCAGATAGCCTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTgg                        >  1:2280789/1‑150 (MQ=255)
                                                                                                                                 tgCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAAt                     <  2:1060737/150‑1 (MQ=255)
                                                                                                                                 tgCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAAt                     <  1:4492304/150‑1 (MQ=255)
                                                                                                                                 tgCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAAt                     <  2:4159077/150‑1 (MQ=255)
                                                                                                                                   cGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATgg                   <  2:4971503/150‑1 (MQ=255)
                                                                                                                                     gTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATgggg                 <  1:1166634/150‑1 (MQ=255)
                                                                                                                                        aTCGAAGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCCGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATGGGGGAg              <  1:3492341/150‑1 (MQ=255)
                                                                                                                                          cGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATGGGGGAGcc            <  1:1102745/150‑1 (MQ=255)
                                                                                                                                            aCGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATGGGGGAGCCGc          <  1:674899/150‑1 (MQ=255)
                                                                                                                                                 aGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATGGGGGAGCCGCAATAg     <  1:4889002/150‑1 (MQ=255)
                                                                                                                                                  gATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATGGGGGAGCCGCAATAGc    >  2:2466462/1‑150 (MQ=255)
                                                                                                                                                  gATAGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATGGGGGAGCCGCAATAGc    >  2:1364656/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCGACAGGCAATGATTCGGTTAATGCCTGGATTTTGCGCTGAATATCATGCAGATACTCATCAGTAGTGATTTCGATATCCAGCCAGACAGGTGGCTCCTGCGATACATCGCGCCACTGTTCCAGCTGTGCGGTAATCGACGCCAGATCGCCTTTCAGCACTGCCATGGGTTGCGTTACCGGTACGTTCAGGTTTTCCACGCTCTCTAATTTGCCGTTTGAAAATGTCACCAGATGGACATATTTACTCTTACCGCATTCATCAAAACTCAGTGGAATGGGGGAGCCGCAATAGCGA  >  CP009273/411407‑411704

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: