Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,234,543 T→C 11.1% W300R (TGG→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,5430TC11.1% 164.8 / 9.6 63W300R (TGG→CGG) PP_4612fec operon regulator FecR
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base C (2/5);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGCTATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTGCGCCCCAACCCGTTC  >  NC_002947/5234396‑5234690
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                <  3:135010/149‑1 (MQ=255)
                                                    gcATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTc                                                                                                <  3:113687/149‑1 (MQ=255)
                                                    gcATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTc                                                                                                <  4:160829/149‑1 (MQ=255)
                                                    gcATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTc                                                                                                <  2:276902/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                                                                                       <  1:160040/91‑1 (MQ=255)
                                                                                                                       ggTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACt                                                                                       >  2:160040/1‑91 (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)
                                                                                                                                                tACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTGCGCCCCAAcccat    >  6:257991/1‑147 (MQ=255)
                                                                                                                                                  cTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTGCGCCCCAACCcgttc  <  1:155622/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGCGCTATCGCCATGGGCAGTTGGGGTGCAGCGAGCAGGTTGCCGACTTGCGCATTTCTGGGTCGTATCCGTTGGATGACAGTGAGCGCATCCTGCAAATGCTGGAAGTGGCGTTGCCGGTGCGTGTGAAGCAGTTCACCCGCTACTGGGTGACGGTCGAGCCCAAAATCAGTTAGCGCCAGCAAAGGCGAAATATTTTCATATCGACTGAGCCATTTTCCCAACCTCGCGTGACAGACACGGCAGAAGCCCTCTCTCAGGACCCGACCCATGCCGTTGCGCCCCAACCCGTTC  >  NC_002947/5234396‑5234690

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