Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,280,950 G→A 6.8% intergenic (+83/+55) 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,9500GA 6.8% 214.5 / 4.0 73intergenic (+83/+55)trmA/pcaGtRNA/tmRNA (uridine(54)‑C5)‑methyltransferase TrmA/protocatechuate 3,4‑dioxygenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (31/37);  new base A (2/3);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

CATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTGGCAA  >  NC_002947/5280803‑5281096
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                                                   >  5:144395/1‑125 (MQ=255)
                         ccaccaTATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGtt                                                                                                                                                   <  6:144395/125‑1 (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)
                                     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                                                 >  3:338976/1‑139 (MQ=255)
                                                                                                            gCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGat                                                 <  4:338976/139‑1 (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)
                                                                                                                       cATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCtt                            >  6:258395/1‑149 (MQ=255)
                                                                                                                               ggCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCACAGCGGTACGCCGTCTTCCCATcc                     <  8:112913/148‑1 (MQ=255)
                                                                                                                                 ggTCGGTGTACGCAGATCGTTGTCGGCGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACt                  <  6:287994/148‑1 (MQ=255)
                                                                                                                                     caccgaCGCCGATCAATGTCGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGc              <  6:51926/144‑1 (MQ=38)
                                                                                                                                              aGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTc                           <  6:272915/127‑1 (MQ=255)
                                                                                                                                              aGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTc                           >  5:272915/1‑127 (MQ=255)
                                                                                                                                                aTCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTGGc    >  4:252433/1‑148 (MQ=255)
                                                                                                                                                 tCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTGGCaa  >  6:368953/1‑149 (MQ=255)
                                                                                                                                                  cGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGATTGCGGATATAAAAGCGTTACGCCGGCTTCCCATCCAATTCGCAACGCTTGGCaa  >  7:127434/1‑148 (MQ=39)
                                                                                                                                                   |                                                                                                                                                   
CATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTGGCAA  >  NC_002947/5280803‑5281096

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: