Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,234,948 G→C 6.0% A39A (GCC→GCG)  PP_3711 ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,234,9480GC 6.0% 215.7 / 4.0 67A39A (GCC→GCG) PP_3711diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (27/36);  new base C (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GCGTGCGCCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGAAGCGCATGAAA  >  NC_002947/4234801‑4235096
                                                                                                                                                   |                                                                                                                                                    
gCGTGCGCCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTggg                                                                                                                                                     >  3:183784/1‑149 (MQ=255)
      gcCCGTTTCTTCACAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCGCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCCCAGCGCGCCaa                                                                                                                                               <  4:83710/149‑2 (MQ=255)
        ccGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAggg                                                                                                                                             >  8:153074/1‑149 (MQ=255)
           tttCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTc                                                                                                                                          <  7:21829/149‑1 (MQ=255)
              cttcAACGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCa                                                                                                                                        <  1:325609/148‑1 (MQ=255)
                tcAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGtt                                                                                                                                     >  7:104684/1‑149 (MQ=255)
                  aaaGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTgg                                                                                                                                   >  5:260321/1‑149 (MQ=255)
                  aaaGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTgg                                                                                                                                   <  4:183784/149‑1 (MQ=255)
                      cGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGcgac                                                                                                                               >  2:194880/1‑149 (MQ=255)
                          cGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGAt                                                                                                                            <  5:354281/148‑1 (MQ=255)
                          cGGTGCGGCTCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATa                                                                                                                           <  7:240243/149‑1 (MQ=255)
                            gTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATAcc                                                                                                                         <  2:14328/149‑1 (MQ=255)
                               ctgcGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCCCAGCGCGCCaag                                                                                                                                              <  8:334196/125‑3 (MQ=255)
                               ctgcGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCCCAGCGCGCCaag                                                                                                                                              >  7:334196/1‑123 (MQ=255)
                                  cgcgcCTGCTCCAGGTTGAATGCCGGCAACGCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAg                                                                                                                    <  7:340701/148‑1 (MQ=255)
                                   gcgcCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAggg                                                                                                                  <  6:260321/149‑1 (MQ=255)
                                   gcgcCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAggg                                                                                                                  <  2:138990/149‑1 (MQ=255)
                                       cTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGTGTTCCAGTTGGCGACGATACCTTCAGGGGTGa                                                                                                              <  1:27469/149‑1 (MQ=255)
                                       cTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGt                                                                                                                >  4:349750/1‑147 (MQ=255)
                                       cTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGa                                                                                                              <  1:21668/149‑1 (MQ=255)
                                          cTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAg                                                                                                             <  3:28669/147‑1 (MQ=255)
                                          cTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGGGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCa                                                                                                           <  6:242162/149‑1 (MQ=255)
                                                 ttGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGAt                                                                                                    <  6:199720/149‑1 (MQ=255)
                                                         cGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAAt                                                                                            >  2:142364/1‑149 (MQ=255)
                                                         cGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAAt                                                                                            >  4:18236/1‑149 (MQ=255)
                                                             aaCCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATccac                                                                                        >  1:107636/1‑149 (MQ=255)
                                                               cccGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATccacca                                                                                      >  5:62556/1‑149 (MQ=255)
                                                               cccGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATccacca                                                                                      >  3:193979/1‑149 (MQ=255)
                                                                ccGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCcaccac                                                                                     >  4:182963/1‑149 (MQ=255)
                                                                      gcgcgcTCGGCATCGGAGTAGAACAACGAGTAGTGCTGGCCGACGATATCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCACGGGTCAGCACGTAGATGGAATAATCCAACACGctctc                                                                               >  2:185777/1‑148 (MQ=255)
                                                                       cgcgcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGTGGTGAGCAGGTAGATGGCATAATCCACCACGctct                                                                                <  8:104684/147‑1 (MQ=255)
                                                                       cgcgcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATTATCCACCACGctct                                                                                <  6:8928/147‑1 (MQ=255)
                                                                        gcgcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGa                                                                             <  1:142364/149‑1 (MQ=255)
                                                                         cgcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGAt                                                                            <  6:62556/149‑1 (MQ=255)
                                                                           cTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATACTCCACCACGCTCTGGATc                                                                           >  7:321343/1‑148 (MQ=255)
                                                                           cTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCa                                                                          <  7:91710/149‑1 (MQ=255)
                                                                                  tCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCg                                                                   <  3:182963/149‑1 (MQ=255)
                                                                                    ggTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCgg                                                                  <  3:303219/148‑1 (MQ=255)
                                                                                      tgtAGAACAACGAGTAGTGCTGGCCGACGATTTCATCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGTGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGAAGGTAGATGGCATAATCCACAACGCTCTGGATCAACAGACGGTa                                                                >  6:241888/1‑148 (MQ=255)
                                                                                       gtAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACa                                                              <  3:349750/149‑1 (MQ=255)
                                                                                        tAGAACTACGAGTAGTGCTGGCCGTCGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGTTACa                                                              <  5:47486/148‑1 (MQ=255)
                                                                                          gAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACatga                                                           <  7:75855/149‑1 (MQ=255)
                                                                                           aacaacGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACatgat                                                          >  1:269505/1‑149 (MQ=255)
                                                                                              aacGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTc                                                       >  1:359862/1‑149 (MQ=255)
                                                                                                 gAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCggg                                                    <  3:152852/149‑1 (MQ=255)
                                                                                                 gAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATAAACAGCCGGTACATGATGTCggg                                                    >  7:161681/1‑149 (MQ=255)
                                                                                                   gtagtGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGctct                                                                                <  4:346706/119‑1 (MQ=255)
                                                                                                   gtagtGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGctct                                                                                >  3:346706/1‑119 (MQ=255)
                                                                                                         cTGGCGGACGATTTCCTCGGCAGGGTACCCCTTGGCGCGCTGGGCGCGAGGGTTCCAGTTGGCGACGATACCTTCTGGGGTGAGCAGGTAGTTGGCATAATCCACCACGCTGTGGATCATCAGCCGGTACATGATGGCGGGGTTtgct                                             <  1:185777/148‑1 (MQ=255)
                                                                                                          tGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCgg                                                                  <  1:42234/126‑1 (MQ=255)
                                                                                                          tGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCgg                                                                  >  2:42234/1‑126 (MQ=255)
                                                                                                           tgCCGAGGAAATCGTCGGCCCGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTtgctgc                                           >  6:277347/2‑148 (MQ=255)
                                                                                                            gCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTtgctgc                                           <  2:107636/147‑1 (MQ=255)
                                                                                                            gCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcga                                         <  4:193979/149‑1 (MQ=255)
                                                                                                            gCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCGGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcga                                         <  4:108734/149‑1 (MQ=255)
                                                                                                              cgacgaTTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcgagc                                       >  5:9892/1‑149 (MQ=255)
                                                                                                                  gaTTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACa                                   >  3:3312/1‑149 (MQ=255)
                                                                                                                     ttCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTACAGTTGGCGACGATACCTTGAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCAGGATCAACAGCCGGTAAATGATGTCGGGGTTTGCTGGGAGCGACATg                                 >  5:258151/1‑148 (MQ=255)
                                                                                                                        ctcggctcggTACCCCTTGGCGCGTTGGGCGCCAGTGTTCCAGTTGGCGACGATTCCTTCAGGGGTGAGCAGGTTGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAAc                              <  5:241888/148‑1 (MQ=255)
                                                                                                                           ggCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGc                                                                     >  3:288856/1‑106 (MQ=255)
                                                                                                                           ggCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGc                                                                     <  4:288856/106‑1 (MQ=255)
                                                                                                                                      ccctGGCGCCCAGCGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGCTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTc               >  5:367196/4‑149 (MQ=255)
                                                                                                                                      ccTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTc               <  6:9892/149‑1 (MQ=255)
                                                                                                                                      ccTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTc               <  2:359862/149‑1 (MQ=255)
                                                                                                                                       cTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTc               >  3:220111/1‑148 (MQ=255)
                                                                                                                                           tcgtgcTGGGCGCCATGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCATGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTAGATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGTAgc          <  6:186634/145‑1 (MQ=255)
                                                                                                                                             gcgcTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGAAgcg         >  2:366299/1‑148 (MQ=255)
                                                                                                                                               gcTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGAAGCGCAt      >  1:99784/1‑149 (MQ=255)
                                                                                                                                                   ggCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGAAGCGCATGaaa  <  3:18236/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGTGCGCCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGAAGCGCATGAAA  >  NC_002947/4234801‑4235096

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