Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,515,107 C→G 8.2% E307Q (GAG→CAG) ‡ cbpA ← curved DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,515,1070CG 8.2% 180.1 / 6.2 61E307Q (GAG→CAG) ‡cbpAcurved DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base G (3/2);  total (33/28)
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.55e-01

CCGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTGGCTGT  >  NC_002947/5514969‑5515253
                                                                                                                                           |                                                                                                                                                   
ccGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCtt                                                                                                                                            <  5:93508/148‑1 (MQ=255)
ccGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTg                                                                                                                                           <  2:281154/149‑1 (MQ=255)
ccGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTg                                                                                                                                           <  7:57228/149‑1 (MQ=255)
ccGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCCGAGAAAGCtg                                                                                                                                            <  8:58797/148‑2 (MQ=255)
 cGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTATTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTgg                                                                                                                                          <  1:68041/149‑1 (MQ=255)
 cGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTgg                                                                                                                                          >  2:59552/1‑149 (MQ=255)
 cGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTgg                                                                                                                                          <  1:78711/149‑1 (MQ=255)
  gTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGt                                                                                                                                         <  8:141945/149‑1 (MQ=255)
      tCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGaaa                                                                                                                                                >  8:61427/1‑138 (MQ=255)
      tCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGaaa                                                                                                                                                <  7:61427/138‑1 (MQ=255)
        gacgaTTTCGATCACTCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTTCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCaga                                                                                                                                    <  7:269714/148‑1 (MQ=255)
        gacgaTTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCagag                                                                                                                                   <  6:119818/149‑1 (MQ=255)
        gacgaTTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCagag                                                                                                                                   >  5:119818/1‑149 (MQ=255)
        gacgaTTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCagag                                                                                                                                   >  5:115707/1‑149 (MQ=255)
             tttCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTcgc                                                                                                                              <  2:137542/149‑1 (MQ=255)
             tttCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTcgc                                                                                                                              >  7:292151/1‑149 (MQ=255)
                 gATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCggg                                                                                                                           >  2:17566/1‑148 (MQ=255)
                     aCGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGa                                                                                                                       <  1:17566/148‑1 (MQ=255)
                          ttCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTc                                                                                                                 <  8:178101/149‑1 (MQ=255)
                               ccGTGACATGGGCTTCCTGACCCAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTAGTCCATTGTGTCCTCGGGTTGAACGCAGGTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGAcg                                                                                                            >  6:194377/1‑149 (MQ=255)
                                            ttCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTcgcg                                                                                                                             <  6:62840/119‑1 (MQ=255)
                                            ttCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTcgcg                                                                                                                             >  5:62840/1‑119 (MQ=255)
                                                 gACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTTCTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCt                                                                                          >  1:163972/1‑149 (MQ=255)
                                                  acacAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCtt                                                                                         >  2:202466/1‑149 (MQ=255)
                                                  acacAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCCGAGAAAGCtgcgtt                                                                                                                                        <  6:330903/102‑6 (MQ=255)
                                                  acacAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCCGAGAAAGCtgcgtt                                                                                                                                        >  5:330903/1‑97 (MQ=255)
                                                      aGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCAt                                                                                                  <  1:240765/136‑1 (MQ=255)
                                                      aGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCAt                                                                                                  >  2:240765/1‑136 (MQ=255)
                                                         gTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgc                                                                                  <  1:59552/149‑1 (MQ=255)
                                                              cATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCAt                                                                              >  5:309318/1‑148 (MQ=255)
                                                               aTGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATa                                                                             >  8:136053/1‑148 (MQ=255)
                                                                  tcaaGTTGAACGATCAGGGTGCTGCTCATGATCACTAACTCCATTGTGTCCTCGGGTTGAACGCAGCATTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAgg                                                                         >  2:273487/4‑149 (MQ=255)
                                                                  tCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAg                                                                          <  6:115707/148‑1 (MQ=255)
                                                                  tCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTAGAAACGCAGATTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAg                                                                          >  1:91688/1‑149 (MQ=255)
                                                                         gAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGAAGATTCGTCCGACGCTGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCcgc                                                                  >  3:296582/1‑149 (MQ=255)
                                                                          aaCGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCcgc                                                                  <  8:292151/148‑1 (MQ=255)
                                                                            cGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGct                                                               >  5:259468/1‑149 (MQ=255)
                                                                               tCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgcgc                                                                                <  7:200068/129‑1 (MQ=255)
                                                                               tCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgcgc                                                                                >  8:200068/1‑129 (MQ=255)
                                                                                 aGGGTGCTGCTCATGATCACTTTCTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTcc                                                           <  2:163972/148‑1 (MQ=255)
                                                                                   ggTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTg                                                                                                                                           <  8:17/66‑1 (MQ=255)
                                                                                   ggTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTg                                                                                                                                           >  7:17/1‑66 (MQ=255)
                                                                                       ctgctCATCCTCACTTACTCCATTGTGTCCGCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGTGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCtgttt                                                     <  4:219398/148‑1 (MQ=255)
                                                                                        tgctCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGt                                                                        >  4:251583/1‑128 (MQ=255)
                                                                                        tgctCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGt                                                                        <  3:251583/128‑1 (MQ=255)
                                                                                                 tCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGc                                         >  6:162494/1‑149 (MQ=255)
                                                                                                     ttACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCtt                                     >  4:123054/1‑149 (MQ=255)
                                                                                                          ccATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCg                                >  8:129058/1‑149 (MQ=255)
                                                                                                           cATTGTGTCCTCGGGTTGAACTCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCTGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAATGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCgg                               <  2:91688/149‑1 (MQ=255)
                                                                                                           cATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGccc                                 <  6:259468/147‑1 (MQ=255)
                                                                                                           cATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCgg                               <  7:136053/149‑1 (MQ=255)
                                                                                                             ttGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTa                             >  2:133452/1‑149 (MQ=255)
                                                                                                                     tCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAgg                     >  7:175132/1‑149 (MQ=255)
                                                                                                                         gTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCGGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCt                 >  8:125375/1‑149 (MQ=255)
                                                                                                                            gAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCCAGGCCCTTGCCCGGTACACGCAGGCGCTg                <  1:202466/148‑1 (MQ=255)
                                                                                                                                cGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGc           >  4:134293/1‑149 (MQ=255)
                                                                                                                                 gCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGct          >  6:290360/1‑149 (MQ=255)
                                                                                                                                 gCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGct          >  7:116016/1‑149 (MQ=255)
                                                                                                                                 gCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCAGAAAGGTTGCCGCGCTCTCCCGGTTTGGTCGCCAGGCCCTTGCCCGGTCCCCGCAGGCGATGGCCGct          >  6:10159/1‑149 (MQ=255)
                                                                                                                                 caaGCTTTCCGAGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGct          >  1:240147/3‑149 (MQ=255)
                                                                                                                                 caaGCTTTCCGAGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGct          >  5:137961/3‑149 (MQ=255)
                                                                                                                                        tctcGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTGGCt    <  3:134293/148‑1 (MQ=255)
                                                                                                                                           cGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTGGCTGt  >  6:304078/1‑147 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CCGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTGGCTGT  >  NC_002947/5514969‑5515253

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: