Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,104,299 T→A 9.8% L276Q (CTG→CAG)  PP_2714 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,104,2990TA 9.8% 164.5 / 7.8 61L276Q (CTG→CAG) PP_2714sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base A (4/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

GCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGA  >  NC_002947/3104153‑3104445
                                                                                                                                                  |                                                                                                                                                   
gCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCt                                                                                                                                                     <  5:51805/147‑1 (MQ=255)
gCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCt                                                                                                                                                     <  8:15454/147‑1 (MQ=255)
gCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCt                                                                                                                                                     <  6:310788/147‑1 (MQ=255)
gCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTgg                                                                                                                                                   >  2:173608/1‑149 (MQ=255)
 cTGTCGCAGATCGACGATCTGCTGCGGCGTGAACATAGCTTCATTGCCGATTCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTgg                                                                                                                                                   <  7:260546/148‑1 (MQ=255)
  tGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGc                                                                                                                                                 >  6:63507/1‑149 (MQ=255)
     cgcAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCttt                                                                                                                                              <  3:351417/149‑1 (MQ=255)
      gcAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTc                                                                                                                                             <  8:70697/149‑1 (MQ=255)
      gcAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTc                                                                                                                                             >  7:352935/1‑149 (MQ=255)
          aTCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCGGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTtgt                                                                                                                                         <  4:166024/149‑1 (MQ=255)
             gacgaTCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTggg                                                                                                                                      >  5:314759/1‑149 (MQ=255)
                gaTCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCgg                                                                                                                                   >  3:99695/1‑149 (MQ=255)
                 aTCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCggg                                                                                                                                  >  4:117820/1‑149 (MQ=255)
                    tgctgcGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTc                                                                                                                               >  6:198586/1‑149 (MQ=255)
                     gctgcGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTc                                                                                                                               >  1:133587/1‑148 (MQ=255)
                      ctgcGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCggg                                                                                                                                  <  8:286587/144‑1 (MQ=255)
                             gTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCgg                                                                                                                                   >  2:159712/1‑136 (MQ=255)
                             gTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCgg                                                                                                                                   <  1:159712/136‑1 (MQ=255)
                                    tCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACggg                                                                                                              <  3:117820/149‑1 (MQ=255)
                                           aTTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGc                                                                                                                        >  7:27512/1‑133 (MQ=255)
                                           aTTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGc                                                                                                                        <  8:27512/133‑1 (MQ=255)
                                            ttGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAAc                                                                                                      <  6:346654/149‑1 (MQ=255)
                                               ccGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTc                                                                                                         >  6:34654/1‑143 (MQ=255)
                                               ccGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTc                                                                                                         <  5:34654/143‑1 (MQ=255)
                                                   tgctgcCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGc                                                                                               <  7:172320/149‑1 (MQ=255)
                                                       gcCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGa                                                                                            >  7:342056/1‑148 (MQ=255)
                                                       gcCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGa                                                                                            <  8:342056/148‑1 (MQ=255)
                                                         ccACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCtt                                                                                                                                           <  4:57689/100‑1 (MQ=255)
                                                           aCGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTg                                                                                       <  8:352935/149‑1 (MQ=255)
                                                                   cGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCcgcg                                                                                <  4:66887/148‑1 (MQ=255)
                                                                       cGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAg                                                                           <  6:314759/149‑1 (MQ=255)
                                                                                 cGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTgg                                                                 <  5:63507/149‑1 (MQ=255)
                                                                                  gTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTggg                                                                <  1:173608/149‑1 (MQ=255)
                                                                                   tATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTgggg                                                               >  2:166574/1‑149 (MQ=255)
                                                                                        cgcCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGc                                                                          >  5:32337/1‑133 (MQ=255)
                                                                                        cgcCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGc                                                                          <  6:32337/133‑1 (MQ=255)
                                                                                           cTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGc                                                       >  3:320036/1‑149 (MQ=255)
                                                                                               aCGCGCAGAACGCGCTGGGCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGc                                                   <  8:210407/149‑1 (MQ=255)
                                                                                                   gcAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCAccc                                               >  3:36150/1‑149 (MQ=255)
                                                                                                     aGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCgg                                             <  1:166574/149‑1 (MQ=255)
                                                                                                        aCGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGAt                                           <  2:133587/148‑1 (MQ=255)
                                                                                                          gcgcTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGa                                        >  6:173530/1‑149 (MQ=255)
                                                                                                          gcgcTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGa                                        >  8:305678/1‑149 (MQ=255)
                                                                                                            gcTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTgg                                                                 <  8:166681/122‑1 (MQ=255)
                                                                                                            gcTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTgg                                                                 >  7:166681/1‑122 (MQ=255)
                                                                                                                gccgccAGCAATGATGCCGAGCGCCAGAAAGCCATGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGc                                  >  1:134734/1‑149 (MQ=255)
                                                                                                                    ccAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCAccc                                               >  7:113684/1‑132 (MQ=255)
                                                                                                                    ccAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCAccc                                               <  8:113684/132‑1 (MQ=255)
                                                                                                                    ccAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGATCTGGCCGAg                              <  4:99695/149‑1 (MQ=255)
                                                                                                                            gatgCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCg                                <  6:296457/139‑1 (MQ=255)
                                                                                                                            gatgCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCg                                >  5:296457/1‑139 (MQ=255)
                                                                                                                              tgCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTCGGCCAGgtcgatccgg                    <  7:214154/149‑7 (MQ=255)
                                                                                                                                  gAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAAc                >  5:33211/1‑149 (MQ=255)
                                                                                                                                     cgcCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCt             >  8:20262/1‑149 (MQ=255)
                                                                                                                                     cacaaGAAAGCCCAGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCt             >  6:379215/5‑149 (MQ=255)
                                                                                                                                      acaaGAAAGCCCAGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAgg                             >  2:83135/4‑132 (MQ=255)
                                                                                                                                      acaaGAAAGCCCAGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAgg                             <  1:83135/129‑1 (MQ=255)
                                                                                                                                       ccAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGGCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTgg           <  5:198586/149‑1 (MQ=255)
                                                                                                                                       caaGAAAGCCCAGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGgtc                           >  3:119869/3‑133 (MQ=255)
                                                                                                                                       caaGAAAGCCCAGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGgtc                           <  4:119869/131‑1 (MQ=255)
                                                                                                                                        ctgAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCAcc                     <  2:263200/136‑1 (MQ=255)
                                                                                                                                        cAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGAACACGGGTGGACAACCAGTTGCTGACATTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCAcc                     >  1:263200/1‑139 (MQ=255)
                                                                                                                                        aaGAAAGCCCAGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGc          >  3:227343/2‑149 (MQ=255)
                                                                                                                                                ccTGGGCTTTCTTGTGGGCGGGGGAGATCTGCTG‑ACACTGGTGGTCATCCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGTGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGa  <  4:36150/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTG‑ACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGA  >  NC_002947/3104153‑3104445

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