Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,917,195 T→A 7.0% V2072E (GTG→GAG)  PP_2561 → hemolysin‑type calcium‑binding bacteriocin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,917,1950TA 7.0% 183.6 / 3.4 67V2072E (GTG→GAG) PP_2561hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base A (3/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

AGCGGTTGTTCCAGGCGGCCAAGTTCGGCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGACCAGCAACTGGGGCTGATCGCGGCCTTCCT  >  NC_002947/2917048‑2917339
                                                                                                                                                   |                                                                                                                                                
aGCGGTTGTTCCAGGCGGCCAAGTTCGGCAGCGAAATGCCGTACCAGCTCGTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGGCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTg                                                                                                                                                 <  1:362259/149‑1 (MQ=255)
aGCGGTTGTTCCAGGCGGCCAAGTTCGGCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTg                                                                                                                                                 <  8:208656/149‑1 (MQ=255)
          ccagccagCAAAGTTCGGAACCGAAATACAGAACCAGCACATCGTCTTCTACGAGTTAGCCCCAACCCTGCAACCACGGCTAGCCATGTTCTAAGCCCAGACCCAGCTATACCACGTCGACCTCGATGCCTCGATCGATCCCCAGTTCg                                                                                                                                       >  8:103978/1‑149 (MQ=255)
          ccAGGCGGCCAAGTTCGGCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCg                                                                                                                                       >  2:113498/1‑149 (MQ=255)
           cAGGCGGCCATGTTCGGCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGc                                                                                                                                      <  7:103978/149‑1 (MQ=255)
           cAGGCGGCCAAGTTCGGCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGc                                                                                                                                      >  7:82310/1‑149 (MQ=255)
                    aaGTTCGGCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATTGTGGCCGAGTTCGCCCACACCg                                                                                                                              >  6:88006/1‑148 (MQ=255)
                       ttCGGCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGtata                                                                                                                          >  6:107433/1‑149 (MQ=255)
                           gCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGc                                                                                                                      >  3:109983/1‑149 (MQ=255)
                           gCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGc                                                                                                                      >  3:263629/1‑149 (MQ=255)
                            cACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCACGATCGAGGcatc                                                                                                                                            >  7:52759/1‑123 (MQ=255)
                            cACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCACGATCGAGGcatc                                                                                                                                            >  7:54515/1‑123 (MQ=255)
                            cACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCACGATCGAGGcatc                                                                                                                                            <  8:52759/126‑4 (MQ=255)
                            cACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCACGATCGAGGcatc                                                                                                                                            <  8:54515/126‑4 (MQ=255)
                             aCCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCtt                                                                                                                    >  7:180036/1‑149 (MQ=255)
                                 aaaTGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGc                                                                                                                >  2:174678/1‑149 (MQ=255)
                                 aaaTGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGc                                                                                                                <  8:82310/149‑1 (MQ=255)
                                    tGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTg                                                                                                                                                 >  4:313205/1‑113 (MQ=255)
                                    tGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTg                                                                                                                                                 <  3:313205/113‑1 (MQ=255)
                                     gCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACt                                                                                                            >  1:62848/1‑149 (MQ=255)
                                       aGTACCAGCACCTGGAGTTCGAGGATTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATAGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTc                                                                                                           >  8:335731/1‑148 (MQ=255)
                                         tACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGAt                                                                                                        >  6:57305/1‑149 (MQ=255)
                                           ccAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATTGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATg                                                                                                       <  5:88006/148‑1 (MQ=255)
                                              gCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGa                                                                                                   >  2:160517/1‑149 (MQ=255)
                                                aCCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGAcc                                                                                                 <  5:57305/149‑1 (MQ=255)
                                                  cTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGAcc                                                                                                 >  4:175086/1‑147 (MQ=255)
                                                  cTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCga                                                                                               <  1:160517/149‑1 (MQ=255)
                                                    ggTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGCCCTCGCTGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCgaga                                                                                             >  7:159688/1‑149 (MQ=255)
                                                     gtgtTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGAc                                                                                            >  2:273687/1‑149 (MQ=255)
                                                          tgaggagTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTc                                                                                        <  7:60648/147‑1 (MQ=255)
                                                                 ttCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGTTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCtt                                                                                >  2:98528/1‑149 (MQ=255)
                                                                 ttCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCtt                                                                                >  3:289815/1‑149 (MQ=255)
                                                                      ccGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGta                                                                                                                            >  1:210996/1‑100 (MQ=255)
                                                                      ccGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGta                                                                                                                            <  2:210996/100‑1 (MQ=255)
                                                                       cGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACAt                                                                          <  3:314440/149‑1 (MQ=255)
                                                                       cGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACAt                                                                          >  4:314440/1‑149 (MQ=255)
                                                                          aCCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGa                                                                       <  2:68197/149‑1 (MQ=255)
                                                                            cGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGTTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACt                                                                     <  1:98528/149‑1 (MQ=255)
                                                                                   ccTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTc                                                                                                           >  2:15549/1‑104 (MQ=255)
                                                                                   ccTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTc                                                                                                           <  1:15549/104‑1 (MQ=255)
                                                                                       gggTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCa                                                            <  3:175086/147‑1 (MQ=255)
                                                                                         gTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCg                                                                                                >  1:302636/1‑109 (MQ=255)
                                                                                         gTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCg                                                                                                <  2:302636/109‑1 (MQ=255)
                                                                                            gACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGa                                                     <  2:62848/149‑1 (MQ=255)
                                                                                             aCCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGAc                                                    <  5:107433/149‑1 (MQ=255)
                                                                                                  ttCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACGGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCgg                                                <  4:289815/148‑1 (MQ=255)
                                                                                                  ttCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTc                                                             <  2:266233/135‑1 (MQ=255)
                                                                                                  ttCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTc                                                             >  1:266233/1‑135 (MQ=255)
                                                                                                  ttCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCg                                                 <  8:180036/147‑1 (MQ=255)
                                                                                                     tttGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGc                                                                                                                                               <  2:97127/50‑1 (MQ=255)
                                                                                                     tttGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGc                                                                                                                                               >  1:97127/1‑50 (MQ=255)
                                                                                                      ttGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGCCTTCAACGTCATCCAGGACCCGGCCAgc                                           >  6:141638/1‑149 (MQ=255)
                                                                                                          cccGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGcc                                        >  4:337904/1‑148 (MQ=255)
                                                                                                           ccGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCaa                                      >  7:27040/1‑149 (MQ=255)
                                                                                                               cccAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAAcccc                                  <  1:273687/149‑1 (MQ=255)
                                                                                                                ccAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCg                                 <  4:109983/149‑1 (MQ=255)
                                                                                                                ccAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCg                                 <  4:263629/149‑1 (MQ=255)
                                                                                                                  aGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGAc                               <  1:113498/149‑1 (MQ=255)
                                                                                                                    gTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGACCa                             <  7:335731/149‑1 (MQ=255)
                                                                                                                       taCGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGt                                                              <  8:321409/113‑1 (MQ=255)
                                                                                                                       taCGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGt                                                              >  7:321409/1‑113 (MQ=255)
                                                                                                                           acgTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGACCAGCAACTg                      >  6:292889/1‑149 (MQ=255)
                                                                                                                              tCGACCTCGCTTCCTCCATCGTGGCCTCGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGACCAGCAACTgggg                   <  8:159688/149‑1 (MQ=39)
                                                                                                                                 aCCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGACCAGCAACTGGGGCt                 >  7:322092/1‑148 (MQ=255)
                                                                                                                                      gATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGACCAGCAACTGGGGCTGATcgc           >  5:17836/1‑149 (MQ=255)
                                                                                                                                              gATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGACCAGCAACTGGGGCTGATCGCGGCCTTcc   >  7:146986/1‑149 (MQ=255)
                                                                                                                                               aTCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGACCAGCAACTGGGGCTGATCGCGGCCTTCCt  <  5:292889/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AGCGGTTGTTCCAGGCGGCCAAGTTCGGCACCGAAATGCAGTACCAGCACCTGGTGTTCGAGGAGTTCGCCCGCACCGTGCAACCTCGGGTCGACCTGTTCTTTGCCCCGACCCAGGTATACGACGTCGACCTCGATGCCTCGATCGTGGCCGAGTTCGCCCACACCGTATACCGCTTCGGCCACTCGATGCTGACCGAGACTGTCGACCGCTTCGACATCGACTTCAACGTCATCCAGGACCCGGCCAGCGCCAACCCCGACCAGCAACTGGGGCTGATCGCGGCCTTCCT  >  NC_002947/2917048‑2917339

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: