Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,460,517 G→A 6.8% A316A (GCC→GCT)  algJ ← alginate O‑acetylase AlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,460,5170GA 6.8% 174.6 / 2.4 59A316A (GCC→GCT) algJalginate O‑acetylase AlgJ
Reads supporting (aligned to +/- strand):  ref base G (30/25);  new base A (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TCCACTGCGGATCGAAGCTGCTGAGGTCGTTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGA  >  NC_002947/1460369‑1460665
                                                                                                                                                    |                                                                                                                                                     
tCCACTGCGGATCGAAGCTGCTGAGGTCGTTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCg                                                                                                                                                       <  1:305477/149‑1 (MQ=255)
     tGCGGATCGAAGCTGCTGAGGTCGTTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCg                                                                                                                                                    >  1:144308/1‑147 (MQ=255)
       cGGATCGAAGCTGCTGAGGTCGTTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGt                                                                                                                                                >  6:306370/1‑149 (MQ=255)
           tCGAAGCTGCTGAGGTCGTTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGc                                                                                                                                            <  6:19834/149‑1 (MQ=255)
             gAAGCTGCTGAGGTCGTTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCTGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCg                                                                                                                                           <  8:97217/148‑1 (MQ=255)
                   gctgAGGTCGTTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGc                                                                                                                                    >  5:187166/1‑149 (MQ=255)
                    ctgAGGTCGTTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCt                                                                                                                                   >  7:250927/1‑149 (MQ=255)
                               tttCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGcc                                                                                                                        <  2:101470/148‑1 (MQ=255)
                                 tCATTGGCAGATAACGTTCTGGGAATTCCCACACCAGCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGTGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCTGGGActg                                                                                                                      <  8:239935/148‑1 (MQ=255)
                                  cATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCctgc                                                                                                                     <  1:17840/148‑1 (MQ=255)
                                       gCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAg                                                                                                                <  3:202140/148‑1 (MQ=255)
                                        cAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTgc                                                                                                                                 >  2:171339/1‑131 (MQ=255)
                                        cAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTgc                                                                                                                                 <  1:171339/131‑1 (MQ=255)
                                         aGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCAGCAGGCCCTGC‑GCAgcga                                                                                                             <  5:306370/149‑2 (MQ=255)
                                                 ttCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAgg                                                                                                      >  1:190893/1‑148 (MQ=255)
                                                        aaTTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCa                                                                                              <  2:36338/149‑1 (MQ=255)
                                                          ttCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGt                                                                                            >  1:42295/1‑149 (MQ=255)
                                                            cccACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTg                                                                                           <  6:152889/148‑1 (MQ=255)
                                                             ccACACCACCACTTGTGGTGCGGCGTTCTTGACGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCg                                                                                         <  6:227109/149‑1 (MQ=255)
                                                             ccACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGGGTAGTTGGCGATGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCg                                                                                         >  3:60444/1‑149 (MQ=255)
                                                             ccACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGc                                                                                          <  7:207765/148‑1 (MQ=255)
                                                                   caccacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCctg                                                                                                                      <  1:125866/114‑1 (MQ=255)
                                                                   caccacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCctg                                                                                                                      >  2:125866/1‑114 (MQ=255)
                                                                   caccacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTgg                                                                                   >  8:261421/1‑149 (MQ=255)
                                                                      cacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGgcgc                                                                                >  3:24883/1‑149 (MQ=255)
                                                                         ttGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTg                                                                              >  7:36185/1‑148 (MQ=255)
                                                                         ttGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCAGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTg                                                                              >  3:313170/1‑148 (MQ=255)
                                                                         ttGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCAGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTg                                                                              >  3:308110/1‑148 (MQ=255)
                                                                           gtggtgCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAg                                                                           >  8:167240/1‑149 (MQ=255)
                                                                           gtggtgCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGCAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAg                                                                           >  2:251302/1‑149 (MQ=255)
                                                                              gtgCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTg                                                                        <  4:24883/149‑1 (MQ=255)
                                                                                 cggcgTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTg                                                                     >  8:182775/1‑149 (MQ=255)
                                                                                     gTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTa                                                                 >  3:13229/1‑149 (MQ=255)
                                                                                     gTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTa                                                                 <  6:187166/149‑1 (MQ=255)
                                                                                      ttCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTAc                                                                >  3:210441/1‑149 (MQ=255)
                                                                                         ttGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAg                                                             >  4:35089/1‑149 (MQ=255)
                                                                                             aGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGc                                                          >  7:177700/1‑148 (MQ=255)
                                                                                                cATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCGAGTGCGGGTTGGCGCTGTAGCTTGTGCCTACCAGTGCAAcc                                                      <  6:25239/149‑1 (MQ=255)
                                                                                                 aTCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTg                                                                     >  5:78168/1‑133 (MQ=255)
                                                                                                 aTCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTg                                                                     <  6:78168/133‑1 (MQ=255)
                                                                                                  tCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGGGTAGTTGGCGATGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAATGCAACCgg                                                    <  4:60444/149‑1 (MQ=255)
                                                                                                     cTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGAt                                                 >  3:172194/1‑149 (MQ=255)
                                                                                                      tttGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATc                                                >  3:33317/1‑149 (MQ=255)
                                                                                                          cAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGc                                             <  3:35089/148‑1 (MQ=255)
                                                                                                           aGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCtt                                           <  8:250927/149‑1 (MQ=255)
                                                                                                                 ttGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTCCCAGTGCAACCGGGATCTGCTTGTCGGc                                     >  8:245286/1‑149 (MQ=255)
                                                                                                                 ttGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGc                                     >  8:125477/1‑149 (MQ=255)
                                                                                                                 ttGAGCATCGGCAGCAGCGGGCCATGGCCGACTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTCCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGc                                     >  2:135064/1‑149 (MQ=255)
                                                                                                                   gAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGc                                     <  7:167240/147‑1 (MQ=255)
                                                                                                                        tCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGaacaac                              <  8:177700/149‑1 (MQ=255)
                                                                                                                          ggCAGCAGCGGGCCATGGCCGTCTTCAGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAAcgc                            <  4:313170/149‑1 (MQ=255)
                                                                                                                          ggCAGCAGCGGGCCATGGCCGTCTTCAGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAAcgc                            <  4:308110/149‑1 (MQ=255)
                                                                                                                           gcagcaGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAAcgcg                           >  7:288591/1‑149 (MQ=255)
                                                                                                                                 gcGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAAAAAAGCGGAGCCg                     >  6:225662/1‑149 (MQ=255)
                                                                                                                                   gggCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGc                   <  4:172194/149‑1 (MQ=255)
                                                                                                                                             cGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCg         >  2:166395/1‑149 (MQ=255)
                                                                                                                                             cGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCg         <  8:36185/149‑1 (MQ=255)
                                                                                                                                             cGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCg         >  2:257317/1‑149 (MQ=255)
                                                                                                                                               tCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGc       <  4:13229/149‑1 (MQ=255)
                                                                                                                                                 ttCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCAt     >  2:167949/1‑149 (MQ=255)
                                                                                                                                                    ggCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGa  <  4:210441/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCACTGCGGATCGAAGCTGCTGAGGTCGTTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGC‑GCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGA  >  NC_002947/1460369‑1460665

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑