Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 654,192 T→A 6.0% pseudogene (512/594 nt) ybeM → pseudogene, C‑N hydrolase superfamily homology;putative enzyme; Not classified; putative amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273654,1920TA 6.0% 1136.8 / 47.6 332pseudogene (512/594 nt)ybeMpseudogene, C‑N hydrolase superfamily homology;putative enzyme; Not classified; putative amidase
Reads supporting (aligned to +/- strand):  ref base T (152/160);  new base A (10/10);  total (162/170)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

Only 100 of 337 total aligned reads displayed.

GCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATTGTGTGCGCGCGTCGAATGACC  >  CP009273/654043‑654341
                                                                                                                                                     |                                                                                                                                                     
 cttgcCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTg                                                                                                                                                      <  1:4574880/150‑1 (MQ=255)
  ttgcCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGa                                                                                                                                                     <  2:1337244/150‑1 (MQ=255)
  ttgcCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGa                                                                                                                                                     <  1:4305545/150‑1 (MQ=255)
    gcCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACg                                                                                                                                                   <  2:2723807/150‑1 (MQ=255)
    gcCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACg                                                                                                                                                   >  1:4707939/1‑150 (MQ=255)
      cGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGcc                                                                                                                                                 <  1:4134533/150‑1 (MQ=255)
      cGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGcc                                                                                                                                                 <  1:1682861/150‑1 (MQ=255)
       gCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGccc                                                                                                                                                >  1:3062290/1‑150 (MQ=255)
       gCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGAGTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGccc                                                                                                                                                >  2:4930319/1‑150 (MQ=255)
          cGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGaa                                                                                                                                             >  1:1845789/1‑150 (MQ=255)
             gcgcTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACgt                                                                                                                                          >  2:3020451/1‑150 (MQ=255)
              cgcTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACgtg                                                                                                                                         >  1:4496627/1‑150 (MQ=255)
               gcTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGTGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACgtgt                                                                                                                                        <  2:1405486/150‑1 (MQ=255)
                 tGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGc                                                                                                                                      <  1:1584862/150‑1 (MQ=255)
                     tACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCa                                                                                                                                  <  1:5019380/150‑1 (MQ=255)
                      accaccTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGAGACGCCCGAACGTGTGCGTCAg                                                                                                                                 >  1:2982221/1‑150 (MQ=255)
                        caccTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGt                                                                                                                               >  2:4838175/1‑150 (MQ=255)
                         accTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTg                                                                                                                              <  1:734671/150‑1 (MQ=255)
                         accTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTg                                                                                                                              >  1:1363635/1‑150 (MQ=255)
                           cTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTgcg                                                                                                                            >  1:1108049/1‑150 (MQ=255)
                            tGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTgcgc                                                                                                                           <  1:2849734/150‑1 (MQ=255)
                             gTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTgcgcg                                                                                                                          <  1:2370839/150‑1 (MQ=255)
                                    ggtggCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACt                                                                                                                   >  2:2236719/1‑150 (MQ=255)
                                    ggtggCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACt                                                                                                                   <  1:396604/150‑1 (MQ=255)
                                       ggcggcGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGcc                                                                                                                <  2:2670971/150‑1 (MQ=255)
                                          ggcggGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGt                                                                                                             >  1:199144/1‑150 (MQ=255)
                                            cggGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCt                                                                                                           >  1:3108380/1‑150 (MQ=255)
                                             gggggAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCtt                                                                                                          <  1:2655116/150‑1 (MQ=255)
                                                  aGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTTCGCGCGCAACTGCCCGTCTTAaaca                                                                                                     <  2:2889341/150‑1 (MQ=255)
                                                  aGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAaaca                                                                                                     >  1:3563164/1‑150 (MQ=255)
                                                  aGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAaaca                                                                                                     <  1:944674/150‑1 (MQ=255)
                                                    tGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAaacaac                                                                                                   >  1:4605631/1‑150 (MQ=255)
                                                        ggAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTc                                                                                               <  2:1755706/150‑1 (MQ=255)
                                                        ggAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTc                                                                                               <  2:4315987/150‑1 (MQ=255)
                                                               aaaTATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTgc                                                                                        <  1:2201640/150‑1 (MQ=255)
                                                               aaaTATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTgc                                                                                        <  1:4042061/150‑1 (MQ=255)
                                                                 atatCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTgcgc                                                                                      >  2:559508/1‑150 (MQ=255)
                                                                 atatCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTgcgc                                                                                      <  2:3255684/150‑1 (MQ=255)
                                                                  tatCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGcgcc                                                                                     <  1:1020993/150‑1 (MQ=255)
                                                                   atCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGcgccg                                                                                    <  2:2324123/150‑1 (MQ=255)
                                                                    tCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGAGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGcgccgc                                                                                   >  2:2104432/1‑150 (MQ=255)
                                                                      ggTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGcgccgccg                                                                                 <  2:1608843/150‑1 (MQ=255)
                                                                        tCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCa                                                                               <  1:1304824/150‑1 (MQ=255)
                                                                        tCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCa                                                                               <  2:3682575/150‑1 (MQ=255)
                                                                        tCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGAGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCa                                                                               <  2:3920236/150‑1 (MQ=255)
                                                                         cAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCaa                                                                              <  1:275692/150‑1 (MQ=255)
                                                                          aaaGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAat                                                                             >  2:3278052/1‑150 (MQ=255)
                                                                          aaaGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAat                                                                             >  2:4888933/1‑150 (MQ=255)
                                                                                ggATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATtattat                                                                       >  2:336468/1‑150 (MQ=255)
                                                                                ggATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATtattat                                                                       <  1:3472886/150‑1 (MQ=255)
                                                                                 gATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAAGTATTatg                                                                      >  2:3084232/1‑150 (MQ=255)
                                                                                   ttATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTatgat                                                                    <  2:2350091/150‑1 (MQ=255)
                                                                                     ataGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGt                                                                  >  1:1951048/1‑150 (MQ=255)
                                                                                        gATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGtttt                                                               <  1:2232119/150‑1 (MQ=255)
                                                                                        gATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGtttt                                                               >  1:4025881/1‑150 (MQ=255)
                                                                                         aTCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGttttt                                                              <  2:2545317/150‑1 (MQ=255)
                                                                                          tCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGtttttt                                                             >  1:4462213/1‑150 (MQ=255)
                                                                                              tttGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGAGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACt                                                         >  2:346071/1‑150 (MQ=255)
                                                                                               ttGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTc                                                        <  2:4306639/150‑1 (MQ=255)
                                                                                                tGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCg                                                       <  1:670896/150‑1 (MQ=255)
                                                                                                   cGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGgcg                                                    <  1:559508/150‑1 (MQ=255)
                                                                                                     tCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCt                                                  >  2:3536030/1‑150 (MQ=255)
                                                                                                     tCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCt                                                  <  1:1255098/150‑1 (MQ=255)
                                                                                                     tCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGAGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCt                                                  <  1:3309522/150‑1 (MQ=255)
                                                                                                      cACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCtt                                                 <  1:1722133/150‑1 (MQ=255)
                                                                                                        ccATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGa                                               >  1:3775891/1‑150 (MQ=255)
                                                                                                         cATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGAt                                              >  2:4099379/1‑150 (MQ=255)
                                                                                                         cATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGAt                                              >  2:3146427/1‑150 (MQ=255)
                                                                                                         cATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGAt                                              >  1:1192197/1‑150 (MQ=255)
                                                                                                         cATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGAt                                              >  2:4951903/1‑150 (MQ=255)
                                                                                                         cATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGAt                                              >  1:3765087/1‑150 (MQ=255)
                                                                                                           gtGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTc                                            <  2:277456/149‑1 (MQ=255)
                                                                                                              cGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCAcc                                         <  2:106426/150‑1 (MQ=255)
                                                                                                               ggCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCt                                        >  1:4952305/1‑150 (MQ=255)
                                                                                                                 cAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTg                                      >  1:3225305/1‑150 (MQ=255)
                                                                                                                   gCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGtt                                    <  2:1486432/150‑1 (MQ=255)
                                                                                                                    cGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTa                                   >  2:2705978/1‑150 (MQ=255)
                                                                                                                     gTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTAc                                  >  2:4383638/1‑150 (MQ=255)
                                                                                                                      tCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACa                                 <  2:1784854/150‑1 (MQ=255)
                                                                                                                         gAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGAt                              >  1:630604/1‑150 (MQ=255)
                                                                                                                         gAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGAt                              >  2:1987175/1‑150 (MQ=255)
                                                                                                                         gAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGAt                              <  1:4981018/150‑1 (MQ=255)
                                                                                                                           aaTGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCACTTATTCGGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGAGTg                            >  2:2952608/1‑150 (MQ=255)
                                                                                                                             tGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCt                          <  2:4199613/150‑1 (MQ=255)
                                                                                                                              gCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTa                         >  1:490218/1‑150 (MQ=255)
                                                                                                                              gCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTa                         <  2:4506193/150‑1 (MQ=255)
                                                                                                                              gCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTAGTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGcga                         >  1:3653907/1‑148 (MQ=255)
                                                                                                                               ccTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTAt                        >  2:4800714/1‑150 (MQ=255)
                                                                                                                                cTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTAtt                       <  1:2596106/150‑1 (MQ=255)
                                                                                                                                cTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTAtt                       >  2:3848086/1‑150 (MQ=255)
                                                                                                                                 tGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATtg                      <  1:2274103/150‑1 (MQ=255)
                                                                                                                                 tGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATtg                      <  1:5076431/150‑1 (MQ=255)
                                                                                                                                    aCTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATtgtgt                   >  1:1938251/1‑150 (MQ=255)
                                                                                                                                     cTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTCTtgtgtg                  >  1:1499684/1‑150 (MQ=255)
                                                                                                                                     cTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTAAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCATGATTCACCTTGTTACAGATTGCTATtgtgtg                  >  2:743262/1‑150 (MQ=255)
                                                                                                                                      tCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATTGTGTgc                 >  1:816871/1‑150 (MQ=255)
                                                                                                                                      tCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATTGTGTgc                 <  2:1157077/150‑1 (MQ=255)
                                                                                                                                                 gAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATTGTGTGCGCGCGTCGAAt      <  1:4109544/150‑1 (MQ=255)
                                                                                                                                                 gAAGAGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATTGTGTGCGCGCGTCGAAt      <  1:1580240/150‑1 (MQ=255)
                                                                                                                                                     tGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATTGTGTGCGCGCGTCGAATGAcc  <  2:2507443/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTCATTATGGCGGAAGTGACGCCCGAACGTGTGCGTCAGGTGCGCGCGCAACTGCCCGTCTTAAACAACCGTCGCTTTGCGCCGCCGCAATTATTATGATGTTTTTTTACTCGGCGCTTGATTCACCTTGTTACAGATTGCTATTGTGTGCGCGCGTCGAATGACC  >  CP009273/654043‑654341

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: