Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,844,245 G→T 12.1% V384V (GTG→GTT sppA → protease IV (signal peptide peptidase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,844,2450GT12.1% 683.4 / 46.0 206V384V (GTG→GTTsppAprotease IV (signal peptide peptidase)
Reads supporting (aligned to +/- strand):  ref base G (111/70);  new base T (9/16);  total (120/86)
Fisher's exact test for biased strand distribution p-value = 1.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 207 total aligned reads displayed.

CGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTATCGG  >  CP009273/1844098‑1844391
                                                                                                                                                   |                                                                                                                                                  
  cAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTAt                                                                                                                                                  <  2:1977364/148‑1 (MQ=255)
    aTTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTc                                                                                                                                                <  1:1653763/148‑1 (MQ=255)
     ttATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTcg                                                                                                                                               <  2:1053484/148‑1 (MQ=255)
       atggatggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTcgcg                                                                                                                                             <  2:2022195/148‑1 (MQ=255)
       atggatggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTcgcg                                                                                                                                             >  1:1503431/1‑148 (MQ=255)
        tggatggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCgcgc                                                                                                                                            <  1:462121/148‑1 (MQ=255)
           atggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGa                                                                                                                                         <  2:2361258/148‑1 (MQ=255)
           atggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGa                                                                                                                                         >  1:1805189/1‑148 (MQ=255)
                gAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTgg                                                                                                                                    <  1:1329633/148‑1 (MQ=255)
                 aGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGgc                                                                                                                                   <  1:2477358/148‑1 (MQ=255)
                    aaaCTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGgcagc                                                                                                                                >  1:3060688/1‑148 (MQ=255)
                     aaCTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGgcagca                                                                                                                               <  2:1154847/148‑1 (MQ=255)
                     aaCTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGgcagca                                                                                                                               >  1:2046283/1‑148 (MQ=255)
                      aCTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGcagcag                                                                                                                              >  1:654468/1‑148 (MQ=255)
                       cTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCagcagc                                                                                                                             >  1:2938735/1‑148 (MQ=255)
                            gggAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCggg                                                                                                                        >  1:3107474/1‑148 (MQ=255)
                               aaTGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAg                                                                                                                     >  1:860026/1‑148 (MQ=255)
                                 tGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCg                                                                                                                   >  1:3045912/1‑148 (MQ=255)
                                  gTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCgg                                                                                                                  <  1:1155339/148‑1 (MQ=255)
                                      gCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTaa                                                                                                              >  1:2186658/1‑148 (MQ=255)
                                         gTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGcc                                                                                                           <  2:2472317/148‑1 (MQ=255)
                                         gTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGcc                                                                                                           <  2:3429736/148‑1 (MQ=255)
                                          tGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCt                                                                                                          >  1:3200164/1‑148 (MQ=255)
                                           gATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtg                                                                                                         >  2:695701/1‑148 (MQ=255)
                                           gATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtg                                                                                                         >  2:695718/1‑148 (MQ=255)
                                           gATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTGAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtg                                                                                                         >  2:1508605/1‑148 (MQ=255)
                                           gATAACACTGCTGCACAAATCCGCGACGCTCGCCTTGAACCGAAAGTGAAAGCGATTGTCCTTCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtg                                                                                                         >  2:2244189/1‑148 (MQ=255)
                                            aTACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtgt                                                                                                        >  2:3277819/1‑148 (MQ=255)
                                            aTACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtgt                                                                                                        >  2:3422403/1‑148 (MQ=255)
                                             tACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtgtg                                                                                                       >  1:1053729/1‑148 (MQ=255)
                                              aCCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTgg                                                                                                      >  2:1866446/1‑148 (MQ=255)
                                              aCCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTgg                                                                                                      <  1:1655860/148‑1 (MQ=255)
                                               ccACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGt                                                                                                     <  1:1036701/148‑1 (MQ=255)
                                                cACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGtt                                                                                                    >  2:265471/1‑148 (MQ=255)
                                                  cTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGt                                                                                                  >  2:3239719/1‑148 (MQ=255)
                                                  cTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGt                                                                                                  >  2:846612/1‑148 (MQ=255)
                                                  cTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGt                                                                                                  >  2:1218635/1‑148 (MQ=255)
                                                    gCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTAt                                                                                                <  2:1093405/148‑1 (MQ=255)
                                                    gCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTAt                                                                                                >  1:2680812/1‑148 (MQ=255)
                                                       gCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGa                                                                                             <  1:286198/148‑1 (MQ=255)
                                                       gCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGa                                                                                             >  2:308778/1‑148 (MQ=255)
                                                        cacaAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGAt                                                                                            >  2:2504548/1‑148 (MQ=255)
                                                        cacaAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGAt                                                                                            <  2:485383/148‑1 (MQ=255)
                                                        cacaAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTGGTATCGAt                                                                                            >  2:3552638/1‑148 (MQ=255)
                                                           aaaTCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATggg                                                                                         >  2:3463956/1‑148 (MQ=255)
                                                           aaaTCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCGGTGGTTGTATCGATggg                                                                                         >  2:2733685/1‑148 (MQ=255)
                                                             aTCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGACTGTGGTTGTATCGATGggcg                                                                                       <  1:886047/148‑1 (MQ=255)
                                                              tCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcgg                                                                                      >  2:1102807/1‑148 (MQ=255)
                                                              tCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcgg                                                                                      >  2:2329206/1‑148 (MQ=255)
                                                               ccGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggc                                                                                     <  1:3239719/148‑1 (MQ=255)
                                                                   gACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatgg                                                                                 >  1:688260/1‑148 (MQ=255)
                                                                     cGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggca                                                                                    >  1:2972863/1‑143 (MQ=255)
                                                                        tCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcggca                                                                            >  1:680094/1‑148 (MQ=255)
                                                                        tCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcggca                                                                            >  1:2707374/1‑148 (MQ=255)
                                                                         cGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcggcat                                                                           <  1:2377374/148‑1 (MQ=255)
                                                                          gCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATc                                                                          <  2:1493117/148‑1 (MQ=255)
                                                                              tGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggt                                                                      >  2:598289/1‑148 (MQ=255)
                                                                               gACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggtg                                                                     >  1:1145382/1‑148 (MQ=255)
                                                                                 cccGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGTAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTggtggt                                                                   <  1:1286833/148‑1 (MQ=255)
                                                                                  ccGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGtt                                                                  >  1:2536464/1‑148 (MQ=255)
                                                                                  ccGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGtt                                                                  >  2:2718801/1‑148 (MQ=255)
                                                                                   cGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTGACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTa                                                                 <  2:968625/148‑1 (MQ=255)
                                                                                            aaGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTc                                                        <  1:3378012/148‑1 (MQ=255)
                                                                                              gCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCa                                                      <  2:3067098/148‑1 (MQ=255)
                                                                                               cGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAACTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCAGTGGTTGTATCGATCGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCAc                                                     >  1:109901/1‑148 (MQ=255)
                                                                                                 aTTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGc                                                   >  2:1762165/1‑148 (MQ=255)
                                                                                                  ttGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGcc                                                  >  1:561444/1‑148 (MQ=255)
                                                                                                  ttGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGcc                                                  >  1:561573/1‑148 (MQ=255)
                                                                                                  ttGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGcc                                                  >  2:1643137/1‑148 (MQ=255)
                                                                                                  ttGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAAAATGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGcc                                                  >  2:2097038/1‑148 (MQ=255)
                                                                                                   tGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCa                                                 >  2:876350/1‑148 (MQ=255)
                                                                                                   tGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCa                                                 >  2:2369599/1‑148 (MQ=255)
                                                                                                     tCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc                                               <  1:2728783/148‑1 (MQ=255)
                                                                                                       cTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTa                                             <  2:3061276/148‑1 (MQ=255)
                                                                                                       cTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTa                                             >  2:3263666/1‑148 (MQ=255)
                                                                                                        tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                            >  2:52028/1‑148 (MQ=255)
                                                                                                        tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                            <  1:376750/148‑1 (MQ=255)
                                                                                                        tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                            <  1:376445/148‑1 (MQ=255)
                                                                                                        tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                            <  1:670488/148‑1 (MQ=255)
                                                                                                        tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                            >  2:3043855/1‑148 (MQ=255)
                                                                                                           gtgtTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTa                                         <  2:1010740/148‑1 (MQ=255)
                                                                                                            tgtTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTAc                                        >  2:1680040/1‑148 (MQ=255)
                                                                                                                  tAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgtg                                  <  2:3379427/148‑1 (MQ=255)
                                                                                                                  tAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgtg                                  >  2:2937842/1‑148 (MQ=255)
                                                                                                                       cAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTaa                             >  2:1318171/1‑148 (MQ=255)
                                                                                                                        aGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAc                            >  1:3532690/1‑148 (MQ=255)
                                                                                                                         ggcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAcc                           <  1:2284248/148‑1 (MQ=255)
                                                                                                                         ggcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAcc                           >  2:2799693/1‑148 (MQ=255)
                                                                                                                              cAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGc                      >  2:653849/1‑148 (MQ=255)
                                                                                                                                  gTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCAccc                  <  1:3425881/148‑1 (MQ=255)
                                                                                                                                    tACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTg                <  1:199298/148‑1 (MQ=255)
                                                                                                                                     aCCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGa               <  1:52028/148‑1 (MQ=255)
                                                                                                                                      ccGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGAc              >  2:1480046/1‑148 (MQ=255)
                                                                                                                                         cgTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCgg           <  1:141165/148‑1 (MQ=255)
                                                                                                                                            cTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTc        >  2:3225913/1‑148 (MQ=255)
                                                                                                                                            cTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTc        >  2:2893830/1‑148 (MQ=255)
                                                                                                                                            cTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTc        >  2:2893075/1‑148 (MQ=255)
                                                                                                                                               aaGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTAt     <  2:1907613/148‑1 (MQ=255)
                                                                                                                                               aaGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTAt     <  2:1908535/148‑1 (MQ=255)
                                                                                                                                               aaGTTATTCGCGCTGAACGGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTAt     <  2:2983652/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTATCGG  >  CP009273/1844098‑1844391

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: