Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,443,321 C→G 6.3% A487P (GCG→CCG)  PP_4781 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,443,3210CG 6.3% 204.4 / 3.9 63A487P (GCG→CCG) PP_4781sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base G (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAAT  >  NC_002947/5443181‑5443468
                                                                                                                                            |                                                                                                                                                   
cTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAgg                                                                                                                                              <  2:17622/148‑1 (MQ=255)
   cccGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGCCTGCGGGCGCAgccgg                                                                                                                                           >  1:301640/1‑148 (MQ=255)
     cGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGcggccg                                                                                                                                        >  4:30184/1‑149 (MQ=255)
     cGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGCCTGCGGGCGCAgccggcac                                                                                                                                        <  2:301640/149‑3 (MQ=255)
           cGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGc                                                                                                                                  <  4:308197/149‑1 (MQ=255)
           cGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGc                                                                                                                                  <  6:115500/149‑1 (MQ=255)
            gAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTg                                                                                                                                   >  1:359994/1‑147 (MQ=255)
                  cccGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCg                                                                                                                            >  1:142546/1‑148 (MQ=255)
                  cccGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGCCTGCGGGCGCAgccggca                                                                                                                                         <  4:272178/135‑2 (MQ=255)
                  cccGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGCCTGCGGGCGCAgccggca                                                                                                                                         >  3:272178/1‑134 (MQ=255)
                    cGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGAc                                                                                                                         <  6:143002/149‑1 (MQ=255)
                     gTCCGGCAGGTCGGCCGGAATGACCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACa                                                                                                                        <  2:186684/149‑1 (MQ=255)
                        cGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAggg                                                                                                                     <  5:103609/149‑1 (MQ=255)
                        cGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAggg                                                                                                                     >  6:1989/1‑149 (MQ=255)
                        cGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGG‑CTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGt                                                                                                                    <  5:201820/149‑1 (MQ=255)
                         ggCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGcggccg                                                                                                                                        >  6:373206/1‑129 (MQ=255)
                         ggCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGcggccg                                                                                                                                        <  5:373206/129‑1 (MQ=255)
                           cAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTaa                                                                                                                  <  3:166974/149‑1 (MQ=255)
                                           tCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCa                                                                                                   >  8:49706/1‑148 (MQ=255)
                                            ccAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCata                                                                                                 >  5:194155/1‑149 (MQ=255)
                                                  tCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTaa                                                                                                                  >  6:62195/1‑126 (MQ=255)
                                                  tCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTaa                                                                                                                  <  5:62195/126‑1 (MQ=255)
                                                    aGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCcatcat                                                                                         >  7:165210/1‑149 (MQ=255)
                                                      gTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGc                                                                                       <  3:185162/149‑1 (MQ=255)
                                                        aCAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCt                                                                                     >  6:79194/1‑149 (MQ=255)
                                                          aGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTg                                                                                   <  2:142546/149‑1 (MQ=255)
                                                              cTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTggg                                                                                 >  2:294965/1‑147 (MQ=255)
                                                                   cccAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGAccc                                                                          >  4:212977/1‑149 (MQ=255)
                                                                     cAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                        >  4:213843/1‑149 (MQ=255)
                                                                     cAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                        >  1:286565/1‑149 (MQ=255)
                                                                      aGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCcac                                                                       <  5:1989/149‑1 (MQ=255)
                                                                         gcCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGc                                                                    >  2:37522/1‑149 (MQ=255)
                                                                          ccTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGcc                                                                   <  2:286565/149‑1 (MQ=255)
                                                                           cTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGCCCCCACAGCCg                                                                  >  8:62008/1‑149 (MQ=255)
                                                                            tGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgt                                                                 >  3:246339/1‑149 (MQ=255)
                                                                               ccGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgtgt                                                               >  1:21155/1‑148 (MQ=255)
                                                                                cGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGa                                                             >  3:276251/1‑149 (MQ=255)
                                                                                 ggCATAACCGGCGCGCAAACCGTCCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACATCCGTGTGAt                                                            <  3:286360/149‑1 (MQ=255)
                                                                                 ggCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGAt                                                            <  6:194155/149‑1 (MQ=255)
                                                                                  gcttAACGGGCGCGAAACCCTTGCACCTGGGCGCGATGGCTTAAGGTTCCGGGCTGCGCCGGATGGGGGCGGGCTGGATTTTCGGCAAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATg                                                           <  1:294965/146‑1 (MQ=255)
                                                                                  gCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATg                                                           <  7:49706/149‑1 (MQ=255)
                                                                                        ccGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACt                                                      >  3:43749/1‑148 (MQ=255)
                                                                                        ccGGCGCGCAAACCCTGCAACAGCGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTg                                                     >  8:205101/1‑149 (MQ=255)
                                                                                               gcAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCata                                                                                                 >  8:186396/1‑98 (MQ=255)
                                                                                               gcAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCata                                                                                                 <  7:186396/98‑1 (MQ=255)
                                                                                               gcAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCAtt                                              >  5:332951/1‑149 (MQ=255)
                                                                                                cAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCAttt                                             <  1:37522/149‑1 (MQ=255)
                                                                                                 aaaCCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGCGCAGACCCCACAGCCGTGTGATGGTACTGAGCCAttt                                             >  4:341530/1‑148 (MQ=255)
                                                                                                    cccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGc                                         <  3:30184/149‑1 (MQ=255)
                                                                                                    cccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGc                                         <  7:205101/149‑1 (MQ=255)
                                                                                                      cTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                       >  1:196/1‑149 (MQ=255)
                                                                                                      cTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                       <  2:196/149‑1 (MQ=255)
                                                                                                       tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGcccc                                      >  4:361206/1‑149 (MQ=255)
                                                                                                       tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGcccc                                      >  1:45168/1‑149 (MQ=255)
                                                                                                            cAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCg                                 <  3:212977/149‑1 (MQ=255)
                                                                                                                      aGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGCCCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGtgtg                       >  8:180800/1‑149 (MQ=255)
                                                                                                                       ggCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGc                      <  2:45168/149‑1 (MQ=255)
                                                                                                                        gCTGAACGTTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGc                      <  4:43749/148‑1 (MQ=255)
                                                                                                                        gCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGCCAGGGggtgaa                     >  6:256572/1‑147 (MQ=255)
                                                                                                                         cTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCAt                    <  8:165210/149‑1 (MQ=255)
                                                                                                                          tGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCAt                    >  3:275234/1‑148 (MQ=255)
                                                                                                                                     ggCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga        <  5:256572/149‑1 (MQ=255)
                                                                                                                                       gtGCGCCGGCAGGGGGCCGGCGTGGAGTTCGGAGAGGGTAAGAATGGCTTCGACCATATCCATGATGCCTTGGGCAGACCGCACAGCCGTTTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCgagcgagc      <  5:79194/148‑1 (MQ=255)
                                                                                                                                           gcCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAAt  >  3:360033/1‑149 (MQ=255)
                                                                                                                                           gcCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAAt  <  2:21155/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGGCGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAAT  >  NC_002947/5443181‑5443468

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