Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,787,073 C→A 12.9% intergenic (+40/‑142) 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,0730CA12.9% 159.4 / 11.0 60intergenic (+40/‑142)tsf/pyrHelongation factor Ts/uridylate kinase
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base A (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 = 1.00e+00

TGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGG  >  NC_002947/1786925‑1787218
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                                  <  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                                                                                                      <  4:297360/148‑1 (MQ=255)
                                              gagaAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTcc                                                                                                      <  8:274123/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)
                                                                                                                                                    |                                                                                                                                                 
TGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCGCAATGG  >  NC_002947/1786925‑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: