Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,932,241 G→T 7.6% D219Y (GAT→TAT)  rbsR → transcriptional repressor of ribose metabolism

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,932,2410GT 7.6% 682.1 / 18.9 196D219Y (GAT→TAT) rbsRtranscriptional repressor of ribose metabolism
Reads supporting (aligned to +/- strand):  ref base G (81/98);  new base T (7/10);  total (88/108)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

ATATCTGATCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACAGGTTCCGCAGGAT  >  CP009273/3932096‑3932387
                                                                                                                                                 |                                                                                                                                                  
atatCTGATCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTTAt                                                                                                                                                  >  1:2667762/1‑148 (MQ=255)
 tatCTGATCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTTAtt                                                                                                                                                 <  1:1553568/148‑1 (MQ=255)
      gATCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGaa                                                                                                                                            <  2:1348370/148‑1 (MQ=255)
      gATCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGaa                                                                                                                                            <  2:1410965/148‑1 (MQ=255)
      gATCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGaa                                                                                                                                            >  2:2816216/1‑148 (MQ=255)
       aTCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAAt                                                                                                                                           >  2:2571851/1‑148 (MQ=255)
        tCGATACAGGTCATACCCGTATCGCCTGTTTTACCGGCCCGCTGGATAAAACTCCGTCGCGCCTTCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTACTGATGGCTATGAAGTCACTGGTGATTTTGAAtt                                                                                                                                          <  1:124039/148‑1 (MQ=255)
        tCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAAtt                                                                                                                                          <  1:310222/148‑1 (MQ=255)
        tCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAAtt                                                                                                                                          <  1:231073/148‑1 (MQ=255)
         cGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTTATTTTGAAttt                                                                                                                                         >  2:1908911/1‑148 (MQ=255)
         cGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAAttt                                                                                                                                         <  1:261451/148‑1 (MQ=255)
           atatAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGTAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTaa                                                                                                                                       <  2:2684308/148‑1 (MQ=255)
           aTAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTaa                                                                                                                                       <  1:762097/148‑1 (MQ=255)
            tAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAAc                                                                                                                                      <  2:2494817/148‑1 (MQ=255)
            tAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAAc                                                                                                                                      <  2:3408494/148‑1 (MQ=255)
            tAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCGATGAAGTCACTGGTGATTTTGAAGTTAAc                                                                                                                                      >  2:2471402/1‑148 (MQ=255)
             aaaGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAAcg                                                                                                                                     >  1:2239709/1‑148 (MQ=255)
               aGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAAcggc                                                                                                                                   <  1:1819179/148‑1 (MQ=255)
                ggTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAAcggcg                                                                                                                                  <  1:1742659/148‑1 (MQ=255)
                  tCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAAGGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTTATTTTGAATTTAACGGCggg                                                                                                                                <  1:3121045/148‑1 (MQ=255)
                   cATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGATGTTTATCGGGCGGCTATGAAACGTGCGGGTCTAAACATTCCTTATGGCTATTATGTCACTGGTGATTTTGAATTTCACGGCGGGt                                                                                                                               <  2:1518288/148‑1 (MQ=255)
                   cATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGt                                                                                                                               <  2:555473/148‑1 (MQ=255)
                   cATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGt                                                                                                                               <  2:1093525/148‑1 (MQ=255)
                    aTACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGtt                                                                                                                              >  1:925349/1‑148 (MQ=255)
                      aCCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATGGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGTTTTTGAATTTAACGGCGGGTTTg                                                                                                                            <  1:403117/148‑1 (MQ=255)
                       cccGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGa                                                                                                                           <  2:1724537/148‑1 (MQ=255)
                         cGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACg                                                                                                                         >  2:189653/1‑148 (MQ=255)
                           tATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAAGGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCt                                                                                                                       <  1:3322121/148‑1 (MQ=255)
                           tATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCt                                                                                                                       <  2:1226823/148‑1 (MQ=255)
                            aTCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTa                                                                                                                      <  2:495292/148‑1 (MQ=255)
                              cGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATg                                                                                                                    >  2:626414/1‑148 (MQ=255)
                               gCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATgc                                                                                                                   >  2:2983586/1‑148 (MQ=255)
                                ccTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATgcg                                                                                                                  >  2:237486/1‑148 (MQ=255)
                                  tGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGcc                                                                                                                >  1:1993029/1‑148 (MQ=255)
                                   gTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCa                                                                                                               >  2:3291347/1‑148 (MQ=255)
                                   gTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCa                                                                                                               >  2:2625335/1‑148 (MQ=255)
                                    tATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCaa                                                                                                              <  2:1051006/148‑1 (MQ=255)
                                    tATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCaa                                                                                                              >  2:1449017/1‑148 (MQ=255)
                                           ggCCCGCTGGATAAAACTCGGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTAt                                                                                                       <  2:2545436/148‑1 (MQ=255)
                                           ggCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTAt                                                                                                       <  1:1458923/148‑1 (MQ=255)
                                            gCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATc                                                                                                      <  2:3121288/148‑1 (MQ=255)
                                               cGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTAAACGGCGGGTTTGACGCTATGCGCCAACTGCTATcaca                                                                                                   <  1:307442/148‑1 (MQ=255)
                                                  tGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATcc                                                                                                >  2:2552335/1‑148 (MQ=255)
                                                    gATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGc                                                                                              <  1:3291347/148‑1 (MQ=255)
                                                    gATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGc                                                                                              >  2:236261/1‑148 (MQ=255)
                                                      tAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTg                                                                                            >  2:152775/1‑148 (MQ=255)
                                                          aCTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTc                                                                                        <  2:3327523/148‑1 (MQ=255)
                                                            tCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCt                                                                                      <  1:1449017/148‑1 (MQ=255)
                                                             ccGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTc                                                                                     >  2:33996/1‑148 (MQ=255)
                                                             ccGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTc                                                                                     <  1:1063067/148‑1 (MQ=255)
                                                             ccGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTc                                                                                     <  1:1063060/148‑1 (MQ=255)
                                                             ccGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTc                                                                                     <  1:3166378/148‑1 (MQ=255)
                                                             ccGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACGGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTc                                                                                     <  1:1397526/148‑1 (MQ=255)
                                                                 cgcgcCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGc                                                                                 <  1:2010741/148‑1 (MQ=255)
                                                                 cgcgcCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGc                                                                                 <  1:3503471/148‑1 (MQ=255)
                                                                     ccTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTc                                                                             <  1:3484479/148‑1 (MQ=255)
                                                                      cTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCt                                                                            <  2:1839987/148‑1 (MQ=255)
                                                                       tGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCtt                                                                           <  1:2527891/148‑1 (MQ=255)
                                                                             tGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCgg                                                                     <  1:3270118/148‑1 (MQ=255)
                                                                               gAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGaa                                                                   <  2:2239709/148‑1 (MQ=255)
                                                                                aaGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGaaa                                                                  <  1:34281/148‑1 (MQ=255)
                                                                                  ggTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATg                                                                <  2:252611/148‑1 (MQ=255)
                                                                                    ttATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGAc                                                              >  1:2761317/1‑148 (MQ=255)
                                                                                     tATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTTATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACg                                                             <  1:808273/148‑1 (MQ=255)
                                                                                     tATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACg                                                             >  1:2283436/1‑148 (MQ=255)
                                                                                            cggcgATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGc                                                      <  2:1100911/148‑1 (MQ=255)
                                                                                               cgATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGt                                                   <  2:1247393/148‑1 (MQ=255)
                                                                                                 aTGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTTATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTg                                                 >  1:2846420/1‑148 (MQ=255)
                                                                                                 aTGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTg                                                 >  1:217767/1‑148 (MQ=255)
                                                                                                  tGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTgg                                                >  1:3250303/1‑148 (MQ=255)
                                                                                                  tGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTgg                                                >  1:2096295/1‑148 (MQ=255)
                                                                                                   gAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTCTGGCTGTTGGc                                               >  1:1875052/1‑148 (MQ=255)
                                                                                                   gAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGc                                               >  2:3380895/1‑148 (MQ=255)
                                                                                                   gAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGc                                               >  1:2060095/1‑148 (MQ=255)
                                                                                                     aaCGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGt                                             >  2:1980673/1‑148 (MQ=255)
                                                                                                     aaCGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGt                                             <  1:164542/148‑1 (MQ=255)
                                                                                                       cGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGttt                                           >  1:807570/1‑148 (MQ=255)
                                                                                                        gTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTa                                          <  2:3436334/148‑1 (MQ=255)
                                                                                                          gCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTAcc                                        <  1:2805276/148‑1 (MQ=255)
                                                                                                                  cAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTta                                <  2:3167917/148‑1 (MQ=255)
                                                                                                                   aaCATTCCTGATGGATATGAAGTCACTGGTTATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTtat                               <  2:1231955/148‑1 (MQ=255)
                                                                                                                     cATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTatat                             <  2:462814/148‑1 (MQ=255)
                                                                                                                     cATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTatat                             >  2:1652695/1‑148 (MQ=255)
                                                                                                                       ttCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCa                           >  1:243141/1‑148 (MQ=255)
                                                                                                                              tGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCagag                    >  1:1191464/1‑148 (MQ=255)
                                                                                                                              tGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCagag                    <  1:2045976/148‑1 (MQ=255)
                                                                                                                               ggCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGt                   >  1:1284044/1‑148 (MQ=255)
                                                                                                                               ggCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGt                   <  1:852225/148‑1 (MQ=255)
                                                                                                                                gCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGtt                  <  2:1718062/148‑1 (MQ=255)
                                                                                                                                gCTATGAAGTCACGGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGtt                  <  1:318528/148‑1 (MQ=255)
                                                                                                                                   aTGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACa               >  2:213719/1‑148 (MQ=255)
                                                                                                                                   aTGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACa               >  2:213949/1‑148 (MQ=255)
                                                                                                                                    tGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACAg              >  1:746371/1‑148 (MQ=255)
                                                                                                                                     gAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACAgg             >  2:424751/1‑148 (MQ=255)
                                                                                                                                      aaGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACAGGt            >  1:2965134/1‑148 (MQ=255)
                                                                                                                                      aaGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACAGGt            <  1:2552288/148‑1 (MQ=255)
                                                                                                                                        gTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACAGGTTc          >  1:546859/1‑148 (MQ=255)
                                                                                                                                              ggTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTTCTATCACATCCGATGCGTCCTCAGGCCGACTTAAACGGAAATGACGCTATGTCTGTTGGCGTTTAACAGGCGTAATATAACGCACAGTTACAGGTTCCGCAgt    >  1:3384731/1‑147 (MQ=255)
                                                                                                                                               gTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACAGGTTCCGCAGGa   <  1:2571534/148‑1 (MQ=255)
                                                                                                                                                tGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACAGGTTCCGCAGGat  <  1:592129/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATATCTGATCGATAAAGGTCATACCCGTATCGCCTGTATTACCGGCCCGCTGGATAAAACTCCGGCGCGCCTGCGGTTGGAAGGTTATCGGGCGGCGATGAAACGTGCGGGTCTCAACATTCCTGATGGCTATGAAGTCACTGGTGATTTTGAATTTAACGGCGGGTTTGACGCTATGCGCCAACTGCTATCACATCCGCTGCGTCCTCAGGCCGTCTTTACCGGAAATGACGCTATGGCTGTTGGCGTTTACCAGGCGTTATATCAGGCAGAGTTACAGGTTCCGCAGGAT  >  CP009273/3932096‑3932387

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: