Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 999,625 C→G 7.0% G165A (GGC→GCC)  PP_0861 ← outer membrane ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947999,6250CG 7.0% 181.6 / 4.0 57G165A (GGC→GCC) PP_0861outer membrane ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base G (2/2);  total (27/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 = 8.07e-01

GGCGACCGGCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACGGGCGGT  >  NC_002947/999479‑999769
                                                                                                                                                  |                                                                                                                                                
ggCAACCGGCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGgccgcgccg                                                                                                                                                >  8:270484/1‑149 (MQ=255)
  cGACCGGCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGccgcgccgcg                                                                                                                                              >  5:321226/1‑149 (MQ=255)
       ggCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGGCGCGGccgg                                                                                                                                         <  6:47007/149‑3 (MQ=255)
        gCAACGGTATCGCTGAACTCGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCCGCCAcc                                                                                                                                        >  1:103668/1‑149 (MQ=255)
             ggTATCGCTGAACTCGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATgg                                                                                                                                   >  2:234515/1‑149 (MQ=255)
                        aCTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCgg                                                                                                                         <  3:308834/148‑1 (MQ=255)
                            gTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGccc                                                                                                                      <  5:71883/147‑1 (MQ=255)
                            gTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCt                                                                                                                    >  1:26016/1‑149 (MQ=255)
                            gTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCt                                                                                                                    >  3:163168/1‑149 (MQ=255)
                               gTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAACTCGGccatcgg                                                                                                                  <  6:133142/148‑6 (MQ=255)
                                ttGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGc                                                                                                                <  3:159075/149‑1 (MQ=255)
                                     aTCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACt                                                                                                            >  2:98705/1‑148 (MQ=255)
                                        gAAGGTGTAGCGCTGTGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCaa                                                                                                        <  8:264169/149‑1 (MQ=255)
                                        gAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCaa                                                                                                        <  1:231884/149‑1 (MQ=255)
                                         aaGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCaa                                                                                                        >  8:159958/1‑148 (MQ=255)
                                          aGGTGTAGCTCTGGGTCTGGTCGCTGCCCCCGGTCCAGGCGCCGTCCTGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGGTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACa                                                                                                      <  6:309693/149‑1 (MQ=255)
                                             tgtAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGAt                                                                                                   <  4:163168/149‑1 (MQ=255)
                                               tAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGAtt                                                                                                  <  1:234515/148‑1 (MQ=255)
                                                     tGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTg                                                                                                                  >  3:192443/1‑126 (MQ=255)
                                                     tGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTg                                                                                                                  <  4:192443/126‑1 (MQ=255)
                                                      gggTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCg                                                                                          <  7:270484/149‑1 (MQ=255)
                                                        gTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGaa                                                                                        <  1:85895/149‑1 (MQ=255)
                                                        gTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGaa                                                                                        >  4:78955/1‑149 (MQ=255)
                                                           tGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCg                                                                                                                          <  5:62279/112‑1 (MQ=255)
                                                           tGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCg                                                                                                                          >  6:62279/1‑112 (MQ=255)
                                                           tGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGAt                                                                                     >  8:206488/1‑149 (MQ=255)
                                                                      ccAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGttt                                                                          >  6:101231/1‑149 (MQ=255)
                                                                      ccAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGttt                                                                          <  2:26016/149‑1 (MQ=255)
                                                                          gTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGgc                                                                      >  4:285377/1‑149 (MQ=255)
                                                                            ccAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGgcgc                                                                    <  3:11126/149‑1 (MQ=255)
                                                                             cAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGgcgc                                                                    >  3:159239/1‑148 (MQ=255)
                                                                             cAGGCGCCGTCCAGCGAGTTCCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGgcgc                                                                    <  3:6662/148‑1 (MQ=255)
                                                                                    cGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTca                                                            <  4:159239/149‑1 (MQ=255)
                                                                                     gTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTcac                                                           >  4:27805/1‑149 (MQ=255)
                                                                                      tCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATgg                                                                                            <  8:132338/115‑1 (MQ=255)
                                                                                      tCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATgg                                                                                            >  7:132338/1‑115 (MQ=255)
                                                                                        cAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTcacgca                                                        >  7:150510/1‑149 (MQ=255)
                                                                                        cAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTcacgca                                                        <  3:27805/149‑1 (MQ=255)
                                                                                          gCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGACGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTcacgcacg                                                      <  6:321226/149‑1 (MQ=255)
                                                                                              gTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTgg                                                                        >  5:220433/1‑127 (MQ=255)
                                                                                              gTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTgg                                                                        <  6:220433/127‑1 (MQ=255)
                                                                                                  cccAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATgg                                                                                                                                   >  8:43484/1‑64 (MQ=255)
                                                                                                   ccAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTc                                                             >  6:90654/1‑133 (MQ=255)
                                                                                                   ccAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTc                                                             <  5:90654/133‑1 (MQ=255)
                                                                                                    cAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCgg                                                                              >  2:153629/1‑115 (MQ=255)
                                                                                                    cAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCgg                                                                              <  1:153629/115‑1 (MQ=255)
                                                                                                     aGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGt                                            >  2:275937/1‑148 (MQ=255)
                                                                                                                 ttgcttACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTg                               >  1:12262/1‑149 (MQ=255)
                                                                                                                 ttgcttACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTg                               >  5:213722/1‑149 (MQ=255)
                                                                                                                   gcttACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTgg                              <  3:285377/148‑1 (MQ=255)
                                                                                                                   gcttACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTgg                              <  7:206488/148‑1 (MQ=255)
                                                                                                                        ccAGGTTGATGGTACCGCCGGCCGCGCCGCAGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTc                         <  7:159958/148‑1 (MQ=255)
                                                                                                                                aTGGTACCGGCGGCCGCGCCGCGGCCACGGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCcgc                 <  2:103668/148‑1 (MQ=255)
                                                                                                                                 tGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCgg               >  7:161180/1‑149 (MQ=255)
                                                                                                                                     aCCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATg            >  4:213203/1‑148 (MQ=255)
                                                                                                                                         gtggccgCGGCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGaa                      >  7:226801/3‑134 (MQ=255)
                                                                                                                                         gtggccgCGGCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGaa                      <  8:226801/132‑1 (MQ=255)
                                                                                                                                         gtggccgCGGCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACgg       >  8:145918/3‑149 (MQ=255)
                                                                                                                                          cggccgCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACggg      <  1:275937/149‑1 (MQ=255)
                                                                                                                                              cgcgccgcgGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACGGGCGGt  >  1:247874/1‑149 (MQ=255)
                                                                                                                                              cgcgccgcgGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACGGGCGGt  >  5:28176/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGCGACCGGCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACGGGCGGT  >  NC_002947/999479‑999769

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