Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,902,226 G→C 6.5% R926P (CGC→CCC)  PP_1703 → assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,902,2260GC 6.5% 200.2 / 3.8 61R926P (CGC→CCC) PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base C (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CGACAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAACCTGCACCC  >  NC_002947/1902089‑1902371
                                                                                                                                         |                                                                                                                                                 
cGACAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGc                                                                                                                                        >  3:249516/1‑149 (MQ=255)
 gACAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCt                                                                                                                                       >  4:96419/1‑149 (MQ=255)
  aCAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTg                                                                                                                                      <  4:196759/149‑1 (MQ=255)
     gcggcgcggcCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGa                                                                                                                                   <  3:96419/149‑1 (MQ=255)
          gcggcCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCa                                                                                                                              >  6:33154/1‑149 (MQ=255)
                ttCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCct                                                                                                                        >  7:222193/1‑149 (MQ=255)
                       aGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATc                                                                                                                            >  8:245887/1‑138 (MQ=255)
                       aGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATc                                                                                                                            <  7:245887/138‑1 (MQ=255)
                             tGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGt                                                                                                            <  6:291117/148‑1 (MQ=255)
                              gCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTg                                                                                                           >  6:178923/1‑148 (MQ=255)
                               cAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGa                                                                                                         <  4:249516/149‑1 (MQ=255)
                                aGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGAc                                                                                                        <  2:180231/149‑1 (MQ=255)
                                 gggCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACt                                                                                                       >  6:293132/1‑149 (MQ=255)
                                 gggCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACt                                                                                                       >  6:317464/1‑149 (MQ=255)
                                  ggCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTg                                                                                                      >  5:53023/1‑149 (MQ=255)
                                   gCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGc                                                                                                     <  5:237556/149‑1 (MQ=255)
                                      aggaggGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAg                                                                                                  >  1:5680/1‑149 (MQ=255)
                                         aggGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGccc                                                                                               <  7:147723/149‑1 (MQ=255)
                                              gcTCACTGCTCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGTCTCCAGTTAGGACCAGTTCTGCGGTTTTGGCCGCAAGGTCGACCAGCGCCTGGCCGAGCTGGGATGGCATCGCCTGCTGGGGCGGGTGGACTGTGAGCCCGACtt                                                                                          <  7:59139/149‑1 (MQ=255)
                                                  aCTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGGCCAGGCGCTGGTCGAGc                                                                                                                                        <  3:180857/99‑1 (MQ=255)
                                                  aCTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGGCCAGGCGCTGGTCGAGc                                                                                                                                        >  4:180857/1‑99 (MQ=255)
                                                  aCTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGCCCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTcgac                                                                                      >  5:114619/1‑149 (MQ=255)
                                                  aCTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTcgac                                                                                      <  5:178923/149‑1 (MQ=255)
                                                       gcCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGc                                                                                 >  6:20420/1‑149 (MQ=255)
                                                             cTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTc                                                                                         <  1:183686/135‑1 (MQ=255)
                                                             cTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTc                                                                                         >  2:183686/1‑135 (MQ=255)
                                                                  cTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTgctggct                                                                      <  2:5680/149‑1 (MQ=255)
                                                                         gTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCAGTTTTGGCCGCAAGCTCGACCAGCCCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGccc                                                                                               >  1:226110/1‑117 (MQ=255)
                                                                         gTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCAGTTTTGGCCGCAAGCTCGACCAGCCCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGccc                                                                                               <  2:226110/117‑1 (MQ=255)
                                                                               gCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGaa                                                              <  6:209832/144‑1 (MQ=255)
                                                                               gCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGaa                                                              >  5:209832/1‑144 (MQ=255)
                                                                               gCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGt                                                          <  7:102644/148‑1 (MQ=255)
                                                                                     ggCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGAACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGc                                                    <  3:312034/148‑1 (MQ=255)
                                                                                      gCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCt                                                  >  3:2881/1‑149 (MQ=255)
                                                                                      gCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCt                                                  <  5:293132/149‑1 (MQ=255)
                                                                                           tCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGt                                                                                                            >  4:92882/1‑86 (MQ=255)
                                                                                           tCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACt                                             >  7:301117/1‑149 (MQ=255)
                                                                                             cAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTgg                                           >  8:30153/1‑149 (MQ=255)
                                                                                             cAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTCGGAGGGCATCGCCTGCTGGCGCGGGTg                                                                                                           >  1:126928/1‑85 (MQ=255)
                                                                                              aGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTggg                                          >  4:145296/1‑149 (MQ=255)
                                                                                                ttACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCa                                        >  4:315403/1‑149 (MQ=255)
                                                                                                          ttCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCgcggc                              >  4:99386/1‑149 (MQ=255)
                                                                                                             tGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGgcgc                            >  7:179704/1‑148 (MQ=255)
                                                                                                              gCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCa                          >  1:48760/1‑149 (MQ=255)
                                                                                                               cGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGcc                           <  5:33154/147‑1 (MQ=255)
                                                                                                               cGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGcc                           <  6:53023/147‑1 (MQ=255)
                                                                                                                       gCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGc                  <  3:99386/148‑1 (MQ=255)
                                                                                                                        ccGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGc                  <  1:11969/147‑1 (MQ=255)
                                                                                                                        ccGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGc                  >  2:11969/1‑147 (MQ=255)
                                                                                                                        ccGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGa                <  4:2881/149‑1 (MQ=255)
                                                                                                                         cGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCg                                                                <  6:216326/100‑1 (MQ=255)
                                                                                                                         cGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCg                                                                >  5:216326/1‑100 (MQ=255)
                                                                                                                              gCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAAc          <  3:145296/149‑1 (MQ=255)
                                                                                                                              gCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAAc          <  5:317464/149‑1 (MQ=255)
                                                                                                                                tCGACCAGCGCCTGGCCGAGCTGGGCGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGCCTTCGACGAAGCCTTTGCTGGCTGGCTGGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAACCt        >  3:5100/1‑149 (MQ=255)
                                                                                                                                tCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCACGACTTCGAAGAAGCATTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGAACTGCCGCGGCGCCATATCCGGCGAGTGAACCt        >  6:243494/1‑149 (MQ=255)
                                                                                                                                  gACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAACCt        <  3:315403/147‑1 (MQ=255)
                                                                                                                                   aCCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGc                                                            <  5:168604/94‑1 (MQ=255)
                                                                                                                                   aCCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGc                                                            >  6:168604/1‑94 (MQ=255)
                                                                                                                                      aGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAACCTGCAccc  <  2:48760/149‑1 (MQ=255)
                                                                                                                                      aGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAACCTGCAccc  >  8:204301/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
CGACAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGAAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAACCTGCACCC  >  NC_002947/1902089‑1902371

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