Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,412,629 G→A 6.1% G105S (GGC→AGC)  PP_1234 → permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,412,6290GA 6.1% 196.9 / 2.5 65G105S (GGC→AGC) PP_1234permease
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base A (2/2);  total (32/33)
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.68e-01

TCAACGCCCTGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGGCGCTGACCTTCTCGC  >  NC_002947/1412483‑1412773
                                                                                                                                                   |                                                                                                                                                
tCAACGCCCTGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCggg                                                                                                                                                 >  1:270402/1‑148 (MQ=255)
tCAACGCCCTGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCggg                                                                                                                                                 >  6:98343/1‑148 (MQ=255)
tCAACGCCCTGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGc                                                                                                                                                <  5:79221/149‑1 (MQ=255)
     gCCCTGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCt                                                                                                                                           >  6:243907/1‑149 (MQ=255)
       ccTGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTc                                                                                                                                          <  5:106662/148‑1 (MQ=255)
       ccTGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTc                                                                                                                                          >  6:106662/1‑148 (MQ=255)
       ccTGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCg                                                                                                                                         <  6:195571/149‑1 (MQ=255)
         tGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTAATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCgg                                                                                                                                        >  1:257245/1‑148 (MQ=255)
           gAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGc                                                                                                                                       >  2:190345/1‑147 (MQ=255)
             gcgcgtgcgcgtACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCATGCCCGGCAGCTCgttgaac                                                                                                                                   <  4:369885/149‑7 (MQ=255)
               gcgtgcgcgtACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTgg                                                                                                                                 <  5:243907/149‑1 (MQ=255)
                 gtgcgcgtACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCAGGCATGCTCGGCAAGTgg                                                                                                                                 >  1:51803/1‑147 (MQ=255)
                        tACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCgttg                                                                                                                        <  5:98343/149‑1 (MQ=255)
                         aCCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCgttgt                                                                                                                       >  3:323898/1‑149 (MQ=255)
                          ccTACACGGCTGGAGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCgttgt                                                                                                                       <  1:137411/148‑1 (MQ=255)
                           cTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGttgttg                                                                                                                     <  8:25549/149‑1 (MQ=255)
                              cacGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTATTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTtgctg                                                                                                                  >  2:35641/1‑149 (MQ=255)
                                cGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTtgctgct                                                                                                                <  1:190345/149‑1 (MQ=255)
                                 ggCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTtgctgctg                                                                                                               <  7:215354/149‑1 (MQ=255)
                                      gCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGcc                                                                                                          <  4:191987/149‑1 (MQ=255)
                                                       cGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAg                                                                                                      >  3:56048/1‑136 (MQ=255)
                                                       cGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAg                                                                                                      <  4:56048/136‑1 (MQ=255)
                                                                cATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCg                                                                                 <  2:257245/148‑1 (MQ=255)
                                                                  tGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCAGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGAc                                                                              <  2:51803/149‑1 (MQ=255)
                                                                             gCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTAtt                                                                   >  7:111046/1‑149 (MQ=255)
                                                                                 tgtTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCAc                                                               <  4:323898/149‑1 (MQ=255)
                                                                                     cATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCa                                                           >  4:286015/1‑149 (MQ=255)
                                                                                     cATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCa                                                           >  4:315119/1‑149 (MQ=255)
                                                                                      aTGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCAt                                                          >  7:260380/1‑149 (MQ=255)
                                                                                       tgctggtgctggTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCAtt                                                         <  1:35641/149‑1 (MQ=255)
                                                                                            gtgctggtACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGt                                                     <  8:111046/148‑1 (MQ=255)
                                                                                                tggtACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGTTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTg                                                >  7:53868/1‑149 (MQ=255)
                                                                                                tggtACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTg                                                >  6:320094/1‑149 (MQ=255)
                                                                                                  gtACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTa                                                                                                 >  5:337743/1‑98 (MQ=255)
                                                                                                  gtACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTa                                                                                                 <  6:337743/98‑1 (MQ=255)
                                                                                                  gtACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACCAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGAAGGTATTGCACGCCATTGAGTCGCTgcg                                              >  2:351157/1‑149 (MQ=255)
                                                                                                   tACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTgcgc                                             <  3:315119/149‑1 (MQ=255)
                                                                                                     ccGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCAGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGcgcg                                            <  5:320094/148‑1 (MQ=255)
                                                                                                      cGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGa                                                                               <  8:147500/112‑1 (MQ=255)
                                                                                                      cGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGa                                                                               >  7:147500/1‑112 (MQ=255)
                                                                                                      cGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTg                                                                  >  5:292881/1‑125 (MQ=255)
                                                                                                      cGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTg                                                                  <  6:292881/125‑1 (MQ=255)
                                                                                                      cGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGAATTGCACGCCAtt                                                         >  5:315969/1‑134 (MQ=255)
                                                                                                      cGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGAATTGCACGCCAtt                                                         <  6:315969/134‑1 (MQ=255)
                                                                                                            tgtgGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTg                                                                                       >  5:305639/1‑98 (MQ=255)
                                                                                                            tgtgGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTg                                                                                       <  6:305639/98‑1 (MQ=255)
                                                                                                                  aCCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAg                              >  4:100647/1‑149 (MQ=255)
                                                                                                                          aTCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGAc                                                                              <  4:76492/93‑1 (MQ=255)
                                                                                                                          aTCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGAc                                                                              >  3:76492/1‑93 (MQ=255)
                                                                                                                            cACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAAt                    >  1:60979/1‑149 (MQ=255)
                                                                                                                            cACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAAt                    >  3:259485/1‑149 (MQ=255)
                                                                                                                              cGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATgg                  <  2:270402/149‑1 (MQ=255)
                                                                                                                               gCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCAtt                                                         >  5:148789/1‑109 (MQ=255)
                                                                                                                               gCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCAtt                                                         <  6:148789/109‑1 (MQ=255)
                                                                                                                                 tGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGgcg               >  1:191010/1‑149 (MQ=255)
                                                                                                                                  gTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGgcgc              >  5:135450/1‑149 (MQ=255)
                                                                                                                                     cAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGGCGCTg            <  2:60979/148‑1 (MQ=255)
                                                                                                                                      aaCGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCa                                >  1:333247/1‑127 (MQ=255)
                                                                                                                                      aaCGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCa                                <  2:333247/127‑1 (MQ=255)
                                                                                                                                        cGAGCTGCCGGTCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGgcgc              <  1:29616/143‑1 (MQ=255)
                                                                                                                                        cGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGgcgc              >  2:29616/1‑143 (MQ=255)
                                                                                                                                            cttGCCGAGCATGCCCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGGCGCTg            >  5:128268/3‑142 (MQ=255)
                                                                                                                                            cttGCCGAGCATGCCCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTCCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGGCGCTg            <  6:128268/140‑1 (MQ=255)
                                                                                                                                            cTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGGCGCTGACCTtctc    <  4:259485/149‑1 (MQ=255)
                                                                                                                                             ttGCCGAGCATGCCCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGgc                <  4:176548/136‑1 (MQ=255)
                                                                                                                                             ttGCCGAGCATGCCCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGgc                >  3:176548/2‑137 (MQ=255)
                                                                                                                                               gCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGGCGCTGACCTTCTCGc  >  7:238022/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCAACGCCCTGGAGCGCGTGCGCGTACCTACACGGCTGGCGGTGATGCTGGTGTTCGCCTTGTTCATGGGCGCACTGGCGGTGTTCATGCTGGTGCTGGTACCGCTGCTGTGGCACCAGCTGATCACGCTGTTCAACGAGCTGCCGGGCATGCTCGGCAAGTGGCAGTCGTTGTTGCTGCTGCTGCCGGAGCGCTATCCGCATCTGGTGTCGGACGAGCAGGTATTGCACGCCATTGAGTCGGTGCGCGGTGAAATTGGCAAGTTTGGTCAATGGGCGCTGACCTTCTCGC  >  NC_002947/1412483‑1412773

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: