Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,304,058 T→C 6.5% F232L (TTC→CTC)  cysS → cysteine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,304,0580TC 6.5% 176.9 / 2.2 62F232L (TTC→CTC) cysScysteine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base C (2/2);  total (31/31)
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.63e-01

CTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGTTGAACAACTTCTTCACCATCCGCGAC  >  NC_002947/3303911‑3304201
                                                                                                                                                   |                                                                                                                                                
cTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGtt                                                                                                                                                 >  1:227300/1‑149 (MQ=255)
   tgGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTccc                                                                                                                                              >  3:40927/1‑149 (MQ=255)
   tgGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTccc                                                                                                                                              <  4:267346/149‑1 (MQ=255)
   tgGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTccc                                                                                                                                              <  7:151467/149‑1 (MQ=255)
     gAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGc                                                                                                                                            >  7:207374/1‑149 (MQ=255)
       aGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGcac                                                                                                                                         >  6:234414/1‑149 (MQ=255)
        gggCGACAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGcacc                                                                                                                                        <  4:342682/149‑1 (MQ=255)
        gggCGACAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGcacc                                                                                                                                        <  4:348085/149‑1 (MQ=255)
            gTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCacga                                                                                                                                    <  7:327896/149‑1 (MQ=255)
            gTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCacga                                                                                                                                    <  4:40927/149‑1 (MQ=255)
                 gCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTAATGTCCACCTGCTGCCTGGGTGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCacgagaac                                                                                                                                >  5:311660/1‑148 (MQ=255)
                    gggCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCacgagaacgaga                                                                                                                            <  7:157399/149‑1 (MQ=255)
                    gggCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTACCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCacgagaacgaga                                                                                                                            <  7:323292/149‑1 (MQ=255)
                           ccAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGcc                                                                                                                       >  2:111935/1‑147 (MQ=255)
                           ccAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCa                                                                                                                     <  5:234414/149‑1 (MQ=255)
                             aaGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCaga                                                                                                                   <  3:298485/149‑1 (MQ=255)
                                cTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCg                                                                                                                <  8:264064/149‑1 (MQ=255)
                                           gtgccagCCCGGGCGACCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCAcc                                                                                                     >  3:281580/7‑149 (MQ=255)
                                                   ccGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCa                                                                                             <  1:111935/149‑1 (MQ=255)
                                                        tCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACg                                                                                        <  3:242660/149‑1 (MQ=255)
                                                           cccGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGcc                                                                                      >  4:329001/1‑148 (MQ=255)
                                                                  tGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTg                                                                              >  5:247175/1‑149 (MQ=255)
                                                                  tGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTg                                                                              >  8:23570/1‑149 (MQ=255)
                                                                  tGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGCTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTg                                                                              >  2:343219/1‑149 (MQ=255)
                                                                  tGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGCTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTg                                                                              >  2:348669/1‑149 (MQ=255)
                                                                     cacaTCGAGTTCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTTGCGGCAGCGACCTGGATGTCCCGCACCACGAGAACGAGATTTCCCAGAGCGAGGCTGCCACCGGCAAGCAGTACGCCAACGCCTGGAt                                                                           <  6:356360/149‑1 (MQ=255)
                                                                      acaTCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATg                                                                          >  5:245826/1‑149 (MQ=255)
                                                                         tCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCAGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCg                                                                                                    <  2:369217/120‑1 (MQ=255)
                                                                         tCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCAGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCg                                                                                                    >  1:369217/1‑120 (MQ=255)
                                                                            aGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGc                                                                    >  7:182434/1‑149 (MQ=255)
                                                                                cTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCa                                                                >  8:18804/1‑149 (MQ=255)
                                                                                 tCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCAt                                                               <  6:245826/149‑1 (MQ=255)
                                                                                 tCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCAt                                                               >  7:266881/1‑149 (MQ=255)
                                                                                   ggTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATg                                                              <  7:23570/148‑1 (MQ=255)
                                                                                       aTGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCg                                                         >  3:329850/1‑149 (MQ=255)
                                                                                            cACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATc                                                     >  5:15549/1‑148 (MQ=255)
                                                                                            cACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATc                                                     >  3:162150/1‑148 (MQ=255)
                                                                                               ctgctgCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACg                                                 <  2:227300/149‑1 (MQ=255)
                                                                                                tgctgcCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACgg                                                >  2:307432/1‑149 (MQ=255)
                                                                                                         ggCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGagaag                                         <  3:329001/147‑1 (MQ=255)
                                                                                                         ggCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGAt                                       <  4:162150/149‑1 (MQ=255)
                                                                                                         ggCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGAt                                       <  7:18804/149‑1 (MQ=255)
                                                                                                                gCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGCTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCaa                                 <  1:348669/148‑1 (MQ=255)
                                                                                                                gCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGCTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCaa                                 <  1:343219/148‑1 (MQ=255)
                                                                                                                 cTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCaa                                 <  4:355369/147‑1 (MQ=255)
                                                                                                                 cTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCaa                                 >  3:355369/1‑147 (MQ=255)
                                                                                                                  ttCGACATTCACGGTGGCGGCAGCGCCCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTc                              >  8:230025/1‑149 (MQ=255)
                                                                                                                  ttCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTc                              >  7:310388/1‑149 (MQ=255)
                                                                                                                     gACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGtt                           <  8:297304/149‑1 (MQ=255)
                                                                                                                     gACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGtt                           <  8:266881/149‑1 (MQ=255)
                                                                                                                      aCATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCTACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGtt                           >  1:246770/1‑149 (MQ=255)
                                                                                                                        aTTCACGGTGGCGGCAGCGACCTTGATTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCaa                                 <  7:53654/140‑1 (MQ=255)
                                                                                                                        aTTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGt                                                        <  6:271591/117‑1 (MQ=255)
                                                                                                                        aTTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGt                                                        >  5:271591/1‑117 (MQ=255)
                                                                                                                        aTTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGc                                               >  7:2390/1‑126 (MQ=255)
                                                                                                                        aTTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGc                                               <  8:2390/126‑1 (MQ=255)
                                                                                                                            aCGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCCCCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGTTGaacaac                    >  8:65043/1‑149 (MQ=255)
                                                                                                                               gTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACAGGCAAGCAGTACGCCAACGCCTGGATCCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGTTGAACAActtc                 >  1:19502/1‑149 (MQ=255)
                                                                                                                                     gcAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGTTGAACAACTTCTTCAcc           <  1:307432/149‑1 (MQ=255)
                                                                                                                                     gcAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGTTGAACAACTTCTTCAc            >  4:66212/1‑148 (MQ=255)
                                                                                                                                            ccTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGTTGAATAACTTCTTCACCATCcgc     >  3:335138/1‑148 (MQ=255)
                                                                                                                                              tGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGTTGAACAACTTCTTCACCATCCGCGAc  >  4:226963/1‑149 (MQ=255)
                                                                                                                                                gAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGTTGAACAACTTCTTCAGCATCCGCGAc  <  2:19502/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGCCAACGCCTGGATGCATTGCGGCATGATCCGTATCAACGGCGAGAAGATGTCCAAGTCGTTGAACAACTTCTTCACCATCCGCGAC  >  NC_002947/3303911‑3304201

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: