Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 652,767 G→C 9.0% intergenic (‑104/+14) 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,7670GC 9.0% 218.6 / 10.4 76intergenic (‑104/+14)dsbD‑I/ctpLthiol/disulfide interchange protein 1/methyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base C (3/4);  total (36/40)
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.08e-01

CCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCTGTG  >  NC_002947/652627‑652909
                                                                                                                                            |                                                                                                                                              
ccAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCT                                                                                                                                        >  4:88373‑M1/1‑149 (MQ=255)
 cagcagTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTC                                                                                                                                       >  4:267097‑M1/1‑149 (MQ=255)
              aGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCt                                                                                                                          <  1:300876/149‑1 (MQ=255)
                cATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGa                                                                                                                        <  1:289321/149‑1 (MQ=255)
                  tGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAAc                                                                                                                      <  2:119891/149‑1 (MQ=255)
                        aCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCa                                                                                                                                >  6:116280/1‑133 (MQ=255)
                        aCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCa                                                                                                                                <  5:116280/133‑1 (MQ=255)
                            gcggcgCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCg                                                                                                            >  2:5632/1‑149 (MQ=255)
                              ggcgCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGcc                                                                                                          >  8:225852/1‑149 (MQ=255)
                                      gCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAg                                                                                                  <  2:227648/149‑1 (MQ=255)
                                        ggACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCt                                                                                                <  7:225852/149‑1 (MQ=255)
                                         gACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTc                                                                                               <  3:113974/149‑1 (MQ=255)
                                              gcgAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGc                                                                                          >  5:103178/1‑149 (MQ=255)
                                              gcgAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCCCGCGCGGCGCTtcgc                                                                                                                                    <  6:62970‑M1/107‑5 (MQ=255)
                                              gcgAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCCCGCGCGGCGCTtc                                                                                                                                      <  1:286556‑M1/105‑3 (MQ=255)
                                              gcgAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCCCGCGCGGCGCTtc                                                                                                                                      >  2:286556‑M1/1‑103 (MQ=255)
                                                 aGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCa                                                                                       <  3:88373/149‑1 (MQ=255)
                                                     ggAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCt                                                                                   >  1:154244/1‑149 (MQ=255)
                                                              ccGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCT                                                                                                                                  >  1:223435‑M1/1‑93 (MQ=255)
                                                              ccGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCT                                                                                                                                  <  2:223435‑M1/93‑1 (MQ=255)
                                                                 gCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTc                                                                       <  6:103178/149‑1 (MQ=255)
                                                                    tGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCAcc                                                                    >  3:270909/1‑149 (MQ=255)
                                                                       tAACGGATGCGAAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCt                                                                  <  2:230641/148‑1 (MQ=255)
                                                                        aaCGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCT                                                                                                                                  >  4:306790‑M1/1‑83 (MQ=255)
                                                                        aaCGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCT                                                                                                                                  <  3:306790‑M1/83‑1 (MQ=255)
                                                                          cGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGtt                                                               >  5:69807/1‑148 (MQ=255)
                                                                          aaatcccgcaggaatAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgtacacctcccaggcgctgaaaat                                                                                  <  2:183734‑M1/123‑48 (MQ=255)
                                                                          aaatcccgcaggaatAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcggacacctcccaggcgctgaaaat                                                                                  >  1:183734‑M1/7‑82 (MQ=255)
                                                                          aaatcccgcaggaatAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaacactcccggcgaacacctcccaggcgctgaaaaacccccggcgaactcctccca                                                              >  3:191013‑M1/7‑82 (MQ=255)
                                                                           aatcccgcaggaaTAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccag                                                             <  3:291994‑M1/144‑69 (MQ=255)
                                                                                cgcAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgct                                                        <  5:145774‑M1/149‑74 (MQ=255)
                                                                                 gcAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGcccc                                                       >  6:233243/1‑149 (MQ=255)
                                                                                 gcAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGccc                                                        >  5:245717/1‑148 (MQ=255)
                                                                                  cAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaac                                                                      <  6:165569‑M1/133‑60 (MQ=255)
                                                                                  cAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaac                                                                      >  5:165569‑M1/1‑74 (MQ=255)
                                                                                  cAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCa                                                                                              <  4:188577/109‑1 (MQ=255)
                                                                                  cAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCa                                                                                              >  3:188577/1‑109 (MQ=255)
                                                                                   aGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccgccgaacacctcccacgcgctgaaaaacacccgcccaagacctcccagacgctgaa                                                     >  2:167194‑M1/1‑73 (MQ=255)
                                                                                    ggAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaa                                                    <  1:294829‑M1/149‑78 (MQ=255)
                                                                                     gAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGc                                                   >  6:244520/1‑149 (MQ=255)
                                                                                     gAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGc                                                   >  8:235561/1‑149 (MQ=255)
                                                                                               cccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGccc                                                        <  7:225369/134‑1 (MQ=255)
                                                                                               cccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGccc                                                        >  8:225369/1‑134 (MQ=255)
                                                                                                ccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaataccccggcgaacacctcccaggcgctgcaaacctcccgcca                                         >  3:120620‑M1/1‑60 (MQ=255)
                                                                                                ccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCAccc                                                                   <  6:47845/122‑1 (MQ=255)
                                                                                                ccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCAccc                                                                   >  5:47845/1‑122 (MQ=255)
                                                                                                 cGAAACAAACGTAGGAGCGTGTTTACCCGCGAAGCTCCGCGCGTGGGGCGCTCGATCTCccaggcgctgaaacgctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaa                                       <  5:293296‑M1/149‑91 (MQ=255)
                                                                                                  gAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATgcg                                            <  4:223687/143‑1 (MQ=255)
                                                                                                  gAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATgcg                                            >  3:223687/1‑143 (MQ=255)
                                                                                                     acaaacGTAGGAGCGGGTCTCCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGtt                                   <  5:237486/149‑1 (MQ=255)
                                                                                                        aacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCgg                                                                                                       >  5:61816/1‑78 (MQ=255)
                                                                                                        aacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCgg                                                                                                       <  6:61816/78‑1 (MQ=255)
                                                                                                           gTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTcag                                                                                >  6:28196/1‑98 (MQ=255)
                                                                                                            tAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATg                                              >  7:289078/1‑131 (MQ=255)
                                                                                                            tAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATg                                              <  8:289078/131‑1 (MQ=255)
                                                                                                              ggAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCAcc                          <  7:235561/149‑1 (MQ=255)
                                                                                                              ggAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCAcc                          <  5:244520/149‑1 (MQ=255)
                                                                                                                 gCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcacccgccgaacacctcacaggcgcggcaaccctaccgccaaaaacctaaccgccaacg                       >  4:84618‑M1/1‑43 (MQ=255)
                                                                                                                 gCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAg                       >  8:59636/1‑149 (MQ=255)
                                                                                                                  cGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCa                            <  8:294575/143‑1 (MQ=255)
                                                                                                                  cGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCa                            >  7:294575/1‑143 (MQ=255)
                                                                                                                   gggTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCACCAATGCGCTGCAGGTTGTCTTCACCCAGtt                     >  7:129101/1‑149 (MQ=255)
                                                                                                                      tttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGt                  <  3:267097/149‑1 (MQ=255)
                                                                                                                      tttACACGCGAAGCGCCGCGCGGTCGGGGCGCTATCTCcgagccgctgaaactctccgggcgtacatctcccagccgctgaaaatccccgggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcga                  <  2:253170‑M1/149‑112 (MQ=255)
                                                                                                                      tatacCCGCGAAGCGCCGCGCGTGCGGCGCTCTATCTCccaggcgctgaaactctcccggcgaacacctcccagccgctgaaactcccccggccaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcga                  <  5:237346‑M1/147‑112 (MQ=255)
                                                                                                                           ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATc             >  6:302141/1‑149 (MQ=255)
                                                                                                                             gcgAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattccgc           >  3:132062‑M1/1‑31 (MQ=255)
                                                                                                                             atcgaGCGCCGCCCGCGCGGCGCTCGATCTCAGACCCTGAACGCTCGCCCCGCCCGGTTCAGCTCACCGCCCAGGCTCCGCAGCTGTTCCCCCTGGTTACGCCCCTCGCCAATGCGGTGCAGGTTGTCTTCCCCCAGTTCGTGAATCcg           >  5:307367/5‑149 (MQ=255)
                                                                                                                               cgaGCGCCGCCCGCGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGc                                                          <  4:306570/98‑1 (MQ=255)
                                                                                                                               cgaGCGCCGCCCGCGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGc                                                          >  3:306570/3‑100 (MQ=255)
                                                                                                                                aaGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCcgctc        <  6:69807/149‑1 (MQ=255)
                                                                                                                                  gcgcCTAGCGGTCGGCGCTCGATCTCcgaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccaccaaagcgcgaattccgcccatt      <  1:167194‑M1/149‑124 (MQ=255)
                                                                                                                                     ccgcccgcGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgt   >  4:157289/1‑149 (MQ=255)
                                                                                                                                      cgcgcgGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattccgcccattgccc  <  7:153099‑M1/149‑128 (MQ=255)
                                                                                                                                      cgcgcgGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgtg  <  7:59636/149‑1 (MQ=255)
                                                                                                                                      cgcgcgGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgtg  <  6:245717/149‑1 (MQ=255)
                                                                                                                                      cgcccgCGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgtg  <  3:157289/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCTGTG  >  NC_002947/652627‑652909

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: