Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,787,072 A→G 12.7% intergenic (+39/‑143) 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,0720AG12.7% 156.0 / 12.4 63intergenic (+39/‑143)tsf/pyrHelongation factor Ts/uridylate kinase
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (4/4);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

GTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGG  >  NC_002947/1786924‑1787218
                                                                                                                                                    |                                                                                                                                                  
gTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCTCCCCAAAGAGGCTGCCCGCTCa                                                                                                                                                    <  6:39707/149‑1 (MQ=255)
 tGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAc                                                                                                                                                   <  1:260378/149‑1 (MQ=255)
     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                                                                                                                  <  6:43859/148‑1 (MQ=255)
                                   ggCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGt                                                                                                                  <  2:112669/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                     >  6:60211/2‑147 (MQ=39)
                                                                                                                                 acAAAGAGGCTGCGCGCGTGAGCGGGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgccg                     >  3:354548/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)
                                                                                                                                                    |                                                                                                                                                  
GTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGG  >  NC_002947/1786924‑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: