Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,153,312 T→G 8.3% intergenic (‑1389/+25) PP_5696 ← / ← PP_4537 hypothetical protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,153,3120TG 8.3% 172.8 / 7.9 60intergenic (‑1389/+25)PP_5696/PP_4537hypothetical protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base G (2/3);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-01

GCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGTTCCAGATCAA  >  NC_002947/5153166‑5153459
                                                                                                                                                  |                                                                                                                                                   
gCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTg                                                                                                                                                    >  6:177584/1‑148 (MQ=255)
gCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTg                                                                                                                                                    <  7:220878/148‑1 (MQ=255)
gCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTg                                                                                                                                                    >  7:31956/1‑148 (MQ=255)
   gtgtTGCAGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGc                                                                                                                                                >  8:346913/1‑149 (MQ=255)
      ttGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCAGGGCAGcccc                                                                                                                                             <  6:198131/149‑1 (MQ=255)
      ttGCCGGAGTTCGGCTTGAGCCTGTATAAGCAGACGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGACGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGcccc                                                                                                                                             <  8:169550/149‑1 (MQ=255)
          cGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGa                                                                                                                                         <  2:48821/149‑1 (MQ=255)
          cGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGa                                                                                                                                         <  5:153165/149‑1 (MQ=255)
           ggAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCGGGGGCTGCAGGGCAGccccggcaa                                                                                                                                        <  6:115230/149‑6 (MQ=17)
           ggAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCAGGGCAGccccggcaa                                                                                                                                        <  7:354763/149‑6 (MQ=255)
              gTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCa                                                                                                                                      >  6:227589/1‑148 (MQ=255)
               ttCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCaaa                                                                                                                                    <  5:151138/149‑1 (MQ=255)
                 cGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAgg                                                                                                                                  >  5:330947/1‑149 (MQ=255)
                   gCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGc                                                                                                                                >  1:14722/1‑149 (MQ=255)
                     ttGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGa                                                                                                                              <  3:115965/149‑1 (MQ=255)
                       gAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGa                                                                                                                              <  6:330947/147‑1 (MQ=255)
                             gTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGcc                                                                                                                      >  6:391453/1‑149 (MQ=255)
                                tCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTg                                                                                                                   >  5:55422/1‑149 (MQ=255)
                                  agcagcCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGa                                                                                                                              <  5:18921/136‑1 (MQ=255)
                                  agcagcCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGa                                                                                                                              >  6:18921/1‑136 (MQ=255)
                                  agcagcCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGcc                                                                                                                 <  6:49782/149‑1 (MQ=255)
                                        cGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAg                                                                                                           <  7:346913/149‑1 (MQ=255)
                                          aGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAgcg                                                                                                         >  5:296233/1‑149 (MQ=255)
                                              gCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGgcgc                                                                                                     >  5:262681/1‑149 (MQ=255)
                                              gCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGgcgc                                                                                                     <  6:262681/149‑1 (MQ=255)
                                              gCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGgcgc                                                                                                     <  2:14722/149‑1 (MQ=255)
                                                          gcgTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                                               >  7:352814/1‑143 (MQ=255)
                                                          gcgTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                                               <  8:352814/143‑1 (MQ=255)
                                                           cgTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGccca                                                                                        >  5:44621/1‑149 (MQ=255)
                                                            gTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGcccac                                                                                       <  5:57408/149‑1 (MQ=255)
                                                            gTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGcccac                                                                                       >  7:377859/1‑149 (MQ=255)
                                                               gggTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGa                                                                                                                              >  1:104622/1‑107 (MQ=255)
                                                               gggTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGa                                                                                                                              <  2:104622/107‑1 (MQ=255)
                                                               gggTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGcccaccca                                                                                    <  5:391453/149‑1 (MQ=255)
                                                                  ttGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAggc                                                                                 <  6:296233/149‑1 (MQ=255)
                                                                    ggCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAggcgg                                                                               <  5:59906/149‑1 (MQ=255)
                                                                    ggCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAggcgg                                                                               >  6:239464/1‑149 (MQ=255)
                                                                    ggCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAggcg                                                                                <  8:125527/148‑1 (MQ=255)
                                                                        gggCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCa                                                                           <  6:44621/149‑1 (MQ=255)
                                                                         ggCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCCGCCCCCGGATGAAAGGGGGAGTCAGGCCATGCCGGCCCGCGGCGCCTGTGCGCGCCCCCCCCGGCGGCtgcg                                                                          >  3:219231/1‑146 (MQ=255)
                                                                          gcgcACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGc                                                                         <  5:227589/149‑1 (MQ=255)
                                                                          gcgcACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTACAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGc                                                                         >  5:373386/1‑149 (MQ=255)
                                                                               cggcggCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTACAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCacac                                                                    <  6:373386/149‑1 (MQ=255)
                                                                                       gagCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCa                                                            <  5:239464/149‑1 (MQ=255)
                                                                                               gCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATc                                                    >  8:52129/1‑149 (MQ=255)
                                                                                                        cATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCccgccg                                           <  7:52129/149‑1 (MQ=255)
                                                                                                         aTTGAAGATATGAGTGAGGGTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCCCCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGa                                          >  5:109262/1‑149 (MQ=255)
                                                                                                                        gagGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCa                           >  3:83065/1‑149 (MQ=255)
                                                                                                                           gCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACg                        <  4:83065/149‑1 (MQ=255)
                                                                                                                                agaTTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTc                    <  3:209067/148‑1 (MQ=255)
                                                                                                                                 gaTTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGt                  >  3:123937/1‑149 (MQ=255)
                                                                                                                                  aTTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTg                 >  7:362530/1‑149 (MQ=255)
                                                                                                                                       gggggCTGCAGGGCAG‑CCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGtt           >  5:197241/2‑149 (MQ=255)
                                                                                                                                        ggggCTGCAGGGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGtt           >  2:98664/1‑149 (MQ=255)
                                                                                                                                         gggCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGTTc          <  8:362530/149‑1 (MQ=255)
                                                                                                                                             tGCCCTGCAGCCCCCGGATCAAAGGGCCAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGTTc          >  1:18909/1‑145 (MQ=255)
                                                                                                                                               cccTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGTCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGTTCCAGa      <  8:377859/147‑1 (MQ=255)
                                                                                                                                                ccTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGTTCCAGATCa   <  8:210257/149‑1 (MQ=255)
                                                                                                                                                ccTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGACGAACAGCCGCCCCAACGGCTCGTGCTGGTTCCAGATc    >  1:160331/1‑148 (MQ=255)
                                                                                                                                                 cTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGTTCCAGATCaa  <  5:177584/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGTTCCAGATCAA  >  NC_002947/5153166‑5153459

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: ‑