Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,622,188 G→A 6.8% E368K (GAG→AAG)  PP_1421 → sensor histidine kinase TctE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,622,1880GA 6.8% 170.8 / 2.0 59E368K (GAG→AAG) PP_1421sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (2/2);  total (30/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.61e-01

GCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCG  >  NC_002947/1622045‑1622336
                                                                                                                                                |                                                                                                                                                    
gCGGTGGGCAGCGTCTGGACTTGAGTCAGTTGGCCGGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGTGGCCGAGGCGCCGGTGTGGTTGAAGCGGGGAGCCGACATTGCTCAACGAgtt                                                                                                                                                  <  1:3375/149‑1 (MQ=255)
    tGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCGGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGa                                                                                                                                             <  4:90632/149‑1 (MQ=255)
    tGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGGGGCGCTGGCGCGGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGa                                                                                                                                             <  6:302209/149‑1 (MQ=255)
        cAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCaa                                                                                                                                         >  5:185629/1‑149 (MQ=255)
                 gACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGa                                                                                                                                <  5:40096/149‑1 (MQ=255)
                    ttGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACAACTCGTTGAGCaat                                                                                                                                        >  5:125703/1‑137 (MQ=255)
                    ttGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACAACTCGTTGAGCaat                                                                                                                                        <  6:125703/138‑2 (MQ=255)
                      gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGTCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATg                                                                                                                           >  6:73831/1‑149 (MQ=255)
                      gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGAc                                                                                                                               >  5:179415/1‑145 (MQ=255)
                      gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGAc                                                                                                                               <  6:179415/145‑1 (MQ=255)
                      gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACaa                                                                                                                             <  5:224150/147‑1 (MQ=255)
                      gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAAt                                                                                                                            >  5:215893/1‑148 (MQ=255)
                      gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAAt                                                                                                                            <  6:215893/148‑1 (MQ=255)
                      gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATg                                                                                                                           >  3:166306/1‑149 (MQ=255)
                        gTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGcc                                                                                                                         >  3:129419/1‑149 (MQ=255)
                         tCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTAGATCAAG‑GTGGACCCCAAATAAATAAAAAAAAAATTGAGCAACC‑GGTGGACAATGCCtt                                                                                                                       >  2:98/1‑149 (MQ=255)
                            gTTGGGCGGCGTGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTg                                                                                                                      <  4:166306/148‑1 (MQ=255)
                                gCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGgcgc                                                                                                                  <  4:129419/148‑1 (MQ=255)
                                  ccGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGa                                                                                                                                                     >  2:268468/1‑111 (MQ=255)
                                  ccGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGa                                                                                                                                                     <  1:268468/111‑1 (MQ=255)
                                          tGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCg                                                                                                       >  1:206679/1‑149 (MQ=255)
                                            ggCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGTAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACAAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGTATGCGCCaaggca                                                                                                     <  1:72106/149‑6 (MQ=255)
                                            ggCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCTAGGCGCCTGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCggg                                                                                                     >  5:212919/1‑149 (MQ=255)
                                            ggCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCggg                                                                                                     >  5:136096/1‑149 (MQ=255)
                                            ggCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCggg                                                                                                     >  6:222948/1‑149 (MQ=255)
                                              catggccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTg                                                                                                   >  1:131961/1‑149 (MQ=255)
                                                tggccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGc                                                                                                 <  4:257937/149‑1 (MQ=255)
                                                   ccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAt                                                                                              <  2:206679/149‑1 (MQ=255)
                                                            cTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTg                                                                                                   >  7:99345/1‑135 (MQ=255)
                                                            cTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTg                                                                                                   <  8:99345/135‑1 (MQ=255)
                                                            cTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTg                                                                                     >  6:215923/1‑149 (MQ=255)
                                                                 gcgcACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCggg                                                                                 <  5:215923/148‑1 (MQ=255)
                                                                    cacaAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAAt                                                                             <  5:222948/149‑1 (MQ=255)
                                                                                gTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGt                                                                 <  5:259745/149‑1 (MQ=255)
                                                                                    cgctggcgctggAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtg                                                             <  5:73831/149‑1 (MQ=255)
                                                                                          cgctggAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACC‑GGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtgg                                                      <  1:98/149‑1 (MQ=255)
                                                                                           gctggAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtgg                                                      <  6:185629/149‑1 (MQ=255)
                                                                                           gctggAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtg                                                       >  4:63484/1‑148 (MQ=255)
                                                                                                     gAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACTACgg                                            >  7:256551/1‑149 (MQ=255)
                                                                                                          gcCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAg                                       >  7:245202/1‑149 (MQ=255)
                                                                                                          gcCGGTGTGGGGTTAG‑GGTGCGCCGACATTGATCAACGAGTTGTGGAGCAACCTGGTGGACACTGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGAGTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAg                                       <  8:256551/149‑1 (MQ=255)
                                                                                                            cGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtgg                                                      >  3:220887/1‑132 (MQ=255)
                                                                                                            cGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtgg                                                      <  4:220887/132‑1 (MQ=255)
                                                                                                                 tgGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTcg                                >  2:295241/1‑148 (MQ=255)
                                                                                                                   gTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTccc                               <  3:63484/148‑1 (MQ=255)
                                                                                                                   gTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCg                              >  6:234471/1‑149 (MQ=255)
                                                                                                                      gAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtgg                                                                                                                                    <  8:102898/44‑1 (MQ=255)
                                                                                                                      gAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtgg                                                                                                                                    >  7:102898/1‑44 (MQ=255)
                                                                                                                       aaG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGc                          >  3:252778/1‑149 (MQ=255)
                                                                                                                         g‑ggggAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCg                        >  4:103581/1‑149 (MQ=255)
                                                                                                                             ggAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                     <  6:136096/149‑1 (MQ=255)
                                                                                                                             ggAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                     <  3:149052/149‑1 (MQ=255)
                                                                                                                             ggAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                     <  3:144454/149‑1 (MQ=255)
                                                                                                                                     cATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACgg             >  1:142458/1‑149 (MQ=255)
                                                                                                                                           tCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc       >  4:123080/1‑149 (MQ=255)
                                                                                                                                           tCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGAAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc       >  5:42929/1‑149 (MQ=255)
                                                                                                                                           tCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGAAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc       >  5:38392/1‑149 (MQ=255)
                                                                                                                                             aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGa     <  4:252778/149‑1 (MQ=255)
                                                                                                                                             aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCAGGCAGTGCTGGAAGTGGAAgacgt                                               >  8:115364/1‑106 (MQ=255)
                                                                                                                                                gAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCg  >  1:77656/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAG‑GGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCG  >  NC_002947/1622045‑1622336

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