Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,010,999 C→A 16.5% intergenic (+84/‑172) pqiC → / → rmf DUF330 family putative lipoprotein/ribosome modulation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,010,9990CA16.5% 538.1 / 64.9 192intergenic (+84/‑172)pqiC/rmfDUF330 family putative lipoprotein/ribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base C (86/74);  new base A (15/17);  total (101/91)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

Only 100 of 197 total aligned reads displayed.

GCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGTTTTCTTTTCCACC  >  CP009273/1010852‑1011144
                                                                                                                                                    |                                                                                                                                                 
 cTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCct                                                                                                                                                  <  1:2275152/148‑1 (MQ=255)
  tGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCctc                                                                                                                                                 <  2:3090419/148‑1 (MQ=255)
   ggCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCctct                                                                                                                                                >  2:3227284/1‑148 (MQ=255)
     ccGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTc                                                                                                                                              >  2:900167/1‑148 (MQ=255)
     ccGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTc                                                                                                                                              <  1:1674467/148‑1 (MQ=255)
        gtgtCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTcaca                                                                                                                                           <  2:2451570/148‑1 (MQ=255)
         tgtCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAt                                                                                                                                          >  1:1384398/1‑148 (MQ=255)
         tgtCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAt                                                                                                                                          <  1:93135/148‑1 (MQ=255)
           tCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAttt                                                                                                                                        <  1:2550532/148‑1 (MQ=255)
            cTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACAtttt                                                                                                                                       >  2:2114682/1‑148 (MQ=255)
              ggAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtttttt                                                                                                                                     <  1:2978201/148‑1 (MQ=255)
               gAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTg                                                                                                                                    >  1:805456/1‑148 (MQ=255)
               gAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTg                                                                                                                                    >  2:62786/1‑148 (MQ=255)
                  tCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTa                                                                                                                                  <  2:3128827/148‑1 (MQ=255)
                  tCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTAc                                                                                                                                 >  1:1535555/1‑148 (MQ=255)
                      gAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAAcc                                                                                                                             <  1:60327/148‑1 (MQ=255)
                        aGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGAAAAATTTGTTCCTCTTAACATTTTTTTTACAACAGa                                                                                                                           >  1:2256053/1‑148 (MQ=255)
                         gCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGAc                                                                                                                          <  2:1171425/148‑1 (MQ=255)
                         gCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGAc                                                                                                                          >  1:175689/1‑148 (MQ=255)
                          ccGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACa                                                                                                                         >  1:259190/1‑148 (MQ=255)
                               tCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTTTTACTCTTCACATTTTTTGTACAACCGACATGccc                                                                                                                    >  1:2298724/1‑148 (MQ=255)
                               tCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGccc                                                                                                                    <  2:1335906/148‑1 (MQ=255)
                                   ttGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCgtgt                                                                                                                >  1:1946475/1‑148 (MQ=255)
                                     gCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAg                                                                                                              <  2:1635476/148‑1 (MQ=255)
                                     gCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAg                                                                                                              <  2:439966/148‑1 (MQ=255)
                                          agagATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTcac                                                                                                         >  1:2206978/1‑148 (MQ=255)
                                           gagaTAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTcaca                                                                                                        >  2:2369852/1‑148 (MQ=255)
                                            agaTAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACaa                                                                                                       >  1:1106469/1‑148 (MQ=255)
                                            agaTAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACaa                                                                                                       >  1:659476/1‑148 (MQ=255)
                                                aaaGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAatat                                                                                                   <  1:2419724/148‑1 (MQ=255)
                                                 aaGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATg                                                                                                  <  1:39493/148‑1 (MQ=255)
                                                    cGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACa                                                                                               <  2:2859506/148‑1 (MQ=255)
                                                    cGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAAAATGACa                                                                                               >  1:2174294/1‑148 (MQ=255)
                                                        tACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAAACGACATGCCCGTGTAGCTCACAAATATGACAGTgg                                                                                           >  1:646914/1‑148 (MQ=255)
                                                         aCCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGc                                                                                          >  1:791277/1‑148 (MQ=255)
                                                           cTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGt                                                                                        <  1:2427846/148‑1 (MQ=255)
                                                             tAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGa                                                                                      <  2:2316417/148‑1 (MQ=255)
                                                              aaTTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGaa                                                                                     >  1:2813721/1‑148 (MQ=255)
                                                                ttATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAAtt                                                                                   <  2:1497237/148‑1 (MQ=255)
                                                                ttATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAAtt                                                                                   <  2:1872078/148‑1 (MQ=255)
                                                                   taAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTg                                                                                <  1:3227284/148‑1 (MQ=255)
                                                                   taAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTg                                                                                <  1:2584843/148‑1 (MQ=255)
                                                                    aaaGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTgc                                                                               >  2:2915956/1‑148 (MQ=255)
                                                                        aTTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCAt                                                                           >  2:2708817/1‑148 (MQ=255)
                                                                        aTTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCAt                                                                           >  2:1055092/1‑148 (MQ=255)
                                                                        aTTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCAt                                                                           >  1:829768/1‑148 (MQ=255)
                                                                           tGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGa                                                                        <  1:73860/148‑1 (MQ=255)
                                                                           tGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGa                                                                        <  1:74585/148‑1 (MQ=255)
                                                                               aaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGc                                                                    <  1:2240600/148‑1 (MQ=255)
                                                                               aaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGc                                                                    <  1:61255/148‑1 (MQ=255)
                                                                                 tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAg                                                                  >  2:1589945/1‑148 (MQ=255)
                                                                                  ataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATTACAGTGGCGTGAATTTTTCGCCTTGACGGCAGt                                                                 >  1:1449413/1‑148 (MQ=255)
                                                                                  ataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGt                                                                 >  2:1211420/1‑148 (MQ=255)
                                                                                   taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGtt                                                                <  1:2647908/148‑1 (MQ=255)
                                                                                   taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGtt                                                                >  2:556917/1‑148 (MQ=255)
                                                                                     aCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTAt                                                              >  2:2737878/1‑148 (MQ=255)
                                                                                      ccGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATg                                                             <  1:3154735/148‑1 (MQ=255)
                                                                                       cGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGa                                                            <  2:2035547/148‑1 (MQ=255)
                                                                                        gTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGAt                                                           <  1:1202268/148‑1 (MQ=255)
                                                                                         tctcCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGAtt                                                          >  2:2681947/1‑148 (MQ=255)
                                                                                         tctcCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGAtt                                                          >  1:3313437/1‑148 (MQ=255)
                                                                                          ctcCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTc                                                         >  1:1816374/1‑148 (MQ=255)
                                                                                            ccGGTATGTTGCCTGAGGCGGTTTTTTTTTATATAACCTGCGGAAAAATTTGTTCCTCTTAACATTTTTTGTAAAACCCACATTCCCGGGTATCTCAAAAATATTACAGTGGCGTGTAATTTGGGCATTGTAGGCAGGTATGATTcgc                                                       >  1:1231685/1‑148 (MQ=255)
                                                                                            ccGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTcgc                                                       >  1:2460417/1‑148 (MQ=255)
                                                                                                tATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGTCGGCAGTTATGATTCGCGGTa                                                   >  2:2626171/1‑148 (MQ=255)
                                                                                                tATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTa                                                   >  2:1978494/1‑148 (MQ=255)
                                                                                                  tGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGAAAAATTTGTTGATCTTCACAATTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGTCGGCAGTTATGATTGGCGGTAtt                                                 >  1:2785521/1‑148 (MQ=255)
                                                                                                  tGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTAtt                                                 >  2:1099681/1‑148 (MQ=255)
                                                                                                  tGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTAtt                                                 >  2:1099171/1‑148 (MQ=255)
                                                                                                   gTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTg                                                <  2:2997013/148‑1 (MQ=255)
                                                                                                    ttGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGc                                               >  2:447926/1‑148 (MQ=255)
                                                                                                        cTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTaa                                           <  1:481861/148‑1 (MQ=255)
                                                                                                             gCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACtgtg                                      <  2:1810402/148‑1 (MQ=255)
                                                                                                             gCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTACTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTTACGGCAGTTTTTATTCGCGGTATTGCTTAAAtgtg                                      >  1:2516432/1‑148 (MQ=255)
                                                                                                              cGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGa                                     <  1:556917/148‑1 (MQ=255)
                                                                                                               ggTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGAt                                    <  1:2976729/148‑1 (MQ=255)
                                                                                                                       ttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACAt                             <  2:1697404/148‑1 (MQ=255)
                                                                                                                           ctctAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAg                         <  2:833725/148‑1 (MQ=255)
                                                                                                                            tctAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGt                        >  1:3236839/1‑148 (MQ=255)
                                                                                                                             ctAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTa                       <  2:2108000/148‑1 (MQ=255)
                                                                                                                             ctAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTa                       <  2:382151/148‑1 (MQ=255)
                                                                                                                                aCGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATc                    >  1:1717080/1‑148 (MQ=255)
                                                                                                                                 cGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCa                   >  1:1631301/1‑148 (MQ=255)
                                                                                                                                 cGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCa                   >  1:1503018/1‑148 (MQ=255)
                                                                                                                                 cGTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCa                   >  2:1560807/1‑148 (MQ=255)
                                                                                                                                  gTGCGGAAATATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCAc                  <  2:1730438/148‑1 (MQ=255)
                                                                                                                                  gTGCGGAAAAATTTGTTCATCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCAc                  <  1:1327728/148‑1 (MQ=255)
                                                                                                                                     cGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGCCATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGt               <  1:900167/148‑1 (MQ=255)
                                                                                                                                     cGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGt               <  1:444836/148‑1 (MQ=255)
                                                                                                                                     cGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGt               <  1:2109985/148‑1 (MQ=255)
                                                                                                                                       gAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGttt             <  1:13757/148‑1 (MQ=255)
                                                                                                                                       gAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGttt             >  2:3262524/1‑148 (MQ=255)
                                                                                                                                       gAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGttt             >  2:3262619/1‑148 (MQ=255)
                                                                                                                                        aaaaaTTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGtttt            <  2:1674206/148‑1 (MQ=255)
                                                                                                                                        aaaaaTTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGtttt            <  2:1384398/148‑1 (MQ=255)
                                                                                                                                         aaaaTTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGttttc           <  2:2934685/148‑1 (MQ=255)
                                                                                                                                          aaaTTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGttttct          >  1:2360748/1‑148 (MQ=255)
                                                                                                                                             tttGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGttttctttt       <  2:791277/148‑1 (MQ=255)
                                                                                                                                               tGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACCTTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGTTTTCTTTTcc     <  2:1946475/148‑1 (MQ=255)
                                                                                                                                                  tCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGTTTTCTTTTccacc  <  2:1310256/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCACAAGAGATAAAGCGTCTACCTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAACTGTGATTGCACATTTAGTAATCACTGTTTTCTTTTCCACC  >  CP009273/1010852‑1011144

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: