Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,575,640 C→A 7.4% pseudogene (993/3567 nt) hsdR ← pseudogene, premature stop codon, host restriction; endonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,575,6400CA 7.4% 1061.8 / 37.4 279pseudogene (993/3567 nt)hsdRpseudogene, premature stop codon, host restriction; endonuclease R
Reads supporting (aligned to +/- strand):  ref base C (129/129);  new base A (10/11);  total (139/140)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

CAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGT  >  CP009273/4575491‑4575789
                                                                                                                                                     |                                                                                                                                                     
cAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCTGTGAGTAGTGATCAAGCAAGGTTTCCCGCAGGAAGCAATTATAGAAACATTTACTGTAGCGATACGAATCATTGAGCCTGGCGGGAACGACGATATTTAGACGTTTc                                                                                                                                                       >  1:492364/1‑150 (MQ=255)
     ggTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCctc                                                                                                                                                  >  2:367248/1‑150 (MQ=255)
      gTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCctct                                                                                                                                                 <  1:4056709/150‑1 (MQ=255)
      gTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCctct                                                                                                                                                 <  2:1286242/150‑1 (MQ=255)
        ttCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTac                                                                                                                                               >  1:1220053/1‑150 (MQ=255)
          cATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTAGCCTCTacca                                                                                                                                             >  1:1972979/1‑150 (MQ=255)
           aTACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTaccac                                                                                                                                            <  2:4586385/150‑1 (MQ=255)
           aTACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTaccac                                                                                                                                            >  1:1888496/1‑150 (MQ=255)
             aCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACcg                                                                                                                                          >  1:2545275/1‑150 (MQ=255)
                   gCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATgg                                                                                                                                    <  2:2522025/150‑1 (MQ=255)
                    cACTGCTTCATGCTCTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGTGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATggg                                                                                                                                   <  1:56449/150‑1 (MQ=255)
                    cACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATggg                                                                                                                                   <  1:523383/150‑1 (MQ=255)
                     aCTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGt                                                                                                                                  >  1:74014/1‑150 (MQ=255)
                      cTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGACACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGtt                                                                                                                                 <  2:2377373/150‑1 (MQ=255)
                       tGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGttt                                                                                                                                <  1:100528/150‑1 (MQ=255)
                             aTGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGcc                                                                                                                          >  1:2103708/1‑150 (MQ=255)
                                cACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGAc                                                                                                                       >  2:3860937/1‑150 (MQ=255)
                                 aCTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACa                                                                                                                      <  2:1137505/150‑1 (MQ=255)
                                 aCTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACa                                                                                                                      >  1:526670/1‑150 (MQ=255)
                                   ttCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACaaa                                                                                                                    >  1:1695096/1‑150 (MQ=255)
                                     cATGCGGTGAGTAGTGCTCAAGCAATGTTTCCAGAAGGATGCCATTTTCGAAACATTTACTGTTGAGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTAGCCTCTACCACCGCGATGGGTTTGAGGCCGTCAATCa                                                                                                                  <  1:3245972/150‑1 (MQ=255)
                                     cATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACa                                                                                                                  <  1:3245994/150‑1 (MQ=255)
                                      aTCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAg                                                                                                                 <  1:4171339/150‑1 (MQ=255)
                                         cGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGcac                                                                                                              >  2:3584154/1‑150 (MQ=255)
                                          ggTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGcaca                                                                                                             <  1:1064023/150‑1 (MQ=255)
                                           gTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACAt                                                                                                            >  1:2161787/1‑150 (MQ=255)
                                           gTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACAt                                                                                                            <  2:4140195/150‑1 (MQ=255)
                                           gTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTAGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACAt                                                                                                            <  2:4307457/150‑1 (MQ=255)
                                            tGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATa                                                                                                           <  2:1540400/150‑1 (MQ=255)
                                            tGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATa                                                                                                           <  1:1735696/150‑1 (MQ=255)
                                              agtagtGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAAt                                                                                                         <  1:825333/150‑1 (MQ=255)
                                              agtagtGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAAt                                                                                                         <  1:1961580/150‑1 (MQ=255)
                                                      tCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAg                                                                                                 <  2:2699393/150‑1 (MQ=255)
                                                      tCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAg                                                                                                 <  2:632680/150‑1 (MQ=255)
                                                      tCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAg                                                                                                 <  2:4025643/150‑1 (MQ=255)
                                                             agGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGAt                                                                                          >  2:1080807/1‑150 (MQ=255)
                                                             agGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGAt                                                                                          <  1:33815/150‑1 (MQ=255)
                                                               gTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTa                                                                                        >  1:3825950/1‑150 (MQ=255)
                                                               gTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTa                                                                                        >  2:219550/1‑150 (MQ=255)
                                                                tttCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTAc                                                                                       >  1:2037264/1‑150 (MQ=255)
                                                                  tCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTAccc                                                                                     >  1:3783753/1‑150 (MQ=255)
                                                                  tCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTAccc                                                                                     >  2:2739994/1‑150 (MQ=255)
                                                                    ccGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGt                                                                                   >  1:3089393/1‑150 (MQ=255)
                                                                      gCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGttt                                                                                 <  1:971328/150‑1 (MQ=255)
                                                                       cAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCAACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTc                                                                                >  1:3047063/1‑150 (MQ=255)
                                                                           aaGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCt                                                                            <  1:2967008/150‑1 (MQ=255)
                                                                            aGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCtt                                                                           >  2:1678240/1‑150 (MQ=255)
                                                                             gCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCttt                                                                          >  1:3847917/1‑150 (MQ=255)
                                                                              ccATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCtttt                                                                         <  1:2400294/150‑1 (MQ=255)
                                                                                aTTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTcc                                                                       >  2:381302/1‑150 (MQ=255)
                                                                                aTTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTAGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTcc                                                                       <  1:3892225/150‑1 (MQ=255)
                                                                                 ttATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCg                                                                      <  2:4624719/150‑1 (MQ=255)
                                                                                  tATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCgg                                                                     >  1:864048/1‑150 (MQ=255)
                                                                                   aTCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGt                                                                    >  1:4514136/1‑150 (MQ=255)
                                                                                   aTCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGt                                                                    >  2:4416324/1‑150 (MQ=255)
                                                                                   aTCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGt                                                                    <  2:3908469/150‑1 (MQ=255)
                                                                                    tCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTc                                                                   <  1:234927/150‑1 (MQ=255)
                                                                                     cGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCg                                                                  <  1:4583036/150‑1 (MQ=255)
                                                                                       aaaCATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGc                                                                >  1:2366441/1‑150 (MQ=255)
                                                                                         aCATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCa                                                              >  2:3072781/1‑150 (MQ=255)
                                                                                           aTTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCAtt                                                            <  1:28850/150‑1 (MQ=255)
                                                                                              tACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCgg                                                         <  2:3307469/150‑1 (MQ=255)
                                                                                                  ggaGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAAt                                                     <  1:1859261/148‑1 (MQ=255)
                                                                                                  gTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAAt                                                     >  1:2936162/1‑150 (MQ=255)
                                                                                                    aGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTAGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCg                                                                                    >  2:177121/1‑117 (MQ=255)
                                                                                                      cGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCt                                                 >  2:1266096/1‑150 (MQ=255)
                                                                                                       gATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTAGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCtg                                                >  1:469526/1‑149 (MQ=255)
                                                                                                           cGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTAtt                                            >  2:1502028/1‑150 (MQ=255)
                                                                                                           cGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTAtt                                            >  2:4433033/1‑150 (MQ=255)
                                                                                                            gACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTg                                           >  1:1233866/1‑150 (MQ=255)
                                                                                                            gACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTAGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTg                                           >  1:2276406/1‑150 (MQ=255)
                                                                                                              cTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGAc                                         <  1:2762772/150‑1 (MQ=255)
                                                                                                                cATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGc                                       >  2:559750/1‑150 (MQ=255)
                                                                                                                 aTTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGcc                                      >  2:3467030/1‑150 (MQ=255)
                                                                                                                   tGAGCCTGGCGGGAACGTCGATATTGTTCCGTTTCGCCTCTACAACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTCCCCGTTTCATCTTTTCAGGTCTGCCAATCGTCACTGGCTTTATTGACGCCgg                                    >  2:2478413/1‑150 (MQ=255)
                                                                                                                   tGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCgg                                    >  1:3702218/1‑150 (MQ=255)
                                                                                                                       ccTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATACTCCGCAACGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTc                                <  2:260444/150‑1 (MQ=255)
                                                                                                                       ccTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTc                                <  2:3418519/150‑1 (MQ=255)
                                                                                                                        cTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTcc                               >  1:266653/1‑150 (MQ=255)
                                                                                                                        cTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTcc                               <  1:653625/150‑1 (MQ=255)
                                                                                                                         tGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCg                              >  1:2815346/1‑150 (MQ=255)
                                                                                                                         tGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCg                              >  1:641147/1‑150 (MQ=255)
                                                                                                                          ggCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCgg                             <  1:1651685/150‑1 (MQ=255)
                                                                                                                            cGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGAc                           >  2:2923529/1‑150 (MQ=255)
                                                                                                                                 aCGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTgcg                      >  1:3655820/1‑150 (MQ=255)
                                                                                                                                  cGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTgcgc                     >  1:180362/1‑150 (MQ=255)
                                                                                                                                  cGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTgcgc                     >  2:824506/1‑150 (MQ=255)
                                                                                                                                    tCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCt                   >  2:190514/1‑150 (MQ=255)
                                                                                                                                      gATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACCTAATCCGCAAAGCCCTGATAACACGTTTCATCTTTTACGGTCTGCCATTCGGCAATGGCTGTATTGACGCCGGGTTCCGGACGTgcgcatgt                 >  2:2036774/1‑146 (MQ=255)
                                                                                                                                        tatTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTgg               >  1:4492014/1‑150 (MQ=255)
                                                                                                                                          ttGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGgag             <  2:1956055/150‑1 (MQ=255)
                                                                                                                                           tGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGgaga            >  1:4423023/1‑150 (MQ=255)
                                                                                                                                           tGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGgaga            >  2:4344839/1‑150 (MQ=255)
                                                                                                                                             ttACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAg          <  1:1427662/150‑1 (MQ=255)
                                                                                                                                             ttACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAg          >  1:3399849/1‑150 (MQ=255)
                                                                                                                                              tACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAgc         <  1:236738/150‑1 (MQ=255)
                                                                                                                                                cGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAgcgc       <  1:1286195/150‑1 (MQ=255)
                                                                                                                                                 gTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCa      >  1:3872451/1‑150 (MQ=255)
                                                                                                                                                 gTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCa      <  1:652784/150‑1 (MQ=255)
                                                                                                                                                     cGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGt  >  2:2840623/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGT  >  CP009273/4575491‑4575789

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: