Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,280,920 G→T 7.3% intergenic (+53/+85) trmA → / ← pcaG tRNA/tmRNA (uridine(54)‑C5)‑methyltransferase TrmA/protocatechuate 3,4‑dioxygenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,280,9200GT 7.3% 229.4 / 2.6 79intergenic (+53/+85)trmA/pcaGtRNA/tmRNA (uridine(54)‑C5)‑methyltransferase TrmA/protocatechuate 3,4‑dioxygenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (32/40);  new base T (3/4);  total (35/44)
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

GCAGGACACCCACCGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTC  >  NC_002947/5280774‑5281068
                                                                                                                                                  |                                                                                                                                                     
gCAGGACACCCACCGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGa                                                                                                                                                      >  4:229779/1‑148 (MQ=255)
  aGGACACCCACCGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACAt                                                                                                                                                   <  4:6393/149‑1 (MQ=255)
   ggACACCCACCGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACAtt                                                                                                                                                  >  7:337342/1‑149 (MQ=255)
       acccaccGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATc                                                                                                                                              <  1:357242/149‑1 (MQ=255)
        cccagcGCATAGAGCGCTGTGCATTGTTTGACCAGTTCCGTTATACCCACCATATGGAGAGTGGTGTTCTGCTGGTTAGGCTCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCTCATTCGCGGGCAAGCCAGGTCCTACTTTTATCg                                                                                                                                             <  4:140828/149‑1 (MQ=255)
         ccaccGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAACGATctg                                                                                                                                            >  2:80888/1‑147 (MQ=255)
         ccaccGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCgg                                                                                                                                            >  8:154826/1‑149 (MQ=255)
            ccGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACAATGATCGGCGt                                                                                                                                         >  4:231557/1‑149 (MQ=255)
              gCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCg                                                                                                                                       <  8:95735/149‑1 (MQ=255)
              gCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCg                                                                                                                                       <  3:119394/149‑1 (MQ=255)
                 tCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGgtg                                                                                                                                    <  6:296387/149‑1 (MQ=255)
                   gAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGgtgt                                                                                                                                   <  1:145241/148‑1 (MQ=255)
                    aGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTAc                                                                                                                                 >  3:292641/1‑149 (MQ=255)
                     gctcTGTGCATTGTTTGACCAGTTCCCGTAAACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACg                                                                                                                                <  8:106323/146‑1 (MQ=255)
                             cATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCAGCTCCGACATTTATAGGCGTCGGTGTAAGCAGATCGt                                                                                                                        >  7:196088/1‑149 (MQ=255)
                              aTTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGTGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGtt                                                                                                                       <  3:342080/149‑1 (MQ=255)
                              aTTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGTTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGTGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGtt                                                                                                                       <  5:105426/149‑1 (MQ=255)
                              aTTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGTCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGtt                                                                                                                       >  5:205376/1‑149 (MQ=255)
                               ttGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTg                                                                                                                      >  7:112913/1‑149 (MQ=255)
                                  tttGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAg                                                                                                                   <  5:122507/149‑1 (MQ=255)
                                  tttGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAg                                                                                                                   <  3:229779/149‑1 (MQ=255)
                                   ttGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAg                                                                                                                   <  3:231557/148‑1 (MQ=255)
                                       ccAGTTGCCGTATACCCACCATATGGAGTGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCg                                                                                                              <  3:219406/149‑1 (MQ=255)
                                       ccAGTTCCCGTATACCCACCATATGGAGAGAGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGAAGCCCATTCGCGGGCAAGCCCGCACCGAAATTGATCGGCGTCGGTGTACGCAGATCGATGTAGGACCg                                                                                                              >  4:219406/1‑149 (MQ=39)
                                        cAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCgg                                                                                                             <  5:48034/149‑1 (MQ=255)
                                        cAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCgg                                                                                                             >  3:116282/1‑149 (MQ=255)
                                        cAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCgg                                                                                                             <  3:215285/149‑1 (MQ=255)
                                        cAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGATCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCgg                                                                                                             >  1:72377/1‑149 (MQ=255)
                                          gTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGc                                                                                                           <  3:288020/149‑1 (MQ=255)
                                          gTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATAGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGc                                                                                                           >  3:58105/1‑149 (MQ=255)
                                           ttCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCCGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCt                                                                                                          >  2:181785/1‑149 (MQ=255)
                                             cccGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCtt                                                                                                         <  7:176612/148‑1 (MQ=17)
                                                gTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCGGCCCATTCGCGGGCAAGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCTTGccc                                                                                                     >  3:329478/1‑149 (MQ=14)
                                                gTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGAAAGTTGTAGGAGCGGGCTTGccc                                                                                                     >  1:375173/1‑149 (MQ=255)
                                                 tataCCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCCGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCcg                                                                                                    <  1:181785/149‑1 (MQ=255)
                                                     cccACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCcgc                                                                                                   <  8:287482/146‑1 (MQ=255)
                                                     cccACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGaa                                                                                                <  4:292641/149‑1 (MQ=255)
                                                     cccACCATAGGCAGAGGGGGGTTCTGCTGGTGCGGCGCTGATCGTTGTCTTCCCGGGGCGGCTTTGCAGCCCATTCGCGGCCAAGCCCGCTCCGACATTGATCGGCGGCGGGGTACCCAGATCGTTGTAGCAGCGGGCTTGCCCGCGaa                                                                                                <  2:375173/149‑1 (MQ=255)
                                                      ccaccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCTTGCCCGCGAAt                                                                                              <  1:80888/149‑1 (MQ=14)
                                                      ccaccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGtt                                                                                                                       <  6:144395/125‑1 (MQ=255)
                                                      ccaccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGtt                                                                                                                       >  5:144395/1‑125 (MQ=255)
                                                      ccaccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGCGCGGGCTTGCCCGCGAAt                                                                                              >  6:29734/1‑149 (MQ=255)
                                                      ccaccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCACCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGAAGGAGCGGCCTTGCCCGGGact                                                                                               >  8:284167/1‑147 (MQ=255)
                                                            tatGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCg                                                                                        >  2:132415/1‑149 (MQ=255)
                                                                gagaGGGGGGTTCTGCGGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCaga                                                                                    <  4:108173/149‑1 (MQ=255)
                                                                 agagTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGgt                                                                                                                                     <  5:112473/100‑1 (MQ=255)
                                                                 agagTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGgt                                                                                                                                     >  6:112473/1‑100 (MQ=255)
                                                                  gagTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCgg                                                                                                             <  7:322988/123‑1 (MQ=255)
                                                                        ggTTCTGCTGGTTCGGCGGTTCTCGTTGTATTGCCGGGGGTGCTTTCCTGCCCTTTCGCGGGCAACCCCTCTCCTACATTGATCGGCGTCGGTGTACGAATTTCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGcccc                                                                            <  4:58105/149‑1 (MQ=37)
                                                                               ctgGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGCAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCtt                                                                      <  7:284167/148‑1 (MQ=255)
                                                                               ctgGTGCGGCGCTGCTCGTGGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTCATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGc                                                                        <  3:256886/146‑1 (MQ=255)
                                                                                            gATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGCGCGGCCCCAAGCTTCCCCTCAGAAGTCg                                                        >  7:198599/1‑149 (MQ=255)
                                                                                              tCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGCCGGTGTACGCAGATCGTTGTATGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCgaa                                                      <  7:328211/149‑1 (MQ=255)
                                                                                                    tATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGt                                                                                                                     >  7:288708/1‑81 (MQ=255)
                                                                                                    tATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGt                                                                                                                     <  8:288708/81‑1 (MQ=255)
                                                                                                     aTTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAAGTGGGCCGCAGAGCGGCCCCAAGCTTCAcc                                                                  <  1:28065/131‑1 (MQ=255)
                                                                                                     aTTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCAcc                                                                  >  2:28065/1‑130 (MQ=255)
                                                                                                           ggggCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGGGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTcc                                         <  4:116282/149‑1 (MQ=255)
                                                                                                             ggCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTcccc                                       <  3:61476/149‑1 (MQ=255)
                                                                                                                tgctTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCGTCACCTCAGAAGTCGAAGAAGACGGGCTCCCCTTc                                    >  5:271666/1‑149 (MQ=255)
                                                                                                                  ctTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGGGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTc                                    <  2:121782/147‑1 (MQ=255)
                                                                                                                      gaaGCCAATTCGCGGGCAAGCCCGATCCGAATTTGATGGGCGTCGGTGTACGCAGATCGTTGTTGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCt                               <  6:205376/146‑1 (MQ=255)
                                                                                                                      gCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCt                               >  5:116643/1‑148 (MQ=255)
                                                                                                                          cccATTCGCGGGCAAGCCCGCTCCGACATTGACCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCGTCACCGCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATg                          >  5:287994/1‑149 (MQ=255)
                                                                                                                           cctttCGGGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACg                                               <  1:3812/124‑1 (MQ=255)
                                                                                                                           ccATTCGCGGGCAAGCCCGCTCCGACATTTATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAACCTTCACCTCAGCAGTCGAAGAACACg                                               >  2:3812/1‑127 (MQ=255)
                                                                                                                            cATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAATTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCg                        >  6:163081/1‑149 (MQ=255)
                                                                                                                            cATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCg                        >  7:199564/1‑149 (MQ=255)
                                                                                                                             aTTCGCGGGCAAGTCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCgg                       <  1:132415/149‑1 (MQ=255)
                                                                                                                             aTTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCgg                       <  5:163081/149‑1 (MQ=255)
                                                                                                                              ttCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGa                      >  7:158982/1‑149 (MQ=255)
                                                                                                                                      cAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAg              <  5:29734/149‑1 (MQ=255)
                                                                                                                                       aaGCCCGCTCCTACAACGATCTGCGTACACCGACGCCGATCAATGTCGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGc             >  5:51926/1‑149 (MQ=17)
                                                                                                                                         gCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGat                     <  4:338976/139‑1 (MQ=255)
                                                                                                                                         gCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGat                     >  3:338976/1‑139 (MQ=255)
                                                                                                                                           ccGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTa         >  1:150107/1‑149 (MQ=255)
                                                                                                                                               tCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGa                             <  6:358105/125‑1 (MQ=255)
                                                                                                                                               tCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGa                             >  5:358105/1‑125 (MQ=255)
                                                                                                                                                 cGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGcc     <  6:116643/147‑1 (MQ=255)
                                                                                                                                                  gACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTc  <  2:72377/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCAGGACACCCACCGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTC  >  NC_002947/5280774‑5281068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: