Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,200,194 G→A 5.3% pseudogene (335/1134 nt) aaaE → pseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,200,1940GA 5.3% 262.4 / 3.7 75pseudogene (335/1134 nt)aaaEpseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (31/40);  new base A (2/2);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CCAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGT  >  CP009273/1200058‑1200343
                                                                                                                                        |                                                                                                                                                     
ccAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAAc                                                                                                                                          <  1:1670212/150‑1 (MQ=255)
ccAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAAc                                                                                                                                          <  2:1694405/150‑1 (MQ=255)
ccAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAAc                                                                                                                                          <  2:1694402/150‑1 (MQ=255)
   aaCGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCATGCGAAAAATAACGcc                                                                                                                                       >  2:3443592/1‑150 (MQ=255)
   aaCGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCATGCGAAAAATAACGcc                                                                                                                                       >  2:3443587/1‑150 (MQ=255)
    aCGCTGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCAGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGccc                                                                                                                                      <  2:1304632/150‑1 (MQ=255)
    aCGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGccc                                                                                                                                      <  2:1305723/150‑1 (MQ=255)
      gATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGt                                                                                                                                    >  2:127352/1‑150 (MQ=255)
       aTGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGACAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTc                                                                                                                                   <  2:1140378/150‑1 (MQ=255)
       aTGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTc                                                                                                                                   <  2:1140356/150‑1 (MQ=255)
       aTGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTc                                                                                                                                   <  2:1140377/150‑1 (MQ=255)
         gAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTg                                                                                                                                 >  1:2665744/1‑150 (MQ=255)
          aaCTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTgg                                                                                                                                <  2:4133700/150‑1 (MQ=255)
               ttCGGGATCATCTATACCGTTGAAGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAAc                                                                                                                           >  1:4458789/1‑150 (MQ=255)
                  gggATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTc                                                                                                                        >  1:2149701/1‑150 (MQ=255)
                   ggATCTTGTATACCGTTGATGAAGGTGACGACTGGACCGTCCCGCTGGTGCTGGATTAAGCCAATCCATATATTGGCGTGTAGGTTTATCGCGTATTTTTGTTAAGTCAGCTGCAGCGTGCGAAAAATTACGCCCGTCTGGCAAACGTCt                                                                                                                       <  2:2919072/150‑1 (MQ=255)
                    gATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCtt                                                                                                                      >  1:935427/1‑150 (MQ=255)
                    gATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCtt                                                                                                                      >  1:358268/1‑150 (MQ=255)
                          cTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTaaaac                                                                                                                >  1:2185986/1‑150 (MQ=255)
                               ccGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTaaaacaaaac                                                                                                           <  2:1409428/150‑1 (MQ=255)
                               ccGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCAAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCACTCGTCTTTaacacaatac                                                                                                           >  1:3718565/1‑150 (MQ=255)
                                cGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAaaacaaaaca                                                                                                          >  2:3828154/1‑150 (MQ=255)
                                cGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAaaacaaaaca                                                                                                          >  2:3828166/1‑150 (MQ=255)
                                 gTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAacac                                                                                                         >  2:3270957/1‑150 (MQ=255)
                                        aaGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAata                                                                                                  <  2:3706211/150‑1 (MQ=255)
                                        aaGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAata                                                                                                  <  2:3978965/150‑1 (MQ=255)
                                        aaGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAata                                                                                                  <  2:3705293/150‑1 (MQ=255)
                                           gTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCt                                                                                               <  1:3722166/150‑1 (MQ=255)
                                             gacgacTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTgg                                                                                             <  1:4351725/150‑1 (MQ=255)
                                             gacgacTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTgg                                                                                             <  1:4354193/150‑1 (MQ=255)
                                             gacgacTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTgg                                                                                             <  1:975510/150‑1 (MQ=255)
                                               cgacTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGc                                                                                           <  2:4217395/150‑1 (MQ=255)
                                               cgacTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGc                                                                                           <  2:4217100/150‑1 (MQ=255)
                                               cgacTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGc                                                                                           <  2:4217323/150‑1 (MQ=255)
                                                gacTGGACCGACCCGCTGGTGCTGGAAAAAGCCAATCCAAATCTTGGCGTGTCGGTTTATCGCTAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCg                                                                                          <  1:17776/150‑1 (MQ=255)
                                                gacTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCg                                                                                          <  1:17728/150‑1 (MQ=255)
                                                   tGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCg                                                                                       <  1:2800648/150‑1 (MQ=255)
                                                    ggACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCgg                                                                                      >  1:4030165/1‑150 (MQ=255)
                                                    ggACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCgg                                                                                      >  1:4030376/1‑150 (MQ=255)
                                                             cGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCgg                                                                             <  1:3739459/150‑1 (MQ=255)
                                                             cGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCgg                                                                             <  1:3739585/150‑1 (MQ=255)
                                                                  gTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTAt                                                                        <  1:2625206/150‑1 (MQ=255)
                                                                   tGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTAtt                                                                       <  2:100497/150‑1 (MQ=255)
                                                                     cTGGAAAAAGCCAATCCAAATATTGGCGTTTCGGTTTATCGCTAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTc                                                                     <  2:969096/150‑1 (MQ=255)
                                                                          aaaaGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCt                                                                >  2:4126756/1‑150 (MQ=255)
                                                                          aaaaGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCt                                                                >  2:4125878/1‑150 (MQ=255)
                                                                                     cAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGc                                                     <  1:1988365/150‑1 (MQ=255)
                                                                                     cAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGc                                                     <  1:288743/150‑1 (MQ=255)
                                                                                     cAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGc                                                     <  1:4239081/150‑1 (MQ=255)
                                                                                     cAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGc                                                     <  1:4239024/150‑1 (MQ=255)
                                                                                     cAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGc                                                     <  1:2586572/150‑1 (MQ=255)
                                                                                      aaaTATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCa                                                    >  1:3171904/1‑150 (MQ=255)
                                                                                      aaaTATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCa                                                    <  1:2092817/150‑1 (MQ=255)
                                                                                      aaaTATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCa                                                    >  1:3171913/1‑150 (MQ=255)
                                                                                      aaaTATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCa                                                    >  1:3171948/1‑150 (MQ=255)
                                                                                        atatTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCaga                                                  <  1:3270957/150‑1 (MQ=255)
                                                                                                   tCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGAt                                       >  1:1369162/1‑150 (MQ=255)
                                                                                                     ggTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATaa                                     >  1:4175077/1‑150 (MQ=255)
                                                                                                      gTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATaaa                                    >  1:559697/1‑150 (MQ=255)
                                                                                                      gTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATaaa                                    >  1:2845166/1‑150 (MQ=255)
                                                                                                             gcgAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGa                             >  1:3853413/1‑150 (MQ=255)
                                                                                                                aaTTTTTGTTAAGTCAGCAGCAGCATGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGAccc                          <  1:3443592/150‑1 (MQ=255)
                                                                                                                aaTTTTTGTTAAGTCAGCAGCAGCATGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGAccc                          <  1:3443587/150‑1 (MQ=255)
                                                                                                                 aTTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCt                         >  2:1296690/1‑150 (MQ=255)
                                                                                                                 aTTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCt                         >  2:1296528/1‑150 (MQ=255)
                                                                                                                    tttGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGa                      <  1:1184564/150‑1 (MQ=255)
                                                                                                                           aGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTc               >  2:1490925/1‑150 (MQ=255)
                                                                                                                                  agcagcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGc        >  1:994898/1‑150 (MQ=255)
                                                                                                                                  agcagcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGc        >  1:994902/1‑150 (MQ=255)
                                                                                                                                   gcagcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCa       >  2:3582536/1‑150 (MQ=255)
                                                                                                                                   gcagcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACATGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCa       >  2:3582554/1‑150 (MQ=255)
                                                                                                                                    cagcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAg      <  2:1153221/150‑1 (MQ=255)
                                                                                                                                    cagcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAg      <  2:1153179/150‑1 (MQ=255)
                                                                                                                                    cagcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGCCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAg      <  2:3426659/150‑1 (MQ=255)
                                                                                                                                     agcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGc     >  1:1917891/1‑150 (MQ=255)
                                                                                                                                        gTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGt  >  2:2503651/1‑150 (MQ=255)
                                                                                                                                        |                                                                                                                                                     
CCAAACGATGAACTGTTCGGGATCATCTATACCGTTGATGAAGGTGACGACTGGACCGACCCGCAGGTGCTGGAAAAAGCCAATCCAAATATTGGCGTGTCGGTTTATCGCGAATTTTTGTTAAGTCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGT  >  CP009273/1200058‑1200343

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: