Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,491,549 G→A 5.2% intergenic (+143/+371) yjgZ → / ← insG uncharacterized protein/IS4 transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,491,5490GA 5.2% 1201.0 / 33.0 330intergenic (+143/+371)yjgZ/insGuncharacterized protein/IS4 transposase
Reads supporting (aligned to +/- strand):  ref base G (183/129);  new base A (12/5);  total (195/135)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 335 total aligned reads displayed.

GGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTT  >  CP009273/4491401‑4491697
                                                                                                                                                    |                                                                                                                                                    
gggTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTgg                                                                                                                                                     >  1:5063986/1‑150 (MQ=255)
 ggTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTggg                                                                                                                                                    <  1:3983523/150‑1 (MQ=255)
 ggTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTggg                                                                                                                                                    <  1:5149941/150‑1 (MQ=255)
  gTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGATGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTgggg                                                                                                                                                   <  2:4567874/150‑1 (MQ=255)
   tAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGt                                                                                                                                                  >  1:2755413/1‑150 (MQ=255)
    aaGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTc                                                                                                                                                 >  2:3599604/1‑150 (MQ=255)
        gATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTa                                                                                                                                             >  1:707205/1‑150 (MQ=255)
         aTACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTAc                                                                                                                                            >  1:2504227/1‑150 (MQ=255)
         aTACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTAc                                                                                                                                            >  1:2504190/1‑150 (MQ=255)
           aCCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACtt                                                                                                                                          >  1:3286825/1‑150 (MQ=255)
              tGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAg                                                                                                                                       <  1:2426826/150‑1 (MQ=255)
               gCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGa                                                                                                                                      >  2:4985486/1‑150 (MQ=255)
                       gtgtTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCt                                                                                                                              >  2:2032686/1‑150 (MQ=255)
                        tgtTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTc                                                                                                                             <  2:4046234/150‑1 (MQ=255)
                         gtTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGACATAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTcc                                                                                                                            >  2:1663063/1‑150 (MQ=255)
                            tgtgCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGt                                                                                                                         <  2:1576994/150‑1 (MQ=255)
                                cACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCa                                                                                                                     >  1:2614109/1‑150 (MQ=255)
                                    gggCCATAACCTGCAGTGGGGAATGGAACGCCTGACCACAGGAACGGAGCGCGAGCAGTAACCACGCGGGCTGACCATGGGACAGAATATGAGTTAGCCCGGCGTGCTGAGAGGGGTCAGACCGACAGAGCTAGCCTACTGGACCAggcg                                                                                                                 >  2:2769932/1‑147 (MQ=255)
                                     ggCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGActa                                                                                                                <  2:4490338/150‑1 (MQ=255)
                                         aTAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGActagcta                                                                                                            <  1:2948627/150‑1 (MQ=255)
                                         aTAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGActagcta                                                                                                            <  1:310871/150‑1 (MQ=255)
                                          tAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGActagctag                                                                                                           >  1:4572762/1‑150 (MQ=255)
                                           aaCCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACtagctagc                                                                                                          >  2:1320667/1‑150 (MQ=255)
                                           aaCCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACtagctagc                                                                                                          <  2:3747098/150‑1 (MQ=255)
                                             ccTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGccc                                                                                                        <  2:1156888/150‑1 (MQ=255)
                                             ccTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGccc                                                                                                        >  2:3675525/1‑150 (MQ=255)
                                                gCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTAGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTa                                                                                                     <  2:1041020/150‑1 (MQ=255)
                                                 cAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTAc                                                                                                    <  2:2786819/150‑1 (MQ=255)
                                                 cAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTAc                                                                                                    <  2:3919659/150‑1 (MQ=255)
                                                  aGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACg                                                                                                   <  2:2706017/150‑1 (MQ=255)
                                                     tGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCaa                                                                                                >  1:2138756/1‑150 (MQ=255)
                                                         aaTTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCaa                                                                                            <  2:2813449/150‑1 (MQ=255)
                                                           ttGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCaaat                                                                                          >  1:3637790/1‑150 (MQ=255)
                                                            tGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCaaata                                                                                         >  1:1838959/1‑150 (MQ=255)
                                                            tGAACGCCTGACAACAGGAAATGAGCTAGAGCACTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCaaata                                                                                         <  2:4251257/150‑1 (MQ=255)
                                                               aCGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCaaataaat                                                                                      >  2:2098449/1‑150 (MQ=255)
                                                               aCGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCaaataaat                                                                                      <  2:2329680/150‑1 (MQ=255)
                                                               aCGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCaaataaat                                                                                      >  2:490528/1‑150 (MQ=255)
                                                                     gACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCa                                                                                <  1:2874162/150‑1 (MQ=255)
                                                                        aacaGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCACATAAATAAACCATCg                                                                             >  2:2724607/1‑150 (MQ=255)
                                                                         acaGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGTGTCAGAATATGAGTTAGCCCGGCTTTATTCCTGGGTTCAGTACTTCAGATCTAGCCTCCGGTACCAGTCTAGCTAGCCCATACGCAACCCAATAAATCAACCATctt                                                                            >  2:4128299/1‑148 (MQ=255)
                                                                         acaGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGt                                                                            >  2:529391/1‑150 (MQ=255)
                                                                            ggAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAACTAAATCAACCATCGTCag                                                                         <  1:4963714/150‑1 (MQ=255)
                                                                            ggAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCag                                                                         >  1:535664/1‑150 (MQ=255)
                                                                            ggAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCag                                                                         >  1:4093531/1‑150 (MQ=255)
                                                                            ggAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCag                                                                         >  1:3889423/1‑150 (MQ=255)
                                                                            ggAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCag                                                                         >  1:3167453/1‑150 (MQ=255)
                                                                               aaTGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGt                                                                      >  2:464064/1‑150 (MQ=255)
                                                                                aTGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTg                                                                     <  2:1406298/150‑1 (MQ=255)
                                                                                  gAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGCTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCa                                                                   <  2:4762320/150‑1 (MQ=255)
                                                                                   aGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCaa                                                                  >  1:2060597/1‑150 (MQ=255)
                                                                                     cTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCACTc                                                                <  2:949108/150‑1 (MQ=255)
                                                                                       agagCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGt                                                              >  1:1635314/1‑150 (MQ=255)
                                                                                           cAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCat                                                          >  1:927866/1‑150 (MQ=255)
                                                                                              tAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCAtata                                                       >  1:4256995/1‑150 (MQ=255)
                                                                                               aaaCACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATac                                                      >  2:844135/1‑150 (MQ=255)
                                                                                               aaaCACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATac                                                      >  2:3467159/1‑150 (MQ=255)
                                                                                               aaaCACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATac                                                      >  1:564838/1‑150 (MQ=255)
                                                                                               aaaCACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTAGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATac                                                      >  2:1394032/1‑150 (MQ=255)
                                                                                               aaaCACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCAGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATac                                                      >  1:1121625/1‑150 (MQ=255)
                                                                                                 acacGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATacac                                                    <  1:2238761/150‑1 (MQ=255)
                                                                                                  cacGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACt                                                   >  2:1298542/1‑150 (MQ=255)
                                                                                                  cacGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACt                                                   >  2:3833877/1‑150 (MQ=255)
                                                                                                      cgGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCg                                               >  2:4301355/1‑150 (MQ=255)
                                                                                                      cgGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCg                                               >  1:3813163/1‑150 (MQ=255)
                                                                                                          cTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGa                                           >  1:3680147/1‑150 (MQ=255)
                                                                                                           tGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGAt                                          >  1:1776900/1‑150 (MQ=255)
                                                                                                           tGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGAt                                          >  1:3204073/1‑150 (MQ=255)
                                                                                                           tGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCGAGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGAt                                          <  2:2224113/150‑1 (MQ=255)
                                                                                                            gACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATg                                         >  1:1295129/1‑150 (MQ=255)
                                                                                                              cAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATgcg                                       <  1:4729556/150‑1 (MQ=255)
                                                                                                               aaTGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATgcgc                                      >  1:3628018/1‑150 (MQ=255)
                                                                                                               aaTGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATgcgc                                      >  1:4415562/1‑150 (MQ=255)
                                                                                                                aTGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCt                                     >  1:1542600/1‑150 (MQ=255)
                                                                                                                   ggTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTaaa                                  <  2:2867339/150‑1 (MQ=255)
                                                                                                                        gAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATg                             >  2:1475839/1‑150 (MQ=255)
                                                                                                                        gAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCAGCGTCAGAGTGCAATCGTTCAGATACACTGTCGGTGATGCGCTAAAAAATg                             <  2:2489268/150‑1 (MQ=255)
                                                                                                                         aaTATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGc                            >  2:918485/1‑150 (MQ=255)
                                                                                                                         aaTATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGc                            >  2:668245/1‑150 (MQ=255)
                                                                                                                         aaTATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGc                            <  1:977825/150‑1 (MQ=255)
                                                                                                                          atatGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCg                           >  2:706775/1‑150 (MQ=255)
                                                                                                                             tGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGAAATCGTTCATATAAACTGTCGGTGATGCGCTAAAAAATGCGTTa                        >  2:4333375/1‑150 (MQ=255)
                                                                                                                             tGAGTAAGCCCGGCGGTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGAGCAATCGTTCATATACACTGTCGGTGATGAGCTAAAAAATGCGTTa                        >  2:2204665/1‑150 (MQ=255)
                                                                                                                              gAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTaa                       >  1:3233751/1‑150 (MQ=255)
                                                                                                                                   aGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCTTCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTg                  <  2:1100396/150‑1 (MQ=255)
                                                                                                                                   aGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTg                  >  2:4374408/1‑150 (MQ=255)
                                                                                                                                       cGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGtt              <  2:1920173/150‑1 (MQ=255)
                                                                                                                                       cGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGAGGCGCTAAAAAATGCGTTAATGTTGAGtt              >  1:1214368/1‑150 (MQ=255)
                                                                                                                                        ggCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCACCCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTc             <  2:1622286/150‑1 (MQ=255)
                                                                                                                                         gCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTcc            <  1:4122559/150‑1 (MQ=255)
                                                                                                                                          cGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCt           >  2:2939093/1‑150 (MQ=255)
                                                                                                                                          cGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCt           >  2:300503/1‑150 (MQ=255)
                                                                                                                                           gTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTc          <  2:4232820/150‑1 (MQ=255)
                                                                                                                                            ttATTCATAGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCg         >  2:2366924/1‑150 (MQ=255)
                                                                                                                                              aTTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGaa       <  1:788983/150‑1 (MQ=255)
                                                                                                                                              aTTCCTGGGGTCAGTACTTAAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGaa       <  2:5196175/150‑1 (MQ=255)
                                                                                                                                                tCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAAcc     <  1:3493576/150‑1 (MQ=255)
                                                                                                                                                 ccTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAAccc    >  2:1839290/1‑150 (MQ=255)
                                                                                                                                                  cTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCt   >  1:351918/1‑150 (MQ=255)
                                                                                                                                                   tGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCtt  >  1:2385989/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGTAAGTGATACCTGCCGTCAGGTGTTTGTGCACCGGGCCATAACCTGCAGTTGGGAATTGAACGCCTGACAACAGGAAATGAGCTAGAGCAGTAAACACGCGGGCTGACAATGGGTCAGAATATGAGTTAGCCCGGCGTTATTCCTGGGGTCAGTACTTCAGATCTAGCCTCCGGTACCAGACTAGCTAGCCCCTACGCAACCAAATAAATCAACCATCGTCAGAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTT  >  CP009273/4491401‑4491697

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: