Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,280,955 A→C 9.2% intergenic (+88/+50) 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,9550AC 9.2% 149.7 / 6.9 64intergenic (+88/+50)trmA/pcaGtRNA/tmRNA (uridine(54)‑C5)‑methyltransferase TrmA/protocatechuate 3,4‑dioxygenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base C (2/4);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

TTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTGGCAA  >  NC_002947/5280808‑5281096
                                                                                                                                                   |                                                                                                                                              
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)
                    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)
                                                                   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)
                                                                                                                                                   |                                                                                                                                              
TTTGACCAGTTCCCGTATACCCACCATATGGAGAGTGGGGTTCTGCTGGTTCGGCGCTGATCGTTGTATTGCCGGGGCTGCTTTGCAGCCCATTCGCGGGCAAGCCCGCTCCGACATTGATCGGCGTCGGTGTACGCAGATCGTTGTAGGAGCGGGCTTGCCCGCGAATGGGCCGCAGAGCGGCCCCAAGCTTCACCTCAGAAGTCGAAGAAGACGGTCTCCCCTTCCCCCTGAATGCGGATATCAAAGCGGTACGCCGTCTTCCCATCCACTTCGCAACGCTTGGCAA  >  NC_002947/5280808‑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: