Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,454,275 T→C 7.1% N296S (AAC→AGC)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,454,2750TC 7.1% 164.4 / 2.5 56N296S (AAC→AGC) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base C (2/2);  total (28/28)
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.06e-01

CCAGCTCTTCGCTGCTGACCACCACCCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGGTAGTACTCGATACCGGCGGGGATGATGCCGCTGGCCAGGTC  >  NC_002947/2454132‑2454423
                                                                                                                                               |                                                                                                                                                    
ccAGCTCTTCGCTGCTGACCA‑CACCCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGaa                                                                                                                                               >  1:762656/1‑150 (MQ=255)
 cAGCTCTTCGCTGCTGACCACCACCCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGaa                                                                                                                                               >  2:89898/1‑150 (MQ=255)
  aGCTCTTCGCTGCTGACCACCACCCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAg                                                                                                                                              >  2:949312/1‑150 (MQ=255)
   gCTCTTCGCTGCTGACCACCACCCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGt                                                                                                                                             <  2:1527813/150‑1 (MQ=255)
      ctTCGCTGCTGACCACCACCCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCt                                                                                                                                          <  1:1547729/150‑1 (MQ=255)
             gctgACCACCACCCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCggcggc                                                                                                                                   >  1:1261915/1‑150 (MQ=255)
             gctgACCACCACCCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCggcggc                                                                                                                                   <  1:1321664/150‑1 (MQ=255)
                 accaccaccCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGt                                                                                                                               >  2:858855/1‑150 (MQ=255)
                   caccaccCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTc                                                                                                                             <  1:89898/150‑1 (MQ=255)
                   caccaccCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTc                                                                                                                             >  2:1527079/1‑150 (MQ=255)
                   caccaccCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTc                                                                                                                             >  2:1332602/1‑150 (MQ=255)
                         ccGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACAcc                                                                                                                       >  2:1133885/1‑150 (MQ=255)
                          cGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCg                                                                                                                      <  2:1563080/150‑1 (MQ=255)
                             ggCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGt                                                                                                                   <  1:9213/150‑1 (MQ=255)
                                     cTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGaa                                                                                                           <  2:1261915/150‑1 (MQ=255)
                                          gCCGGGCGAAGCAATCTTCGACCGGCAGGGCTAGCCCGGGCGGCGGCAACAGCGGCCGGCGGAGGCCGCCGCGCCGGTCGCCATAACCCCCGCGCACGGCGTTCCAGTAGTGCTGGGCGGCCTGTTCCACACCGGGTAGCGGGACCCCCt                                                                                                      >  2:367997/1‑150 (MQ=255)
                                                  aaGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTgtcg                                                                                              <  1:1133885/150‑1 (MQ=255)
                                                      aaTCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACATGGGAgtcggtcg                                                                                          >  1:1157380/1‑150 (MQ=255)
                                                       aTCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGtcggtcgg                                                                                         >  1:1295707/1‑150 (MQ=255)
                                                        tCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGc                                                                                        <  1:673817/150‑1 (MQ=255)
                                                          ttCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAg                                                                                      <  1:88804/150‑1 (MQ=255)
                                                             cACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCCTGCGGGCGTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTa                                                                                   <  1:1350789/150‑1 (MQ=255)
                                                              aCCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGaa                                                                                                           >  2:629545/1‑125 (MQ=255)
                                                              aCCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGaa                                                                                                           <  1:629545/125‑1 (MQ=255)
                                                               ccGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGt                                                                                 >  1:35674/1‑150 (MQ=255)
                                                                cGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTc                                                                                                                             >  2:1240753/1‑105 (MQ=255)
                                                                cGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTc                                                                                                                             <  1:1240753/105‑1 (MQ=255)
                                                                     ggAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACa                                                                           <  2:1526845/150‑1 (MQ=255)
                                                                      gAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTa                                                                                                                  >  2:517504/1‑110 (MQ=255)
                                                                      gAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTa                                                                                                                  <  1:517504/110‑1 (MQ=255)
                                                                            cTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACa                                                                           <  1:537526/143‑1 (MQ=255)
                                                                            cTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACa                                                                           >  2:537526/1‑143 (MQ=255)
                                                                            cTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGa                                                                    >  1:546654/1‑150 (MQ=255)
                                                                            cTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGa                                                                    <  1:1631368/150‑1 (MQ=255)
                                                                             tCGGCCGGCGGCAACAGCGGCCGGCCGCGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAg                                                                   >  1:458262/1‑150 (MQ=255)
                                                                                    gcggcAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTttctt                                                            >  1:1535462/1‑150 (MQ=255)
                                                                                    gcggcAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTttctt                                                            <  1:363525/150‑1 (MQ=255)
                                                                                        cAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGTTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCgaa                                                        <  1:1332602/150‑1 (MQ=255)
                                                                                         tgttgccgccggccgAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCgaag                                                       >  1:170744/8‑150 (MQ=255)
                                                                                               ggccggcTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGTGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGt                                                 <  2:1064732/150‑1 (MQ=255)
                                                                                               ggccggcTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGt                                                 >  2:281918/1‑150 (MQ=255)
                                                                                                    gcTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAg                                            >  1:148417/1‑150 (MQ=255)
                                                                                                     cTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAgg                                           <  2:1295707/150‑1 (MQ=255)
                                                                                                          gTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGgtagta                                      >  2:87696/1‑150 (MQ=255)
                                                                                                          gTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGCTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTttcttc                                                           <  1:1489795/129‑1 (MQ=255)
                                                                                                          gTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGCTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTttcttc                                                           <  1:1489678/129‑1 (MQ=255)
                                                                                                          gTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGCTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTttcttc                                                           >  2:1489795/1‑129 (MQ=255)
                                                                                                          gTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGCTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTttcttc                                                           >  2:1489678/1‑129 (MQ=255)
                                                                                                                 cggcggTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGATGTTTCTTCGAAGAACAGTGGCAGGTAGTACTCGATa                               <  2:236720/150‑1 (MQ=255)
                                                                                                                    cggTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGGTAGTACTCGATACcg                            <  2:546654/150‑1 (MQ=255)
                                                                                                                        cTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACATGGGTGTCGGTCGGCAGGTAGCCGAACAGCCTCGACGTATCTTTGTAGAACGGTGTCATTTAGTACTCGATACcggcgg                        <  2:1008175/150‑1 (MQ=255)
                                                                                                                         ttCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGACCAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGGTCGTACTCGCTACCGGCggg                       >  2:733084/1‑150 (MQ=255)
                                                                                                                         ttCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGGTAGTACTCGATACCGGCggg                       >  2:461480/1‑150 (MQ=255)
                                                                                                                            aTAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGGTAGTACTCGATACCGGCGGGgat                    <  2:170744/150‑1 (MQ=255)
                                                                                                                            aTAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGGTAGTACTCGATACCGGCGGGgat                    >  1:434327/1‑150 (MQ=255)
                                                                                                                                 gCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGGTAGTACTCGATACCGGCGGGGATGATGc               <  1:461480/150‑1 (MQ=255)
                                                                                                                                           gtcgtTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTgtcggtc                                                                                           >  1:1326412/1‑64 (MQ=255)
                                                                                                                                           gtcgtTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTgtcggtc                                                                                           <  2:1326412/64‑1 (MQ=255)
                                                                                                                                              gtTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGt                                                                                 <  2:1360405/71‑1 (MQ=255)
                                                                                                                                              gtTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGt                                                                                 >  1:1360405/1‑71 (MQ=255)
                                                                                                                                              gtTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGGCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGGTAGTACTCGATACCGGCGGGGATGATGCCGCTGGCCAgggc  >  1:716759/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCAGCTCTTCGCTGCTGACCACCACCCGGGGCCACTGCTTGAGCCGGGCGAAGCAATCTTCCACCGGCAGGAACAGCTCGGCCGGCGGCAACAGCGGCCGGCTGAGGTCGCCGCGGCGGTCTTCATAACGCCCGCGCACGTCGTTCCAGAAGTGCTCGGCGGCCTGTTCCACACCGGGTAGCGAGAACACCTGGGTGTCGGTCGGCAGGTAGTCGAACAGGGTCGAGGTTTCTTCGAAGAACAGTGGCAGGTAGTACTCGATACCGGCGGGGATGATGCCGCTGGCCAGGTC  >  NC_002947/2454132‑2454423

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