Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,280,925 T→C 6.4% intergenic (+58/+80) 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,9250TC 6.4% 233.5 / 3.6 78intergenic (+58/+80)trmA/pcaGtRNA/tmRNA (uridine(54)‑C5)‑methyltransferase TrmA/protocatechuate 3,4‑dioxygenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (31/42);  new base C (3/2);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACACCCACCGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTT  >  NC_002947/5280777‑5281070
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                               <  3:215285/149‑1 (MQ=255)
                                     cAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCgg                                                                                                               >  3:116282/1‑149 (MQ=255)
                                     cAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCgg                                                                                                               <  5:48034/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                                                                                                                       <  8:288708/81‑1 (MQ=255)
                                                                                                 tATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGt                                                                                                                       >  7:288708/1‑81 (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)
                                                                                                                                                    |                                                                                                                                                  
GGACACCCACCGCATCGAGCGCTGTGCATTGTTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTT  >  NC_002947/5280777‑5281070

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: