Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,364,366 G→A 80.5% intergenic (+106/‑107) pspE → / → ycjM thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,364,3660GA80.5% 605.9 / 135.6 274intergenic (+106/‑107)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base G (30/22);  new base A (101/116);  total (132/142)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

CGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACG  >  CP009273/1364218‑1364512
                                                                                                                                                      |                                                                                                                                                  
        tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTGAGTCACACCTTATTTAAATACCCAGACATTATTTTCGTCAGCTCAAAAAAACGCTTTTTACCTAGGaaaaaaa                                                                                                                                           >  1:3661351/1‑149 (MQ=255)
         gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTGTGTCACACCGTATTTAAATCCCCCTAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTCAAAAATg                                                                                                                                          >  1:3728207/1‑150 (MQ=255)
            aGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGAtt                                                                                                                                       <  1:4732876/150‑1 (MQ=255)
              aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCTAGTAAAAAATGATTTc                                                                                                                                     >  1:864812/1‑150 (MQ=255)
               cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTcc                                                                                                                                    <  2:1986783/150‑1 (MQ=255)
               cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATGTCACACCGTATTTACATCCCCAAAGATTATTTGCGGCAGCTCCCAAATACGCTTTTTACCGAGTTAAAAATGATTTcc                                                                                                                                    >  1:955114/1‑150 (MQ=255)
                aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt                                                                                                                                   <  1:2737214/150‑1 (MQ=255)
                aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTAAATACCCATAGATTATTTGCGTCAGCTCACAAATAAGCTTTTTCCCGGGTTAAAAATGATTTCCt                                                                                                                                   >  1:2552870/1‑150 (MQ=255)
                  cGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTGTGTGACACCTTATTTACATACCCATAGATTATTTGCGTCAGCGCACAAATAAGCTTTTTACCTAGTAAAAAATTATTTACTGc                                                                                                                                 >  1:233907/1‑150 (MQ=255)
                    cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGt                                                                                                                               <  1:4210036/150‑1 (MQ=255)
                       tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATGTCACACCATATTTAAATCCCCAAATATTATTTTCGTAAGCACACAAAAAAGCATTTTACCAAGTAAAAAATTATTTCCAGCGTgtc                                                                                                                            >  1:99874/1‑148 (MQ=255)
                       tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTGATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCGAGTAAAAAATGATTTCATGCGTGAc                                                                                                                            >  2:55303/1‑150 (MQ=255)
                        gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACt                                                                                                                           <  2:2254473/150‑1 (MQ=255)
                        gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATGTCACACCGTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACt                                                                                                                           >  2:4464445/1‑150 (MQ=255)
                               gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATAAGCTTTTTCCCTAGTAAAAAATGATTTACTGGGTGACTAAAAccc                                                                                                                    >  1:2855612/1‑150 (MQ=255)
                               gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACAAACCCAAAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAAccc                                                                                                                    >  1:4691334/1‑150 (MQ=255)
                               gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTAAATACCCCTATTTTATTTGCCTCAGCTCACAAATAAGCTTTTTCCCTGGTAAAAAATTATTTACAGCGGGTCTAAAAacc                                                                                                                    >  2:221704/1‑150 (MQ=255)
                                 cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTCTGTCAAACCTTATTTAAATACCCAAATATTATTTGCGTCCGCTCACAAATAAGCATTTTCCCGAGTAAAAAATGATTTACTGCGTGTCTaaaaaac                                                                                                                    >  1:5059050/1‑147 (MQ=255)
                                 cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                                                                                                  <  2:273505/150‑1 (MQ=255)
                                 cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                                                                                                  <  2:3133237/150‑1 (MQ=255)
                                 cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                                                                                                  <  2:309463/150‑1 (MQ=255)
                                   aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt                                                                                                                <  1:5150006/150‑1 (MQ=255)
                                   aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt                                                                                                                <  2:83616/150‑1 (MQ=255)
                                   aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTAAAACCTTATTTACATCCCCATAGATTATTTTCGTCAGCTCACAAAAACGCTTTATCCCAAGTAAAAAATTATTTCCTGCGTGTCGAAAAACCAtgt                                                                                                                >  1:2917870/1‑150 (MQ=255)
                                   aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATGTAACACCATATTTAAAACCCCAAATATTATTTTCGTCACCACACAAAAAACCTTTTTACCAAGTAAAAAAATAATTACGGCGTGACTAAAAACCAtgt                                                                                                                >  1:1776069/1‑150 (MQ=255)
                                     ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATGTCACACCGTATTTACAAACCCAAATAATAATTTCGTCAGCAAACAAAAACCCTTTTTCCCAAGTAAAAAAAGAATTAATGCCTGACAAAAAaccaatttc                                                                                                              >  1:1957553/1‑144 (MQ=255)
                                       ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTATTAAAAAATGATTTCCTGCGTGACTAAAAACCTTGTGCTc                                                                                                            >  2:3846474/1‑150 (MQ=255)
                                       ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATGTCACACCTTATTTACATCCCCAAAGATTATTTGCGTCAGCTCACAAATAAGCTTTTTACCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTc                                                                                                            >  2:1440142/1‑150 (MQ=255)
                                        tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                                                                                                           <  1:2775205/150‑1 (MQ=255)
                                        tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                                                                                                           >  2:4204646/1‑150 (MQ=255)
                                          aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                                                                                                         <  2:4581805/150‑1 (MQ=255)
                                            aGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATTTCACACCATATTTACAAACCCAAATAATATTTGCGTCAGCACACAAAAAAGCTTTTTAACAAGTAAAAAAAGATTTACAGCGTGTCTAAAAACCATGTGGTCAAttt                                                                                                       >  1:1854200/1‑149 (MQ=255)
                                                   aCCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATACCCATATATTATTTGCGTCAGCTCACAAATACGCTTTTTACCTGGTAAAAAATGATTTACTGCGGGACTAAAACCCTTGTGCGCAAATGACAGttt                                                                                                >  2:3462862/1‑150 (MQ=255)
                                                     cATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt                                                                                              <  2:1152344/150‑1 (MQ=255)
                                                       tGAGCGGCGCTTATGCGCCGTTTTTTTTTTATGTCACACCTTATTTAAATCCCCAAATATTATTTTCGTCAGCTCACAAATAACCTTTTTACCAATTAAAAAAATATTTACTGCGTGACTAAAAACCCTGTTCTCAATTGAAACTTTAtt                                                                                             >  2:3661003/1‑150 (MQ=255)
                                                            ggcgCTTATGCGCCG‑‑TTTTTTTTTATTTAAAACCATATTTACAAACCCAAAGATTATTTTCGTCAGCTAACAAAAAACCTTTTTACCAATTAAAAAAATATTTAATGAGTGACTAAAAAACATTTTGTAAATTTACAGTTTAttttattc                                                                                       >  1:4858359/1‑146 (MQ=255)
                                                              cgcTTATGCGCCGTTTTTTTTTTATTTCACACCTTATTTACAACCCCAAATATTATTTTCGTCAGCTCACAAAAAAGCATTTTAACAAGTAAAAAATTAATTACTGCGTGACTAAAACCCTTTTGCTCAAATTACAGTTTAttttaatcg                                                                                      >  1:786271/1‑145 (MQ=255)
                                                               gcTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGa                                                                                    <  1:4994171/150‑1 (MQ=255)
                                                                  tATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGGAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCCATTGGCAGTTTTTTTTCTGCGGagta                                                                                 >  2:2922115/1‑150 (MQ=255)
                                                                  tATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagta                                                                                 <  2:3615004/150‑1 (MQ=255)
                                                                  tATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTGTCCCTAGTAAAAAATGATCTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagta                                                                                 <  2:679967/150‑1 (MQ=255)
                                                                   aTGCGCCG‑‑TTTTTTTTTCTTTCAAACCTTATTTACATCCCCATAGATTATTTTCGTCAGCTCACAAATAAGCTTTTTACCAAGTAAAAAATTATTTCCTGCGTGACTAAAAACCATGTGCTCAATTGACAGTTTATTTTATGCGGagtag                                                                                >  2:2401131/1‑150 (MQ=255)
                                                                      cgcCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGtct                                                                             <  2:1875025/150‑1 (MQ=255)
                                                                      cgcCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACAAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTATtct                                                                             >  2:1862764/1‑150 (MQ=255)
                                                                         cG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGtctctc                                                                          <  1:4296774/150‑1 (MQ=255)
                                                                               tttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTTGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGttt                                                                      <  1:1301150/150‑1 (MQ=255)
                                                                               tttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGttt                                                                      <  1:2122803/150‑1 (MQ=255)
                                                                                 tttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCa                                                                    <  1:2719262/150‑1 (MQ=255)
                                                                                 tttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCa                                                                    <  1:564966/150‑1 (MQ=255)
                                                                                  ttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCTCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCAt                                                                   <  1:1929752/150‑1 (MQ=255)
                                                                                  ttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCAt                                                                   >  2:3631212/1‑150 (MQ=255)
                                                                                  ttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCAt                                                                   >  2:3631268/1‑150 (MQ=255)
                                                                                    ttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATgg                                                                 >  2:4341407/1‑150 (MQ=255)
                                                                                    ttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATgg                                                                 <  2:1609722/150‑1 (MQ=255)
                                                                                      cTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGa                                                               <  1:4535167/150‑1 (MQ=255)
                                                                                      cTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGa                                                               <  1:4104261/150‑1 (MQ=255)
                                                                                       tGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAc                                                              <  2:5150582/150‑1 (MQ=255)
                                                                                       tGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAc                                                              >  1:4893691/1‑150 (MQ=255)
                                                                                        gTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAcc                                                             >  2:1737546/1‑150 (MQ=255)
                                                                                        gTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAcc                                                             <  2:1679902/150‑1 (MQ=255)
                                                                                        gTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAcc                                                             <  2:3957209/150‑1 (MQ=255)
                                                                                        gTCACACCTTATTTACATCCACATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCAATGGTAAAAAATGATTTCCTGCGTGAATAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAcc                                                             <  1:2472999/150‑1 (MQ=255)
                                                                                              ccTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAcc                                                       >  2:781457/1‑150 (MQ=255)
                                                                                                    ttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGaa                                                 <  2:3580855/150‑1 (MQ=255)
                                                                                                     tACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGaaa                                                >  1:1347597/1‑150 (MQ=255)
                                                                                                      aCATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGTAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAg                                               >  1:261106/1‑150 (MQ=255)
                                                                                                      aCATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAg                                               <  2:4624997/150‑1 (MQ=255)
                                                                                                      aCATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAg                                               >  2:2524576/1‑150 (MQ=255)
                                                                                                       cATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAgg                                              >  2:3625177/1‑150 (MQ=255)
                                                                                                          ccccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCaa                                           <  2:302550/150‑1 (MQ=255)
                                                                                                          ccccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCaa                                           <  2:3185958/150‑1 (MQ=255)
                                                                                                          ccccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCaa                                           <  1:4900301/150‑1 (MQ=255)
                                                                                                            ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACa                                         <  1:4010880/150‑1 (MQ=255)
                                                                                                                aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaa                                     >  1:6158/1‑150 (MQ=255)
                                                                                                                aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaa                                     >  1:4672640/1‑150 (MQ=255)
                                                                                                                aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaa                                     <  1:1455159/150‑1 (MQ=255)
                                                                                                                  attattTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                   >  2:931055/1‑150 (MQ=255)
                                                                                                                   ttattaGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACa                                  >  1:3803581/1‑150 (MQ=255)
                                                                                                                   ttattTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACa                                  >  1:2089624/1‑150 (MQ=255)
                                                                                                                     attTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGa                                >  1:1706773/1‑150 (MQ=255)
                                                                                                                        tGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaaa                             <  1:3923182/150‑1 (MQ=255)
                                                                                                                           gTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTa                          >  1:4235808/1‑150 (MQ=255)
                                                                                                                           gTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTa                          >  1:3380407/1‑150 (MQ=255)
                                                                                                                            tCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTAc                         >  2:2825742/1‑150 (MQ=255)
                                                                                                                             cAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTAct                        <  1:4450271/150‑2 (MQ=255)
                                                                                                                              aGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACgg                       <  1:3661003/150‑1 (MQ=255)
                                                                                                                               gCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAAAATGAAACAGAAAATTACGGa                      >  1:3295385/1‑150 (MQ=255)
                                                                                                                                 tCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAtt                    <  2:2450877/150‑1 (MQ=255)
                                                                                                                                 tCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAtt                    >  1:3176974/1‑150 (MQ=255)
                                                                                                                                  cacaAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTa                   >  2:3361321/1‑150 (MQ=255)
                                                                                                                                    caAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGACAGGCAACATGAAACAGAAAATTACGGATTAcc                 >  2:436256/1‑150 (MQ=255)
                                                                                                                                      aaTACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTg               >  2:2756499/1‑150 (MQ=255)
                                                                                                                                       aTACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTgg              >  2:1839060/1‑150 (MQ=255)
                                                                                                                                        tACTCTTTTTCCCTGTTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTATCTCTTTTCATGTGACCGCTACCACGGACAGGCAACATGAAACAGAAAATTACGGATTACCTTGa             <  1:1130243/150‑1 (MQ=255)
                                                                                                                                        tACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGa             <  2:2607085/150‑1 (MQ=255)
                                                                                                                                         aCGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGAc            <  2:3728207/150‑1 (MQ=255)
                                                                                                                                         aCGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGAc            <  1:3104508/150‑1 (MQ=255)
                                                                                                                                         aCGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGAc            <  1:2951680/150‑1 (MQ=255)
                                                                                                                                                  cccTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTAc   <  2:675400/150‑1 (MQ=255)
                                                                                                                                                   ccTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACg  >  2:2063533/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACG  >  CP009273/1364218‑1364512

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