Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,543,355 T→A 44.7% N184Y (AAC→TAC)  yfeX ← porphyrinogen oxidase, cytoplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,543,3550TA44.7% 79.0 / 271.7 190N184Y (AAC→TAC) yfeXporphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base T (51/54);  new base A (42/43);  total (93/97)
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 196 total aligned reads displayed.

GGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCT  >  CP009273/2543208‑2543501
                                                                                                                                                   |                                                                                                                                                  
ggCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGtt                                                                                                                                                    <  2:3094023/148‑1 (MQ=255)
ggCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGtt                                                                                                                                                    <  2:3093306/148‑1 (MQ=255)
ggCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGtt                                                                                                                                                    <  1:873977/148‑1 (MQ=255)
   gAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTagag                                                                                                                                                 <  2:653677/148‑1 (MQ=255)
     aCAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAgct                                                                                                                                               >  1:1069528/1‑148 (MQ=255)
      cAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAgctg                                                                                                                                              >  2:2961039/1‑148 (MQ=255)
      cAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAgctg                                                                                                                                              >  2:2651682/1‑148 (MQ=255)
        aTCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAgctgct                                                                                                                                            >  2:984190/1‑148 (MQ=255)
        aTCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAgctgct                                                                                                                                            >  2:70495/1‑148 (MQ=255)
          cTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTc                                                                                                                                          <  1:3017161/148‑1 (MQ=255)
               gCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGtt                                                                                                                                     >  1:114812/1‑148 (MQ=255)
                cccTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTtg                                                                                                                                    <  1:2929568/148‑1 (MQ=255)
                cccTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTtg                                                                                                                                    >  2:1235782/1‑148 (MQ=255)
                cccTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTtg                                                                                                                                    >  2:606090/1‑148 (MQ=255)
                  cTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTtgtg                                                                                                                                  >  1:2373795/1‑148 (MQ=255)
                   tttGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTgtgt                                                                                                                                 <  1:1398037/148‑1 (MQ=255)
                      gCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTcc                                                                                                                              <  1:2590945/148‑1 (MQ=255)
                      gCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTcc                                                                                                                              >  1:2035841/1‑148 (MQ=255)
                        cATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCa                                                                                                                            <  2:2844991/148‑1 (MQ=255)
                         aTCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCaa                                                                                                                           >  1:1736324/1‑148 (MQ=255)
                           cttcttTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACg                                                                                                                         >  2:290488/1‑148 (MQ=255)
                           cttcttTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACg                                                                                                                         <  1:1478754/148‑1 (MQ=255)
                             tcttTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCt                                                                                                                       >  2:3291813/1‑148 (MQ=255)
                               tttCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTgg                                                                                                                     >  2:1573141/1‑148 (MQ=255)
                                   aGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACaa                                                                                                                 <  1:1344645/148‑1 (MQ=255)
                                   aGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACaa                                                                                                                 <  1:1344601/148‑1 (MQ=255)
                                   aGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGGTGGCCTCGTTGGTGCGCCCGATCACCATCTCATGTACGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGCCaa                                                                                                                 >  1:1284248/1‑148 (MQ=255)
                                    gATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACaaa                                                                                                                >  2:2623359/1‑148 (MQ=255)
                                    gATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACaaa                                                                                                                >  2:2624031/1‑148 (MQ=255)
                                    gATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACaaa                                                                                                                >  1:3084268/1‑148 (MQ=255)
                                     aTCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAac                                                                                                               <  1:3336405/148‑1 (MQ=255)
                                       cAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAacac                                                                                                             >  2:717412/1‑148 (MQ=255)
                                       cAACGCGGGTGAGGTGAGAGGTTTCCGGAAGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAacac                                                                                                             <  1:1769067/148‑1 (MQ=255)
                                        aaCGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAAcaca                                                                                                            <  1:1838948/148‑1 (MQ=255)
                                        aaCGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAAcaca                                                                                                            <  1:341417/148‑1 (MQ=255)
                                         aCGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACAt                                                                                                           <  2:2110710/148‑1 (MQ=255)
                                          cgcgGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATa                                                                                                          >  1:1049384/1‑148 (MQ=255)
                                          cgcgGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATa                                                                                                          >  2:3257307/1‑148 (MQ=255)
                                          cgcgGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATa                                                                                                          >  2:217890/1‑148 (MQ=255)
                                            cgGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGc                                                                                                        <  2:2423207/148‑1 (MQ=255)
                                               gtgaggtgagAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTgc                                                                                                     <  2:2015103/148‑1 (MQ=255)
                                               gtgaggtgagAGGTTTCCGGACGTTCGTCGCCGTCGATCTCGTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCACCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAAAACATAGCTgc                                                                                                     >  1:1516840/1‑148 (MQ=255)
                                                  aggtgagAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTgccgc                                                                                                  <  2:3415732/148‑1 (MQ=255)
                                                    gtgagAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGccgccc                                                                                                >  1:2951388/1‑148 (MQ=255)
                                                      gagaGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGccgcccgc                                                                                              >  1:2557642/1‑148 (MQ=255)
                                                        gagGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCAt                                                                                            >  2:1075792/1‑148 (MQ=255)
                                                          ggTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATcc                                                                                          <  1:1158528/148‑1 (MQ=255)
                                                          ggTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATcc                                                                                          <  1:841093/148‑1 (MQ=255)
                                                                 ggACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGt                                                                                   <  1:1475799/148‑1 (MQ=255)
                                                                 ggACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGt                                                                                   <  2:1637258/148‑1 (MQ=255)
                                                                  gACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTc                                                                                  <  2:1236031/148‑1 (MQ=255)
                                                                  gACGTTCGTCGCCGTCGATCACATCGTTTGCCACTTTGGTGCGCCCCATAACCAACTCCTGATAGTGAAAGCTCATACGGTAGATCTGCTTCATGTTGTTTTACCAACGCTGGAACAACAAATAGCAGCCGCCCGCAACCACGCCgac                                                                                  >  1:1653833/1‑146 (MQ=255)
                                                                   aCGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCt                                                                                 >  2:795860/1‑148 (MQ=255)
                                                                       tcgtcgCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGa                                                                             <  2:333189/148‑1 (MQ=255)
                                                                        cgtcgCCGTCGATCTCTTCGTTGGCCTCTGTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGAt                                                                            <  2:446516/148‑1 (MQ=255)
                                                                         gtcgCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATa                                                                           <  2:726424/148‑1 (MQ=255)
                                                                          tcgCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATaa                                                                          >  2:2143484/1‑148 (MQ=255)
                                                                          tcgCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATaa                                                                          <  1:2472332/148‑1 (MQ=255)
                                                                             ccGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCg                                                                       <  1:510940/148‑1 (MQ=255)
                                                                                 cGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCAc                                                                   >  1:1159246/1‑148 (MQ=255)
                                                                                 cGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCAc                                                                   <  1:714893/148‑1 (MQ=255)
                                                                                   aTCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACtt                                                                 >  2:1527602/1‑148 (MQ=255)
                                                                                    tctcTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTc                                                                >  1:2617422/1‑148 (MQ=255)
                                                                                    tctcTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTc                                                                >  1:1862666/1‑148 (MQ=255)
                                                                                    tctcTTCGTTGGCCTCTTTGGTGCGCCCGATCACAATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTc                                                                >  2:2598840/1‑148 (MQ=255)
                                                                                     ctctTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTcg                                                               <  2:403253/148‑1 (MQ=255)
                                                                                     ctctTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTcg                                                               >  1:752369/1‑148 (MQ=255)
                                                                                     ctctTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTcg                                                               <  2:1680665/148‑1 (MQ=255)
                                                                                       ctTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTcgcg                                                             >  1:2975848/1‑148 (MQ=255)
                                                                                        ttCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGa                                                            >  1:633333/1‑148 (MQ=255)
                                                                                        ttCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGa                                                            >  1:1403098/1‑148 (MQ=255)
                                                                                             tGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGt                                                       >  1:1644499/1‑148 (MQ=255)
                                                                                             tGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGt                                                       >  1:174692/1‑148 (MQ=255)
                                                                                               gCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGtct                                                     <  1:2785109/148‑1 (MQ=255)
                                                                                                  tctTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCtt                                                  <  1:2463346/148‑1 (MQ=255)
                                                                                                    tttGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCa                                                <  2:2234165/148‑1 (MQ=255)
                                                                                                             cccGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCggg                                       <  1:2239563/148‑1 (MQ=255)
                                                                                                             cccGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCggg                                       <  1:3503109/148‑1 (MQ=255)
                                                                                                                  tCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTc                                  >  2:1298500/1‑148 (MQ=255)
                                                                                                                  tCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTc                                  >  2:3236132/1‑148 (MQ=255)
                                                                                                                  tCACCATCTCCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTc                                  >  1:595760/1‑148 (MQ=255)
                                                                                                                     ccATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGt                               <  2:1456596/148‑1 (MQ=255)
                                                                                                                      cATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTa                              <  1:3220925/148‑1 (MQ=255)
                                                                                                                         ctcCTGATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCg                           >  2:180882/1‑148 (MQ=255)
                                                                                                                           ccTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTc                         >  1:2995048/1‑148 (MQ=255)
                                                                                                                              gATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTcaac                      >  2:2372751/1‑148 (MQ=255)
                                                                                                                              gATCGTGAACGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTcaac                      >  2:1698045/1‑148 (MQ=255)
                                                                                                                               aTCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTcaaca                     >  2:2717740/1‑148 (MQ=255)
                                                                                                                                tCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTcaacaa                    >  1:2444269/1‑148 (MQ=255)
                                                                                                                                tCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTcaacaa                    <  1:2624031/148‑1 (MQ=255)
                                                                                                                                 cGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACaaa                   <  1:17433/148‑1 (MQ=255)
                                                                                                                                  gTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAg                  <  1:808530/148‑1 (MQ=255)
                                                                                                                                   tGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGc                 <  2:373252/148‑1 (MQ=255)
                                                                                                                                   tGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGc                 <  2:855339/148‑1 (MQ=255)
                                                                                                                                      aCGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGc              <  1:217890/148‑1 (MQ=255)
                                                                                                                                      aCGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGc              <  1:3257307/148‑1 (MQ=255)
                                                                                                                                       cGCTCATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCt             <  1:606090/148‑1 (MQ=255)
                                                                                                                                           cATCCGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAgg         <  1:1235222/148‑1 (MQ=255)
                                                                                                                                              ccGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCa      <  1:2842964/148‑1 (MQ=255)
                                                                                                                                              ccGGTAGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCa      <  2:1736324/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTGAGAGGTTTCCGGACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAACGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTGATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCT  >  CP009273/2543208‑2543501

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: