Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 652,764 G→C 8.1% intergenic (‑101/+17) dsbD‑I ← / ← ctpL thiol/disulfide interchange protein 1/methyl‑accepting chemotaxis protein CtpL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947652,7640GC 8.1% 193.4 / 7.6 62intergenic (‑101/+17)dsbD‑I/ctpLthiol/disulfide interchange protein 1/methyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base C (2/3);  total (29/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 = 6.83e-01

GGCCGGGCTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCT  >  NC_002947/652616‑652906
                                                                                                                                                    |                                                                                                                                               
ggCCGGGCTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCcgcg                                                                                                                                                 >  7:12924/1‑149 (MQ=255)
    gggCTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCgcggg                                                                                                                                             <  1:185902/149‑1 (MQ=255)
     ggCTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGCCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCgcgggc                                                                                                                                            >  2:63793/1‑149 (MQ=255)
        tGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGggcggc                                                                                                                                         >  2:3096/1‑149 (MQ=255)
           ccAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGATGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGCTTACCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                                      >  1:294045/1‑149 (MQ=255)
            cagcagTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGAAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                                     >  7:177808/1‑149 (MQ=255)
                agTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                 <  2:1946/149‑6 (MQ=255)
                 gTATGAACAGCATGAAACACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTCGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                <  1:249554/149‑1 (MQ=255)
                 gTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGGAGCGCCGCCCGCGCGGCGCttcgcgg                                                                                                                                <  8:279791/149‑7 (MQ=38)
                  tAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                              <  3:332311/149‑1 (MQ=255)
                       aCAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGAcc                                                                                                                         >  1:22313/1‑149 (MQ=255)
                           cATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGCTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGa                                                                                                                     <  2:294045/149‑1 (MQ=255)
                            aTGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGaa                                                                                                                    >  3:342771/1‑149 (MQ=255)
                              gAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                                        >  6:271883/1‑128 (MQ=255)
                              gAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                                        <  5:271883/128‑1 (MQ=255)
                               aaaaaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGCCCCTGAACGc                                                                                                                 >  6:188347/1‑149 (MQ=255)
                               aaaaaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGc                                                                                                                 >  6:370113/1‑149 (MQ=255)
                               aaaaaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGc                                                                                                                 <  2:22313/149‑1 (MQ=255)
                                aaaaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCt                                                                                                                >  4:227143/1‑149 (MQ=255)
                                  aaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCg                                                                                                              <  6:12526/149‑1 (MQ=255)
                                  aaCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCg                                                                                                              >  7:302909/1‑149 (MQ=255)
                                   aCAGCCGTCGCATGGCGGCCGCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCACAGGAATAATTGGCCCTAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGc                                                                                                             <  4:5732/149‑1 (MQ=25)
                                      ggcggcGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCAcc                                                                                                          <  4:313544/149‑1 (MQ=255)
                                       gcggcgCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCg                                                                                                         >  7:344969/1‑149 (MQ=255)
                                         ggcgCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGcc                                                                                                       <  1:115437/149‑1 (MQ=255)
                                            gcATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCgg                                                                                                    <  4:342771/149‑1 (MQ=255)
                                                  cGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGc                                                                                              >  7:197604/1‑149 (MQ=255)
                                                      cTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCAc                                                                                          <  8:12924/149‑1 (MQ=255)
                                                         gcgAGGTGGAACGTGCCCGGCATGATAGCGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCgcgg                                                                                                                                              >  1:17326/1‑95 (MQ=255)
                                                         gcgAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCgcgg                                                                                                                                              <  2:17326/95‑1 (MQ=255)
                                                         gcgAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCa                                                                                           <  1:63793/146‑1 (MQ=255)
                                                               tGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGc                                                                                 >  2:113036/1‑149 (MQ=255)
                                                                 gAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCt                                                                                <  3:227143/148‑1 (MQ=255)
                                                                         ccGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTg                                                                       >  6:30401/1‑149 (MQ=255)
                                                                           ggCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCCCGCGCGGCGCttc                                                                                                                                    >  3:92725/1‑85 (MQ=255)
                                                                           ggCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCCCGCGCGGCGCttc                                                                                                                                    <  4:92725/87‑3 (MQ=255)
                                                                             cATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCa                                                                   <  6:278702/149‑1 (MQ=255)
                                                                                 aTAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCAccc                                                                <  2:131680/148‑1 (MQ=255)
                                                                                   aaCGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCAcc                                                                 >  8:48254/1‑145 (MQ=255)
                                                                                   aaCGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCAcc                                                                 <  7:48254/145‑1 (MQ=255)
                                                                                   aaCGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGt                                                             >  2:144941/1‑149 (MQ=255)
                                                                                     cGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGc                                                                                 >  1:116170/1‑127 (MQ=255)
                                                                                     cGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGc                                                                                 <  2:116170/127‑1 (MQ=255)
                                                                                            gcAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGcccc                                                    >  5:211434/1‑149 (MQ=255)
                                                                                                gcataatTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGAGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCg                                                 >  5:273246/3‑148 (MQ=255)
                                                                                                gAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGc                                                <  5:151171/149‑1 (MQ=255)
                                                                                                     atTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTTTTTACGCCCCTCGCCAATg                                           <  6:54387/149‑1 (MQ=255)
                                                                                                      ttGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATgc                                          <  1:113036/149‑1 (MQ=255)
                                                                                                         gCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATgcgc                                        <  8:344969/148‑1 (MQ=255)
                                                                                                         gCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCt                                       >  5:266462/1‑149 (MQ=255)
                                                                                                             gAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAg                                   <  6:211434/149‑1 (MQ=255)
                                                                                                                  aaacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGt                              <  8:197604/149‑1 (MQ=255)
                                                                                                                        aGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCAc                        >  2:132086/1‑149 (MQ=255)
                                                                                                                        aGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCAc                        >  3:245772/1‑149 (MQ=255)
                                                                                                                         ggAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCAcc                       <  8:177808/149‑1 (MQ=255)
                                                                                                                         ggAGCGGGTGTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCAc                        <  8:302909/148‑1 (MQ=255)
                                                                                                                           aGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCAcc                       >  2:227345/1‑147 (MQ=255)
                                                                                                                            gCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAg                    >  3:325927/1‑149 (MQ=255)
                                                                                                                                    aCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGt               >  2:162572/1‑146 (MQ=255)
                                                                                                                                    aCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGt               <  1:162572/146‑1 (MQ=255)
                                                                                                                                     cccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGaa            <  1:227345/148‑1 (MQ=255)
                                                                                                                                     cccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAAt           <  6:266462/149‑1 (MQ=255)
                                                                                                                                      ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTc                             <  1:358940/130‑1 (MQ=255)
                                                                                                                                      ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTc                             >  2:358940/1‑130 (MQ=255)
                                                                                                                                       cgcgAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATcc         >  1:339058/1‑149 (MQ=255)
                                                                                                                                             gcgcCGCCCGCGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCcgctcgc   >  3:274734/1‑149 (MQ=255)
                                                                                                                                              cgccgcGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCcgctcgct  <  3:105593/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGCCGGGCTGGCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCT  >  NC_002947/652616‑652906

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: