Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,234,539 C→G 9.7% R298R (CGC→CGG) ‡ PP_4612 → fec operon regulator FecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,234,5390CG 9.7% 183.5 / 8.9 62R298R (CGC→CGG) ‡PP_4612fec operon regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base G (2/4);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TCGGGCGCTATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTGCGCCCCAAC  >  NC_002947/5234394‑5234684
                                                                                                                                                 |                                                                                                                                                 
tCGGGCGCTATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTAc                                                                                                                                                >  2:188768/1‑149 (MQ=255)
  gggCGCTATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTg                                                                                                                                              >  5:59914/1‑149 (MQ=255)
   ggCGCTTTCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTTTCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTgg                                                                                                                                             <  7:25382/149‑1 (MQ=255)
   ggCGCTATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTgg                                                                                                                                             <  4:30866/149‑1 (MQ=255)
   ggCGCTATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTgg                                                                                                                                             <  4:34744/149‑1 (MQ=255)
       cTATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGa                                                                                                                                         >  1:283647/1‑149 (MQ=255)
        tATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCAGTAGCGGGTGaa                                                                                                                                        >  4:263198/1‑148 (MQ=255)
         aTCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCAGTAGCGGGTGaac                                                                                                                                       <  3:222358/149‑3 (MQ=255)
          tCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCAGTAGCGGGTGaact                                                                                                                                      <  3:263198/149‑4 (MQ=255)
          tCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCAGTAGCGGGTGaac                                                                                                                                       <  1:269740/148‑3 (MQ=255)
            gCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGt                                                                                                                                     <  6:226740/148‑1 (MQ=255)
             ccAGGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCg                                                                                                                                   >  4:135010/1‑149 (MQ=255)
              cATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGa                                                                                                                                  <  8:97625/149‑1 (MQ=255)
              cATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGa                                                                                                                                  <  6:107761/149‑1 (MQ=255)
               aTGGGCAGTTGGGGTTCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTACAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCAGTAGCGGGTGaac                                                                                                                                       <  1:245045/143‑3 (MQ=255)
               aTGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTACAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCAGTAGCGGGTGaac                                                                                                                                       >  2:245045/1‑141 (MQ=255)
                  ggCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGccc                                                                                                                              <  3:156696/149‑1 (MQ=255)
                         ggggTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATc                                                                                                                        <  6:59914/148‑1 (MQ=255)
                          gggTGCAGCTAGCTGTTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGTCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTTACGGTCGAGCCCAAAATCAg                                                                                                                      <  4:95547/149‑1 (MQ=255)
                          gggTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCa                                                                                                                       <  1:188768/148‑1 (MQ=255)
                               cAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAgc                                                                                                                 >  1:276902/1‑149 (MQ=255)
                                 gCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAgcgc                                                                                                               <  2:144176/149‑1 (MQ=255)
                                    aGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAg                                                                                                            >  3:235048/1‑149 (MQ=255)
                                    aGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAg                                                                                                            >  1:188272/1‑149 (MQ=255)
                                         gTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGATAGCGCCAGCAAAg                                                                                                       >  4:113687/1‑149 (MQ=255)
                                             ccGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCaaaa                                                                                                                          <  1:169019/126‑1 (MQ=255)
                                             ccGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCaaaa                                                                                                                          >  2:169019/1‑126 (MQ=255)
                                                 cTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTTAGCTCATCCTGCATCTGCTGGATGTGGCGTTGCCGGTGCCTGGGAAGCAGTTCACCCTCTACTGTGTGAAGGTCGAGCCCTTAATCAGTTAGCGCCAGCAAAGGCGTAat                                                                                                <  2:286636/148‑1 (MQ=255)
                                                    gcgcATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATAtt                                                                                             <  1:5419/148‑1 (MQ=255)
                                                      gcATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTc                                                                                          <  2:276902/149‑1 (MQ=255)
                                                      gcATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTc                                                                                          <  3:135010/149‑1 (MQ=255)
                                                      gcATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTc                                                                                          <  3:113687/149‑1 (MQ=255)
                                                      gcATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTc                                                                                          <  4:160829/149‑1 (MQ=255)
                                                      gcATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACCGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTc                                                                                          <  3:35662/149‑1 (MQ=255)
                                                        aTTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCat                                                                                        <  2:188272/149‑1 (MQ=255)
                                                              gggTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGa                                                                                   <  4:235048/148‑1 (MQ=255)
                                                                 tCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGa                                                                               <  2:139652/149‑1 (MQ=255)
                                                                   gTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGc                                                                             >  8:26218/1‑149 (MQ=255)
                                                                       ccGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCaaa                                                                                                        <  3:37845/118‑1 (MQ=255)
                                                                       ccGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCaaa                                                                                                        >  4:37845/1‑118 (MQ=255)
                                                                         gTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGGCAgttt                                                                       >  2:125600/1‑149 (MQ=255)
                                                                         gTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCAtttt                                                                       >  1:178529/1‑149 (MQ=255)
                                                                          ttGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTc                                                                      >  5:76367/1‑149 (MQ=255)
                                                                             gATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTccc                                                                    >  4:183463/1‑148 (MQ=255)
                                                                                  cAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTc                                                              <  2:293554/149‑1 (MQ=255)
                                                                                     tGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCACTTAGAGCCAGCAAAGGCGAAATATTTTCATATCTACAGAGCCATTTTCACAACCtaccg                                                           >  4:18448/1‑145 (MQ=255)
                                                                                                 gCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGCTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACgg                                               >  5:215860/1‑149 (MQ=255)
                                                                                                     aTGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAg                                            >  8:95112/1‑148 (MQ=255)
                                                                                                        cTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAg                                         <  2:283647/148‑1 (MQ=255)
                                                                                                                       ccGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGgacc                          <  7:95112/148‑1 (MQ=255)
                                                                                                                         ggTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACt                                                                                 >  2:160040/1‑91 (MQ=255)
                                                                                                                         ggTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACt                                                                                 <  1:160040/91‑1 (MQ=255)
                                                                                                                                gtgAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGCAGCCCTCTCTCAGGACCCGCCCCATGc                >  1:212235/1‑149 (MQ=255)
                                                                                                                                gagAAGCAGTTACCCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCAAGCAAAGGCGAAATATTTTCATAACGAATGACCCATTTTCCCAACCTCGCGTGACAGACACTGCAGACGCCCTCTCTCAGGACCCGACCCAAGc                >  2:127045/3‑149 (MQ=255)
                                                                                                                                  gAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGCAGCGCTCTCTCAGGACCCGACCCATGCCg              >  4:191468/1‑149 (MQ=255)
                                                                                                                                  gAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCg              >  2:155622/1‑149 (MQ=255)
                                                                                                                                        gTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTgcgc        <  6:76367/149‑1 (MQ=255)
                                                                                                                                        gTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTgcgc        <  3:183463/149‑1 (MQ=255)
                                                                                                                                           cACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTGCGccc      >  6:200651/1‑148 (MQ=255)
                                                                                                                                            aCCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTGCGcccc     >  2:238461/1‑148 (MQ=255)
                                                                                                                                            aCCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGAGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCACTCTCTCAGGCCCCGACCCATGCCTTTGCGCCCCa    >  8:280864/1‑149 (MQ=255)
                                                                                                                                              ccGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTGCGCCCCAAc  >  5:174387/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TCGGGCGCTATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTGCGCCCCAAC  >  NC_002947/5234394‑5234684

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