Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,852,397 T→A 6.4% Q312L (CAG→CTG)  PP_2505 ← GAF domain/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,852,3970TA 6.4% 188.2 / 3.1 62Q312L (CAG→CTG) PP_2505GAF domain/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (2/2);  total (32/30)
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.23e-01

CCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAAT  >  NC_002947/2852249‑2852543
                                                                                                                                                    |                                                                                                                                                  
ccGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCt                                                                                                                                                    >  8:335308/1‑149 (MQ=255)
        gctGCGGACCCGCCCGGATGCGGAGATCGTAACTTCCGCAGGGTAACGCTGGCACATAGACAACGTGGTCACGTGCACCATGTTTGGCCGGATAGAGACCGCCGTCAGCGCGCGTTATTGCATCGTCAATCGGCTCACCCTGGGTGAg                                                                                                                                             >  2:178597/3‑148 (MQ=255)
         aTGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCa                                                                                                                                           <  3:91466/149‑1 (MQ=255)
            cGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGTCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCAGGGTGAGCAga                                                                                                                                         <  2:238622/148‑2 (MQ=255)
            cGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAggg                                                                                                                                        <  4:166578/149‑1 (MQ=255)
               aCCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTgg                                                                                                                                                  <  5:265030/136‑1 (MQ=255)
               aCCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTgg                                                                                                                                                  >  6:265030/1‑136 (MQ=255)
                       ggATCCGGAGATCTTAACTTCCTCAGGATAACGCTTGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCTCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACACAATg                                                                                                                             <  4:39539/149‑1 (MQ=255)
                          tCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTc                                                                                                                          <  3:280689/149‑1 (MQ=255)
                            cGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCAt                                                                                                                        <  3:165673/149‑1 (MQ=255)
                                   ctaaaCTACCTCAGTAAAACGCTGGCACTTCTACAACTTGGTCACGTCCACCATGTTTGGCCGCATCCAGACCCCCGCCACCGCGC‑TTAATGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCaa                                                                                                                               <  5:333634/131‑1 (MQ=255)
                                      aaCTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCAt                                                                                                                  <  6:72037/145‑1 (MQ=255)
                                           ccTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAgg                                                                                                         >  3:112463/1‑149 (MQ=255)
                                                 gATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgccgcc                                                                                                   >  1:218011/1‑149 (MQ=255)
                                                      cGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGt                                                                                               >  3:8705/1‑148 (MQ=255)
                                                             acaTAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGaaa                                                                                       <  7:335308/149‑1 (MQ=255)
                                                             acaTAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGaaa                                                                                       >  6:129476/1‑149 (MQ=255)
                                                             acaTAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGAAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGaaa                                                                                       >  7:340052/1‑149 (MQ=255)
                                                                  gACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGgcg                                                                                  <  8:340052/149‑1 (MQ=255)
                                                                     aaCTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACa                                                                                           <  7:183556/137‑1 (MQ=255)
                                                                     aaCTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACa                                                                                           >  8:183556/1‑137 (MQ=255)
                                                                     aaCTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGgcgcgc                                                                               >  7:332496/1‑149 (MQ=255)
                                                                     aaCTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGgcgcgc                                                                               >  3:220221/1‑149 (MQ=255)
                                                                     aaCTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGgcgcgc                                                                               >  1:254690/1‑149 (MQ=255)
                                                                       cTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTg                                                                             >  2:205294/1‑149 (MQ=255)
                                                                            tCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGAt                                                                        >  1:143397/1‑149 (MQ=255)
                                                                            tCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGAt                                                                        >  1:147972/1‑149 (MQ=255)
                                                                               cGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATtga                                                                     <  2:11221/149‑1 (MQ=255)
                                                                                gtgcaccaTGTTTGGCCGCTTAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATtgac                                                                    <  8:37432/146‑1 (MQ=255)
                                                                                gTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATtgac                                                                    <  2:218011/149‑1 (MQ=255)
                                                                                     ccaTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATtgactgac                                                                >  3:308810/1‑148 (MQ=255)
                                                                                     ccaTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATtgactgac                                                                <  4:308810/148‑1 (MQ=255)
                                                                                        tGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGc                                                              >  1:227416/1‑147 (MQ=255)
                                                                                         gTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAg                                                           >  1:288397/1‑149 (MQ=255)
                                                                                                cGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCg                                                    >  6:138683/1‑149 (MQ=255)
                                                                                                cGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCCGCCGATTGACTGACGCAAGCGCTCCg                                                    >  1:61494/1‑149 (MQ=255)
                                                                                                       agaCCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTg                                             >  2:119332/1‑149 (MQ=255)
                                                                                                          ccgccgTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgcgc                                          <  2:254690/149‑1 (MQ=255)
                                                                                                           cgccgTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGAtt                                                                       >  3:162411/1‑119 (MQ=255)
                                                                                                           cgccgTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGAtt                                                                       <  4:162411/119‑1 (MQ=255)
                                                                                                                  aGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCtt                                  <  4:112463/149‑1 (MQ=255)
                                                                                                                  aGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCtt                                  <  8:332496/149‑1 (MQ=255)
                                                                                                                   gcgcgcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCtt                                  <  5:47289/148‑1 (MQ=255)
                                                                                                                   gcgcgcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCtt                                  >  6:47289/1‑148 (MQ=255)
                                                                                                                   gcgcgcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTg                                 <  5:129476/149‑1 (MQ=255)
                                                                                                                     gcgcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCCGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGc                               <  2:61494/149‑1 (MQ=255)
                                                                                                                         tttATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTc                                                                                                                    >  3:195868/1‑60 (MQ=255)
                                                                                                                         tttATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTc                                                                                                                    <  4:195868/60‑1 (MQ=255)
                                                                                                                               gCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTaa                      <  1:205294/148‑1 (MQ=255)
                                                                                                                               gCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAAc                     >  8:28602/1‑149 (MQ=255)
                                                                                                                                catcatCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAAcc                    >  2:188022/1‑149 (MQ=255)
                                                                                                                                catcatCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAAc                     <  4:8705/148‑1 (MQ=255)
                                                                                                                                   catcaATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCg                  >  1:177925/1‑148 (MQ=255)
                                                                                                                                   catcaATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCg                  >  1:105496/1‑148 (MQ=255)
                                                                                                                                        cccTGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGc              >  6:303375/3‑147 (MQ=255)
                                                                                                                                           tGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAACGCTCCGCAATCTGCGCCAATTCTTgg                                >  7:198264/1‑126 (MQ=255)
                                                                                                                                           tGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAACGCTCCGCAATCTGCGCCAATTCTTgg                                <  8:198264/126‑1 (MQ=255)
                                                                                                                                             cTCACCCTGGGTGAGAAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTaa                                                                                             <  2:352525/63‑1 (MQ=255)
                                                                                                                                                aCCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTAc      <  1:188022/147‑1 (MQ=255)
                                                                                                                                                aCCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTAc      <  1:119332/147‑1 (MQ=255)
                                                                                                                                                  ccTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAAt  >  8:220880/1‑149 (MQ=255)
                                                                                                                                                  ccTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAAt  >  4:206014/1‑149 (MQ=255)
                                                                                                                                                  ccTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAAt  >  3:259670/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAAT  >  NC_002947/2852249‑2852543

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