Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,045,280 G→T 19.4% intergenic (‑80/‑9) gfcA ← / → insA O‑antigen capsule production threonine‑rich inner membrane protein/IS1 repressor TnpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,045,2800GT19.4% 730.7 / 122.4 267intergenic (‑80/‑9)gfcA/insAO‑antigen capsule production threonine‑rich inner membrane protein/IS1 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base G (121/94);  new base T (32/19);  total (154/113)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 268 total aligned reads displayed.

GGGGCAGCAAACGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCT  >  CP009273/1045133‑1045429
                                                                                                                                                   |                                                                                                                                                     
ggggCAGCAAACGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGt                                                                                                                                                     <  1:3762317/150‑1 (MQ=255)
 gggCAGCAAACGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGtt                                                                                                                                                    >  2:3419682/1‑150 (MQ=255)
     agcaAACGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTcc                                                                                                                                                <  1:3187047/150‑1 (MQ=255)
     agcaAACGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTcc                                                                                                                                                <  1:1714619/150‑1 (MQ=255)
      gcaAACGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCa                                                                                                                                               >  2:3911973/1‑150 (MQ=255)
       caAACGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAg                                                                                                                                              <  2:113669/150‑1 (MQ=255)
         aaCGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTg                                                                                                                                            >  1:3459553/1‑150 (MQ=255)
          aCGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTgg                                                                                                                                           >  2:3488577/1‑150 (MQ=255)
          aCGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTgg                                                                                                                                           >  2:4644771/1‑150 (MQ=255)
          aCGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTgg                                                                                                                                           >  2:1164843/1‑150 (MQ=255)
           cGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGc                                                                                                                                          >  1:2850906/1‑150 (MQ=255)
            gCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCt                                                                                                                                         <  2:2970623/150‑1 (MQ=255)
             cAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCtt                                                                                                                                        <  2:4360376/150‑1 (MQ=255)
               gCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCt                                                                                                                                      <  1:2415815/150‑1 (MQ=255)
                 tGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGt                                                                                                                                    >  1:2615969/1‑150 (MQ=255)
                  ggCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGtt                                                                                                                                   >  2:4186785/1‑150 (MQ=255)
                    cgtcgtCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTc                                                                                                                                 >  2:1301119/1‑150 (MQ=255)
                     gtcgtcAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCt                                                                                                                                <  2:2877422/150‑1 (MQ=255)
                       cgtcAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTAt                                                                                                                              >  2:2510574/1‑150 (MQ=255)
                                 aaGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTccc                                                                                                                    <  1:3561943/150‑1 (MQ=255)
                                  aGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcct                                                                                                                   <  2:4582946/150‑1 (MQ=255)
                                   ggCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctc                                                                                                                  >  2:4501177/1‑150 (MQ=255)
                                      cATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTg                                                                                                               <  1:4068728/150‑1 (MQ=255)
                                      cATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCGGTCCCTCCTg                                                                                                               >  1:1228431/1‑150 (MQ=255)
                                        tAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACTGTTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGtt                                                                                                             <  2:4089478/150‑1 (MQ=255)
                                        tAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGtt                                                                                                             >  1:3844647/1‑150 (MQ=255)
                                        tAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACCGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGtt                                                                                                             >  1:3016835/1‑150 (MQ=255)
                                           ggATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAg                                                                                                          >  1:4247362/1‑150 (MQ=255)
                                            gATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCAACGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGc                                                                                                         >  1:3775563/1‑150 (MQ=255)
                                             aTTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCt                                                                                                        <  1:2700386/150‑1 (MQ=255)
                                             aTTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCt                                                                                                        <  1:4398148/150‑1 (MQ=255)
                                               tGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTAc                                                                                                      >  2:535601/1‑150 (MQ=255)
                                                 cAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTg                                                                                                    >  1:1871711/1‑150 (MQ=255)
                                                 cAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTg                                                                                                    >  1:2262707/1‑150 (MQ=255)
                                                  aGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGa                                                                                                   <  2:317601/150‑1 (MQ=255)
                                                  aGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACAGa                                                                                                   >  1:4347720/1‑150 (MQ=255)
                                                    aaaGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg                                                                                                 <  2:692423/150‑1 (MQ=255)
                                                     aaGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgg                                                                                                >  1:259037/1‑150 (MQ=255)
                                                     aaGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGAGGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgg                                                                                                >  2:817297/1‑150 (MQ=255)
                                                       gCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgggg                                                                                              <  1:287340/150‑1 (MQ=255)
                                                        cTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggt                                                                                             <  2:3955364/150‑1 (MQ=255)
                                                          tgtgTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggtgg                                                                                           >  1:4434496/1‑150 (MQ=255)
                                                          tgtgTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggtgg                                                                                           <  1:228577/150‑1 (MQ=255)
                                                           gtgtTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggtggt                                                                                          <  1:667875/150‑1 (MQ=255)
                                                           gtgtTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggtggt                                                                                          <  1:1641623/150‑1 (MQ=255)
                                                               ttCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGt                                                                                      >  2:1940810/1‑150 (MQ=255)
                                                                tCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTa                                                                                     <  2:3972329/150‑1 (MQ=255)
                                                                  aTAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAAc                                                                                   <  1:544080/150‑1 (MQ=255)
                                                                  aTAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAAc                                                                                   >  2:4033767/1‑150 (MQ=255)
                                                                     aCTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGc                                                                                >  1:3037795/1‑150 (MQ=255)
                                                                       ttttCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGcta                                                                              <  2:2633858/150‑3 (MQ=255)
                                                                       ttttCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaa                                                                              <  2:3446050/150‑1 (MQ=255)
                                                                         ttCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaaaa                                                                            <  1:2030536/150‑1 (MQ=255)
                                                                          tCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAg                                                                           >  2:2727012/1‑150 (MQ=255)
                                                                           ccTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGc                                                                          >  2:372233/1‑150 (MQ=255)
                                                                           ccTTTATTCATCGCATGGACAATACGGGTGAAGCTGCCAACTTACTGATTTAGTGAAAGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGc                                                                          >  2:1874958/1‑150 (MQ=255)
                                                                              ttATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACt                                                                       <  1:1691058/150‑1 (MQ=255)
                                                                              ttATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACt                                                                       <  1:2656116/150‑1 (MQ=255)
                                                                               tATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTg                                                                      >  1:1186275/1‑150 (MQ=255)
                                                                                 ttCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGcc                                                                    <  1:1567125/150‑1 (MQ=255)
                                                                                 ttCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGcc                                                                    >  2:1478793/1‑150 (MQ=255)
                                                                                    aTCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGa                                                                 >  2:431959/1‑150 (MQ=255)
                                                                                    aTCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTATTGAGGTTCTCCAGTGGCTTATGTTTCTAACAGCAGTCCCTCCTGATCAGCTACATACGGGGTGGTGCGTAACGGCAAAAGCAATGCCGGa                                                                 >  1:3195496/1‑150 (MQ=14)
                                                                                      cGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACa                                                               >  2:4131577/1‑150 (MQ=255)
                                                                                        cATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATc                                                             >  2:2827606/1‑150 (MQ=255)
                                                                                        cATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATc                                                             >  2:2729753/1‑150 (MQ=255)
                                                                                           ggACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAgc                                                          >  2:2252835/1‑150 (MQ=255)
                                                                                              cAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCt                                                       <  1:2122360/150‑1 (MQ=38)
                                                                                                    gggTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCGCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATctct                                                 >  1:1358551/1‑150 (MQ=17)
                                                                                                     ggTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTg                                                <  2:20557/150‑1 (MQ=32)
                                                                                                     ggTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTATAACAACTTTCCCAACTGTTCAGCAACTGAAGGGGTGGTGCGTAACGGCAAAAGAAATGACGGACATCAGAGCTAtaactg                                                >  1:3897589/1‑145 (MQ=11)
                                                                                                       tGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGct                                              >  2:3347527/1‑150 (MQ=32)
                                                                                                       tGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGct                                              <  1:1043537/150‑1 (MQ=32)
                                                                                                        gATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGctc                                             >  2:3160123/1‑150 (MQ=32)
                                                                                                          tgctgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCtctc                                           <  1:3501240/150‑1 (MQ=32)
                                                                                                             tgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACt                                        <  1:2536242/150‑1 (MQ=17)
                                                                                                             tgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACt                                        >  1:1536470/1‑150 (MQ=17)
                                                                                                             tgcCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACt                                        <  1:1009521/150‑1 (MQ=32)
                                                                                                               ccAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGc                                      >  2:1680001/1‑150 (MQ=32)
                                                                                                               ccAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGc                                      >  1:2010872/1‑150 (MQ=32)
                                                                                                                cAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGcc                                     >  1:2743597/1‑150 (MQ=17)
                                                                                                                 aaCTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCg                                    <  2:4595807/150‑1 (MQ=17)
                                                                                                                 aaCTTACTGATTTAGTGTATGATGGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCg                                    <  2:3647213/150‑1 (MQ=17)
                                                                                                                  aCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGt                                   >  2:3734267/1‑150 (MQ=32)
                                                                                                                            ttAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTTCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGc                         >  1:4642963/1‑150 (MQ=17)
                                                                                                                            ttAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGc                         >  1:455234/1‑150 (MQ=32)
                                                                                                                            ttAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGc                         >  2:1605680/1‑150 (MQ=32)
                                                                                                                            ttAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGc                         <  2:2015772/150‑1 (MQ=32)
                                                                                                                             tAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCa                        >  1:1220357/1‑150 (MQ=32)
                                                                                                                                tgtATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACt                     <  1:296722/150‑1 (MQ=32)
                                                                                                                                      atgGTGTTTTTGATGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGtt               <  1:3419682/150‑1 (MQ=17)
                                                                                                                                          tgtTTTTGAATTTCTCCAGACGCATCTGTTTATATCAGCTGTCCATCATGTACAGCTACTGACGGTGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTAACTCTGCTCTCACTGCCGTAAAACACAGCAACTGCAGTTCAca           >  2:1036707/1‑149 (MQ=11)
                                                                                                                                           gtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCACCTGCAGTTCACtt          >  1:2669156/1‑150 (MQ=17)
                                                                                                                                             ttttGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTac        >  2:936251/1‑150 (MQ=32)
                                                                                                                                             ttttGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGGCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTac        <  1:4461309/150‑1 (MQ=17)
                                                                                                                                              tttGATGTTCTCCAGTGTCTTCTGTTTCTATCAGCTGTCCCTCCTTTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAATGCCGGACATCATCGCTATCTCTGCTCTCACTGCCGTTAAACATGGCAACTGCAGTTCACTTaca       <  1:752902/150‑1 (MQ=11)
                                                                                                                                                tGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACAcc     <  2:1405981/150‑1 (MQ=32)
                                                                                                                                                tGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACAcc     <  2:3293926/150‑1 (MQ=32)
                                                                                                                                                tGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACAcc     <  2:2415063/150‑1 (MQ=32)
                                                                                                                                                   ggTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCt  >  1:2919462/1‑150 (MQ=32)
                                                                                                                                                   |                                                                                                                                                     
GGGGCAGCAAACGCAGCTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCT  >  CP009273/1045133‑1045429

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: