Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,844,245 G→T 8.3% 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,2450GT 8.3% 915.9 / 34.0 252V384V (GTG→GTTsppAprotease IV (signal peptide peptidase)
Reads supporting (aligned to +/- strand):  ref base G (118/113);  new base T (5/16);  total (123/129)
Fisher's exact test for biased strand distribution p-value = 2.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 253 total aligned reads displayed.

CAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTATCGGT  >  CP009273/1844100‑1844392
                                                                                                                                                 |                                                                                                                                                   
cAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTAt                                                                                                                                                   <  1:2479090/148‑1 (MQ=255)
cAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGAt                                                                                                                                                   >  2:207396/1‑148 (MQ=255)
    tATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTcgc                                                                                                                                               >  2:2928470/1‑148 (MQ=255)
     atggatggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTcgcg                                                                                                                                              >  2:1375241/1‑148 (MQ=255)
           ggCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAAc                                                                                                                                        >  1:2701694/1‑148 (MQ=255)
                ggAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGgca                                                                                                                                   >  2:1426625/1‑148 (MQ=255)
                 gAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGgcag                                                                                                                                  >  2:3005778/1‑148 (MQ=255)
                 gAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGgcag                                                                                                                                  >  2:3005241/1‑148 (MQ=255)
                     cTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCagcagc                                                                                                                              <  1:185648/148‑1 (MQ=255)
                        aGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCg                                                                                                                           <  1:1242495/148‑1 (MQ=255)
                             aaTGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAg                                                                                                                      <  2:526551/148‑1 (MQ=255)
                             aaTGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAg                                                                                                                      <  2:2059262/148‑1 (MQ=255)
                              aTGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGc                                                                                                                     <  1:1986025/148‑1 (MQ=255)
                                gTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCgg                                                                                                                   >  1:213521/1‑148 (MQ=255)
                                  tGGCGGTGATACCACTGCGGCACAAAGCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGt                                                                                                                 <  2:1838440/148‑1 (MQ=255)
                                      ggTGATACCACTGCGGCACACCTCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGc                                                                                                             <  2:1145591/148‑1 (MQ=255)
                                        tGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCt                                                                                                           >  2:2022730/1‑148 (MQ=255)
                                         gATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtg                                                                                                          >  2:2719236/1‑148 (MQ=255)
                                          aTACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCtgt                                                                                                         >  2:2699797/1‑148 (MQ=255)
                                            aCCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTgg                                                                                                       <  2:2575388/148‑1 (MQ=255)
                                             ccACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGt                                                                                                      >  2:2637359/1‑148 (MQ=255)
                                                 tGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTa                                                                                                  <  1:926539/148‑1 (MQ=255)
                                                    ggCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCg                                                                                               <  1:950841/148‑1 (MQ=255)
                                                    ggCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCg                                                                                               >  1:3154436/1‑148 (MQ=255)
                                                     gCACAACTCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGa                                                                                              <  1:3004527/148‑1 (MQ=255)
                                                     gCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGa                                                                                              <  1:395013/148‑1 (MQ=255)
                                                      cacaAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGAt                                                                                             >  1:80125/1‑148 (MQ=255)
                                                           aTCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATCGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcg                                                                                        >  2:3450523/1‑148 (MQ=255)
                                                           aTCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcg                                                                                        >  2:1980674/1‑148 (MQ=255)
                                                            tCCGCGACGCTCGCCTTGACCCGCAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGAGCGTAAGCCTGTGGGTGTATCGATGggcgg                                                                                       >  1:2990668/1‑148 (MQ=255)
                                                            tCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGACAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcgg                                                                                       >  1:1881888/1‑148 (MQ=255)
                                                             ccGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggc                                                                                      >  2:2471087/1‑148 (MQ=255)
                                                             ccGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggc                                                                                      <  1:2425219/148‑1 (MQ=255)
                                                                     cTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcggc                                                                              >  2:2524486/1‑148 (MQ=255)
                                                                     cTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcggc                                                                              <  1:3255151/148‑1 (MQ=255)
                                                                     cTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcggc                                                                              >  2:308325/1‑148 (MQ=255)
                                                                       cGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcggcat                                                                            <  1:1036550/148‑1 (MQ=255)
                                                                         ccTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCt                                                                          >  1:1451738/1‑148 (MQ=255)
                                                                            tGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggt                                                                       >  2:1868498/1‑148 (MQ=255)
                                                                            tGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggt                                                                       <  2:553343/148‑1 (MQ=255)
                                                                            tGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggt                                                                       >  2:1700449/1‑148 (MQ=255)
                                                                            tGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggt                                                                       <  1:3319638/148‑1 (MQ=255)
                                                                            tGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggt                                                                       <  2:2599142/148‑1 (MQ=255)
                                                                             gACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggtg                                                                      >  1:2770831/1‑148 (MQ=255)
                                                                              aCCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggtgg                                                                     >  2:415510/1‑148 (MQ=255)
                                                                              aCCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggtgg                                                                     <  2:1174057/148‑1 (MQ=255)
                                                                              aCCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCtggtgg                                                                     >  2:2413510/1‑148 (MQ=255)
                                                                                   aaaGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACt                                                                >  2:448810/1‑148 (MQ=255)
                                                                                    aaGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTg                                                               <  2:107169/148‑1 (MQ=255)
                                                                                     aGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTgg                                                              <  2:210810/148‑1 (MQ=255)
                                                                                      gTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGa                                                             <  2:1492852/148‑1 (MQ=255)
                                                                                      gTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGa                                                             <  1:369281/148‑1 (MQ=255)
                                                                                       tGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGAt                                                            <  1:71499/148‑1 (MQ=255)
                                                                                        gAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGAtt                                                           >  1:1422451/1‑148 (MQ=255)
                                                                                          aaGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTc                                                         >  1:1104714/1‑148 (MQ=255)
                                                                                           aGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTcc                                                        >  1:231867/1‑148 (MQ=255)
                                                                                           aGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTcc                                                        >  1:2357240/1‑148 (MQ=255)
                                                                                            gCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCa                                                       >  1:2314224/1‑148 (MQ=255)
                                                                                               atggTCCTGCGTGTTAATAGCCCAGGCGTCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGc                                                    <  2:2129742/145‑1 (MQ=255)
                                                                                               aTTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGc                                                    <  2:498794/148‑1 (MQ=255)
                                                                                               aTTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGc                                                    <  2:2749527/148‑1 (MQ=255)
                                                                                               aTTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGc                                                    >  2:232249/1‑148 (MQ=255)
                                                                                                  gTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAg                                                 <  1:2191887/148‑1 (MQ=255)
                                                                                                   tCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc                                                <  2:1701626/148‑1 (MQ=255)
                                                                                                     cTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGCTTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTa                                              <  2:871192/148‑1 (MQ=255)
                                                                                                      tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                             <  2:3022796/148‑1 (MQ=255)
                                                                                                      tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                             <  2:1253101/148‑1 (MQ=255)
                                                                                                      tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                             >  1:2251512/1‑148 (MQ=255)
                                                                                                        cGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAAtt                                           <  2:2517772/148‑1 (MQ=255)
                                                                                                           gtTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACa                                        >  1:1546586/1‑148 (MQ=255)
                                                                                                           gtTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACa                                        >  1:1684377/1‑148 (MQ=255)
                                                                                                           gtTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACa                                        <  2:3300628/148‑1 (MQ=255)
                                                                                                           gtTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACa                                        <  2:726396/148‑1 (MQ=255)
                                                                                                            ttAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAt                                       >  1:1392384/1‑148 (MQ=255)
                                                                                                            ttAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAt                                       <  1:2330839/148‑1 (MQ=255)
                                                                                                             tAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAtt                                      <  1:479919/148‑1 (MQ=255)
                                                                                                               aTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTTATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgt                                    <  1:2928470/148‑1 (MQ=255)
                                                                                                               aTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTTGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgt                                    >  1:1312465/1‑148 (MQ=255)
                                                                                                                tAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgtg                                   >  2:2810455/1‑148 (MQ=255)
                                                                                                                tAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgtg                                   >  2:2964578/1‑148 (MQ=255)
                                                                                                                tAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgtg                                   >  2:3041054/1‑148 (MQ=255)
                                                                                                                 aGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTgg                                  <  1:2716587/148‑1 (MQ=255)
                                                                                                                   cccAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCt                                <  2:2504594/148‑1 (MQ=255)
                                                                                                                   cccAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCt                                >  2:2315409/1‑148 (MQ=255)
                                                                                                                   cccAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCt                                >  2:555871/1‑148 (MQ=255)
                                                                                                                    ccAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTa                               <  2:381190/148‑1 (MQ=255)
                                                                                                                        gcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAccc                           >  2:346429/1‑148 (MQ=255)
                                                                                                                        gcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAccc                           >  2:2138403/1‑148 (MQ=255)
                                                                                                                           gcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAg                        >  1:3131133/1‑148 (MQ=255)
                                                                                                                           gcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAg                        <  2:2173399/148‑1 (MQ=255)
                                                                                                                            cAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGc                       >  2:3408296/1‑148 (MQ=255)
                                                                                                                               cGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCAcc                    <  2:3412976/148‑1 (MQ=255)
                                                                                                                               cGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCAcc                    <  2:3491798/148‑1 (MQ=255)
                                                                                                                                 ttACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCt                  >  1:1498521/1‑148 (MQ=255)
                                                                                                                                  tACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTg                 >  2:919763/1‑148 (MQ=255)
                                                                                                                                   aCCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGa                <  2:1374540/148‑1 (MQ=255)
                                                                                                                                     cgcgTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAg                        <  1:2616768/138‑1 (MQ=255)
                                                                                                                                          cTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTc         >  1:3554312/1‑148 (MQ=255)
                                                                                                                                            gAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTATTTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTa       <  2:1446619/148‑1 (MQ=255)
                                                                                                                                                 gATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTATCGGt  <  1:544451/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CAATTATGGATGGCGAGGAAACTCAGGGGAATGTTGGCGGTGATACCACTGCGGCACAAATCCGCGACGCTCGCCTTGACCCGAAAGTGAAAGCGATTGTCCTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCGGTTCTATCGGT  >  CP009273/1844100‑1844392

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: