Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,547,696 G→A 7.2% A342T (GCA→ACA)  murP → N‑acetylmuramic acid permease, EIIBC component, PTS system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,547,6960GA 7.2% 151.0 / 2.9 60A342T (GCA→ACA) murPN‑acetylmuramic acid permease, EIIBC component, PTS system
Reads supporting (aligned to +/- strand):  ref base G (32/22);  new base A (2/2);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

GAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCC  >  NC_000913/2547551‑2547839
                                                                                                                                                 |                                                                                                                                               
gAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAg                                                                                                                                              <  1:943399/149‑1 (MQ=255)
gAACAGTAATCCGTTCGGTTGTGCGGTCTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCTACGTTTCATTggagg                                                                                                                                              >  2:25435/1‑146 (MQ=255)
      tAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACCGCTTATTACCCATCCTTTCAATGGCAGGCGCgg                                                                                                                                        >  1:690704/1‑149 (MQ=255)
      tAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGGGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCgg                                                                                                                                        <  1:167738/149‑1 (MQ=255)
        aTCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGc                                                                                                                                      >  1:354401/1‑149 (MQ=255)
           cGTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAg                                                                                                                                   >  1:651591/1‑149 (MQ=255)
            gTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAgg                                                                                                                                  <  1:75941/149‑1 (MQ=255)
              tCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTg                                                                                                                                <  2:243773/149‑1 (MQ=255)
                 gTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGggc                                                                                                                             >  2:824999/1‑149 (MQ=255)
                    gtgCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGggcgcg                                                                                                                          <  1:358556/149‑1 (MQ=255)
                         gTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGggcgcggcgat                                                                                                                     >  1:1078808/1‑147 (MQ=255)
                           tttAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAAAAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTgg                                                                                                                   >  1:175396/1‑149 (MQ=255)
                             ccgCCGGGCTGTGCCGGATCGCCGGCGTCTTTCGCGTGCATCAGGGCCTTATGCCTGAAGACCTCGCGTTAAGGCACAGCCAGGGATTCAACAGCGTATTTCCCATCCTTTCAATGGCACGCGCGGGCCAGGTGGGCGCGGCGCTGGCa                                                                                                                 <  1:222116/147‑1 (MQ=255)
                                ccGGGCTGTTCCTGATCGCCGTGGTTTTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCActc                                                                                                              <  2:422321/149‑1 (MQ=255)
                                     cTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCCCTCTACTg                                                                                                         >  1:544478/1‑149 (MQ=255)
                                         tCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGAAGGCGCGGGCCAGGTGGGCGCGGCGCTGTCACTCTACTggcgg                                                                                                     >  1:10022/1‑149 (MQ=255)
                                           cTGATCGGCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTTTTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTggcgggc                                                                                                   <  2:264841/149‑1 (MQ=255)
                                           cTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGACAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTggcgggc                                                                                                   >  1:1024724/1‑149 (MQ=255)
                                           cTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGACAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTggcgggc                                                                                                   >  1:1029411/1‑149 (MQ=255)
                                                 gCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGGTTACCTCGCGTTAATGGCCAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCCATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACGCTACTGGCGGGCGCAAcc                                                                                             >  2:112239/1‑149 (MQ=255)
                                                  ccGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACCGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGAGCGGCGCTGGCACTCTACTGGCGGGCGCAACCg                                                                                            >  1:693631/1‑149 (MQ=255)
                                                  ccGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCACCAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCg                                                                                            >  1:972471/1‑149 (MQ=255)
                                                    gTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGACAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGca                                                                                          <  2:1024724/149‑1 (MQ=255)
                                                    gTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGACAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGca                                                                                          <  2:1029411/149‑1 (MQ=255)
                                                    gTGGTATTTGGCGTGCAGCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGca                                                                                          <  1:300137/149‑1 (MQ=255)
                                                         aTTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTGAATGGACAGCCAGGGATTCAACAGCTTATTTCCAATCCTTTCAATGGCAGGCGCGGGCCAGTTGGTCGCTGCGCTGGCACTCAACTGGCGGCCGCAACCGCCCAGtc                                                                                     >  1:227526/1‑148 (MQ=255)
                                                                gTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACg                                                                              >  2:1064702/1‑149 (MQ=255)
                                                                   cATCAGGGCTTTATTCCTGTATACCTCGCGCTAATGGACAGCCAGGGAGTCAACAGCTTATTTCCCATCCTTTCAATGCCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAg                                                                           <  1:487440/149‑1 (MQ=255)
                                                                           cTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTAc                                                                   >  1:1042856/1‑149 (MQ=255)
                                                                              tATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCATTGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTAcgcg                                                                <  2:399813/149‑1 (MQ=255)
                                                                               aTTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCgg                                                               >  1:61370/1‑149 (MQ=255)
                                                                                 tCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCgggg                                                             >  1:968892/1‑149 (MQ=255)
                                                                                  ccTGTTTACCTCGCGTTAATGGACAGCCAGGGATACAACAGCTTATTTACCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGc                                                            >  1:1079525/1‑149 (MQ=255)
                                                                                    tGTTTACCTCGCGTTAATGGACAGCCAGGGATTAAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCTCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGa                                                          >  1:1050534/1‑149 (MQ=255)
                                                                                          ccTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTc                                                    >  1:282517/1‑149 (MQ=255)
                                                                                              gcgTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCgg                                                <  2:1042856/149‑1 (MQ=255)
                                                                                               cgTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGc                                               >  2:99608/1‑149 (MQ=255)
                                                                                               cgTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGc                                               <  1:26432/149‑1 (MQ=255)
                                                                                                      ggacccccAGGGAGTCAACAGCTTAATTCCCACCCTTTCAAGGGCAGGCGCGCGCCAGGGGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTgg                                        <  1:84994/143‑1 (MQ=255)
                                                                                                      ggACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTgg                                        <  2:1109914/149‑1 (MQ=255)
                                                                                                      ggACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTgg                                        <  2:121063/149‑1 (MQ=255)
                                                                                                            ccAGGGATTCAACAGGTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCCCTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTACCGGCCTGCTGGGCGtgg                                  >  2:372455/1‑147 (MQ=255)
                                                                                                            ccAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTg                                  >  1:1017574/1‑149 (MQ=255)
                                                                                                            ccAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTg                                  >  2:137405/1‑149 (MQ=255)
                                                                                                            ccAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTg                                  >  1:803089/1‑149 (MQ=255)
                                                                                                               gggATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTg                               <  2:354401/149‑1 (MQ=255)
                                                                                                                     cAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGc                         >  1:476415/1‑149 (MQ=255)
                                                                                                                     cAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGc                         >  1:621419/1‑149 (MQ=255)
                                                                                                                      aaCAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCt                        >  2:665512/1‑149 (MQ=255)
                                                                                                                        cAGCTTTTTTCCCATCCTTTCAATGGCTGGCGCGGGCCAGGTGTGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGa                      <  2:439336/149‑1 (MQ=255)
                                                                                                                        cAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGa                      >  2:237830/1‑149 (MQ=255)
                                                                                                                        cAGCTTATTTCCCATCCTTTCAATGCCACGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCATCCGCACAGTGCGCTTCGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGa                      <  1:126740/149‑1 (MQ=255)
                                                                                                                             ttattCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGTCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTAc                 <  2:690704/146‑1 (MQ=255)
                                                                                                                                ttCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGgt              >  2:741991/1‑149 (MQ=255)
                                                                                                                                   ccATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGTTGAACCGCTGATTTACGGTGTc           >  2:409299/1‑149 (MQ=255)
                                                                                                                                   ccATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTc           <  1:395480/149‑1 (MQ=255)
                                                                                                                                    cATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCa          >  1:42641/1‑149 (MQ=255)
                                                                                                                                     aTCCTTTGAATGGCAGGCGAGGTCCAGGTGGGCGCGGCGCTGGCACTCTACTGTCGGGCGCAACCGCAGAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCTCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCAc         <  2:227526/149‑1 (MQ=255)
                                                                                                                                      tCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCAcc        <  2:112938/149‑1 (MQ=255)
                                                                                                                                        cTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCt      >  1:1036815/1‑149 (MQ=255)
                                                                                                                                            cAATGGAAGGCGCGGGCGAGGCGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGccc  <  1:1093822/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGCTGGCACTCTACTGGCGGGCGCAACCGCACAGTGCGCTACGCAGTCAGGTACGCGGGGCGATTATTCCCGGCCTGCTGGGCGTTGGTGAACCGCTGATTTACGGTGTCACCCTGCCC  >  NC_000913/2547551‑2547839

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