Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,787,077 C→G 13.3% intergenic (+44/‑138) tsf → / → pyrH elongation factor Ts/uridylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,787,0770CG13.3% 164.1 / 15.1 60intergenic (+44/‑138)tsf/pyrHelongation factor Ts/uridylate kinase
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base G (4/4);  total (30/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.80e-01

GAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGG  >  NC_002947/1786929‑1787218
                                                                                                                                                    |                                                                                                                                             
gAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgcgc                                                                                                                                               <  3:343954/149‑1 (MQ=255)
 aaaTCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgcgcg                                                                                                                                              >  8:31802/1‑149 (MQ=255)
        ttCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCctc                                                                                                                                       >  6:285717/1‑149 (MQ=255)
         tCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCGCGCGTGAGCGGGCAGCctc                                                                                                                                       <  3:32435/148‑1 (MQ=255)
         tCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCctc                                                                                                                                       <  8:25155/148‑1 (MQ=255)
         gccTTCGCCTACTTTGAAGTAGGTGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCctct                                                                                                                                      >  2:223669/2‑149 (MQ=255)
          ccTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCGCGCGTGAGCGGGCAGCctct                                                                                                                                      <  4:21829/148‑1 (MQ=255)
           cTTCGCCTACTTTGAAGTAGGTGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCttt                                                                                                                                    >  5:348769/1‑149 (MQ=255)
               aCCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCGCGCGTGAGCGGGCAGCCTCTTTgggg                                                                                                                                >  3:119961/1‑146 (MQ=38)
               aCCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTc                                                                                                                                 >  5:43859/1‑148 (MQ=255)
                ccTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCa                                                                                                                                >  6:26686/1‑148 (MQ=255)
                 cTACTTCAAAGTAGGCAAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCaa                                                                                                                               >  3:21169/1‑148 (MQ=255)
                     ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACgg                                                                                                                          >  6:24695/1‑149 (MQ=255)
                         aaGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGggag                                                                                                                      <  6:221577/149‑1 (MQ=255)
                              ggCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGt                                                                                                                  <  2:112669/148‑1 (MQ=255)
                              ggCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGt                                                                                                                  <  6:43859/148‑1 (MQ=255)
                               gCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCACGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGtt                                                                                                                 >  8:73999/1‑148 (MQ=255)
                                 gAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTta                                                                                                               >  3:297360/1‑148 (MQ=255)
                                  aaGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGAGAGCCCCGTCTGTCGCCCCAATGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTtata                                                                                                             <  7:203715/149‑1 (MQ=255)
                                      cATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGAc                                                                                                         >  6:302676/1‑149 (MQ=255)
                                       aTCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACt                                                                                                        <  5:26686/149‑1 (MQ=255)
                                          gagaAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTccc                                                                                                     <  6:19536/149‑1 (MQ=255)
                                          gagaAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTcc                                                                                                      <  8:274123/148‑1 (MQ=255)
                                          gagaAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTcc                                                                                                      <  4:297360/148‑1 (MQ=255)
                                                     cgacgaCTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTgg                                                                                           >  2:34887/1‑148 (MQ=255)
                                                      gacgacTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGc                                                                                          <  1:223669/148‑1 (MQ=255)
                                                       acgacTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCa                                                                                                                                            <  8:135137/97‑1 (MQ=255)
                                                       acgacTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCa                                                                                                                                            >  7:135137/1‑97 (MQ=255)
                                                       acgacTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCAc                                                                                        <  5:285717/149‑1 (MQ=255)
                                                            ttCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAg                                                                                   >  5:133138/1‑149 (MQ=255)
                                                            tacGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGggagga                                                                                                                    >  4:346141/3‑116 (MQ=255)
                                                            tacGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGggagga                                                                                                                    <  3:346141/114‑1 (MQ=255)
                                                                   ctgAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTg                                                                              <  6:133138/147‑1 (MQ=255)
                                                                   ctgAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGc                                                                             >  1:152444/1‑148 (MQ=255)
                                                                     gAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCAc                                                                          >  6:144310/1‑149 (MQ=255)
                                                                           gCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCgcttcgg                                                                                >  5:218727/1‑131 (MQ=255)
                                                                           gCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCgcttcgg                                                                                <  6:218727/137‑7 (MQ=255)
                                                                            ccGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCTAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACa                                                                    <  5:144310/148‑1 (MQ=255)
                                                                                  aGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGc                                                              <  5:302676/148‑1 (MQ=255)
                                                                                        ctgctgCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGc                                                       <  5:24695/149‑1 (MQ=255)
                                                                                           ctgcCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCCCGCTACAGTTAGCAGGCTGCaaa                                                    >  6:317041/1‑149 (MQ=255)
                                                                                                aaGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACCCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGcc                                               >  5:158456/1‑149 (MQ=255)
                                                                                                               cccGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTa                                <  1:34887/149‑1 (MQ=255)
                                                                                                                ccGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTaa                               <  2:152444/149‑1 (MQ=255)
                                                                                                                      gttttgacAAAGAGGCTGCGCGCGTGAGCGGGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAc                         >  2:2100/8‑149 (MQ=38)
                                                                                                                         gCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgcc                      >  3:83358/1‑149 (MQ=255)
                                                                                                                           cccAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgccg                     >  7:253296/1‑148 (MQ=255)
                                                                                                                           cccAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgccg                     >  6:3675/1‑148 (MQ=255)
                                                                                                                            acAAAGAGGCTGCGCGCGTGAGCGGGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgccg                     >  3:354548/2‑147 (MQ=39)
                                                                                                                            acAAAGAGGCTGCGCGCGTGAGCGGGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgccg                     >  6:60211/2‑147 (MQ=39)
                                                                                                                             cAAAGAGGCTGCGCGCGTGAGCGGGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAg                  <  1:2100/149‑1 (MQ=39)
                                                                                                                             cAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGt                                                                               <  6:315200/88‑1 (MQ=255)
                                                                                                                             cAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGt                                                                               >  5:315200/1‑88 (MQ=255)
                                                                                                                                 gagGCTGCGCGCGTGAGCGGGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGgaga              <  5:60211/149‑1 (MQ=255)
                                                                                                                                 gagGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGgaga              <  5:3675/149‑1 (MQ=255)
                                                                                                                                   ggCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGgagaga            <  4:21169/149‑1 (MQ=255)
                                                                                                                                    gCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGAc           <  8:253296/149‑1 (MQ=255)
                                                                                                                                      tGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTc         >  2:361685/1‑149 (MQ=255)
                                                                                                                                          cGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACTCAGCAGGAGAGACTCGCaa     >  6:215179/1‑149 (MQ=255)
                                                                                                                                           gCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCATGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAAt    <  6:348769/149‑1 (MQ=255)
                                                                                                                                            cTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACCCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATg   <  6:158456/149‑1 (MQ=255)
                                                                                                                                             tCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACACGCAATgg  >  4:283762/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGG  >  NC_002947/1786929‑1787218

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: