Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,153,310 C→A 8.9% intergenic (‑1387/+27) 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,3100CA 8.9% 173.2 / 6.9 56intergenic (‑1387/+27)PP_5696/PP_4537hypothetical protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base A (2/3);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

GCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGTTCCAGATCA  >  NC_002947/5153166‑5153458
                                                                                                                                                |                                                                                                                                                    
gCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTg                                                                                                                                                   >  7:31956/1‑148 (MQ=255)
gCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTg                                                                                                                                                   <  7:220878/148‑1 (MQ=255)
gCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTg                                                                                                                                                   >  6:177584/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)
                                  agcagcCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGcc                                                                                                                                                      >  8:160297/1‑111 (MQ=255)
                                  agcagcCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGcc                                                                                                                                                      <  7:160297/111‑1 (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                                                                                                    <  2:14722/149‑1 (MQ=255)
                                              gCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGgcgc                                                                                                    >  5:262681/1‑149 (MQ=255)
                                              gCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGgcgc                                                                                                    <  6:262681/149‑1 (MQ=255)
                                                          gcgTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                                              <  8:352814/143‑1 (MQ=255)
                                                          gcgTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTgc                                                                                              >  7:352814/1‑143 (MQ=255)
                                                           cgTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGccca                                                                                       >  5:44621/1‑149 (MQ=255)
                                                            gTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGcccac                                                                                      >  7:377859/1‑149 (MQ=255)
                                                            gTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGcccac                                                                                      <  5:57408/149‑1 (MQ=255)
                                                               gggTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGa                                                                                                                             <  2:104622/107‑1 (MQ=255)
                                                               gggTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGa                                                                                                                             >  1:104622/1‑107 (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)
                                                                                                                                                |                                                                                                                                                    
GCGGTGTTGCCGGAGTTCGGCTTGAGCCTGTATCAGCAGCCGAGTGGCAATGACATGCGCGTGGGGTTGGCCGGGCGCACGGCGGCAGAGCCGATGCCGCGTCGCATTGAAGATATGAGTGAGGCTTGAGATTGCCGGGGCTGCCCTGCAGCCCCCGGATCAAAGGGCGAGTCAGGCCTTGCCGGCCAGCGGCGCCTGTGCGCGCCCACCCAGGCGGCTCAGCCACACCGACGCCAGGATGATCCCCCCGCCGATGAACAGCCGCCCCAACGGCTCGTGCTGGTTCCAGATCA  >  NC_002947/5153166‑5153458

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