Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,098,182 T→G 7.7% I404L (ATT→CTT)  acsA‑I ← acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,098,1820TG 7.7% 231.0 / 8.1 77I404L (ATT→CTT) acsA‑Iacetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base G (2/4);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

GATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGG  >  NC_002947/5098036‑5098327
                                                                                                                                                   |                                                                                                                                                  
gaTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATgg                                                                                                                                                  <  5:310000/149‑1 (MQ=255)
 aTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATgg                                                                                                                                                  >  5:139393/1‑148 (MQ=255)
    gggCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAg                                                                                                                                               <  3:182214/148‑1 (MQ=255)
    gggCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCATTCGCGCCATGATGGC‑CGAAGGGCAg                                                                                                                                               <  5:316993/148‑1 (MQ=255)
     ggCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGtt                                                                                                                                             >  3:116326/1‑149 (MQ=255)
        aaCGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑Ggc                                                                                                                                         >  8:28502/1‑149 (MQ=255)
         aCGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑Ggcg                                                                                                                                        <  5:35496/149‑1 (MQ=255)
          cGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑Ggcgc                                                                                                                                       <  3:164977/149‑1 (MQ=255)
          cGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCATTCGCGCCATGATGGC‑CGAAGGGCAggcagctg                                                                                                                                        <  3:26236/149‑8 (MQ=38)
                  gCCGACGGTCTTGTAGTACGAGTTCCAGGCCTCAGGGTTGAGCGGCTCACCCACCGTCCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTag                                                                                                                                <  1:210463/148‑1 (MQ=255)
                   ccGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTaga                                                                                                                               <  4:104675/148‑1 (MQ=255)
                   ccGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTaga                                                                                                                               <  7:83068/148‑1 (MQ=255)
                   ccGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTaga                                                                                                                               <  7:262097/148‑1 (MQ=255)
                          tCTTGTAGTTGCAGGTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGAGGCAGGCTTGAACGGTCGGCGCCTTCATCAGCTGCCTGCCCTTCG‑GCCATCTTGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTg                                                                                                                       <  4:9388/149‑1 (MQ=255)
                            ttGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGAc                                                                                                                     >  6:25071/1‑149 (MQ=255)
                            ttGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTAAG‑GCCAACATGGCGCGAATGGCACTT‑GGCGCAGTGTAGAGGATGTTGAc                                                                                                                     >  8:310069/1‑149 (MQ=255)
                              gtagtaCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGAcc                                                                                                                    <  8:114520/148‑1 (MQ=255)
                                 gtaCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTtgt                                                                                                                >  5:80438/1‑149 (MQ=255)
                                 gtaCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTtgt                                                                                                                <  4:116326/149‑1 (MQ=255)
                                 gtaCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTtg                                                                                                                 <  2:380872/148‑1 (MQ=255)
                                  taCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTtgtg                                                                                                               <  3:110958/149‑1 (MQ=255)
                                  taCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTtgt                                                                                                                <  7:310069/148‑1 (MQ=255)
                                    ccAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCt                                                                                                             >  4:199552/1‑149 (MQ=255)
                                    ccAGTTCCAGGCCTCAGGGGTGATCGGCTCACCGACCGACCCCAGCAGAGGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGc                                                                                                              <  1:223534/148‑1 (MQ=255)
                                     cAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTa                                                                                                                                 <  5:380172/128‑1 (MQ=255)
                                     cAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTa                                                                                                                                 >  6:380172/1‑128 (MQ=255)
                                           cAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTcgac                                                                                                      >  4:282669/1‑149 (MQ=255)
                                            aGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTcgac                                                                                                      >  2:116361/1‑148 (MQ=255)
                                               ccTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATc                                                                                                  >  8:154613/1‑149 (MQ=255)
                                                cTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCt                                                                                                 >  4:143333/1‑149 (MQ=255)
                                                cTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCt                                                                                                 >  1:377580/1‑149 (MQ=255)
                                                   aGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCtt                                                                                                <  6:139393/147‑1 (MQ=255)
                                                      gTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATa                                                                                           <  8:201002/149‑1 (MQ=255)
                                                         gATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATAcgc                                                                                        >  4:109982/1‑149 (MQ=255)
                                                             ggCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCggg                                                                                     <  3:282669/148‑1 (MQ=255)
                                                                tCACCCACCGACCCCAGCAGACGCAGGCTGGAAACGTCGGCGCCTTAAACACCTGCCTGCCCATCG‑GCCAAAATGGCGCGAATGGGAGTT‑GGCCCTGTGTAGAGGATGTTGACCTTGTGCTTGTCCACGATCATAGCAACCCGGGTAAt                                                                                 >  1:249419/1‑149 (MQ=255)
                                                                    ccaccGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTc                                                                                               <  3:372389/131‑1 (MQ=255)
                                                                    ccaccGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTc                                                                                               >  4:372389/1‑131 (MQ=255)
                                                                    ccaccGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTcc                                                                             >  5:24158/1‑149 (MQ=255)
                                                                      accGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCgg                                                                           <  3:143333/149‑1 (MQ=255)
                                                                      accGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCgg                                                                           >  2:277490/1‑149 (MQ=255)
                                                                      accGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCgg                                                                           <  5:302855/149‑1 (MQ=255)
                                                                       cccaCCCCACCATCCGCAGGGGGGAACCGGGGGGGCCTTAAAAACGTGCCTGCCCTTCG‑GCCATAATGGGGGGAATGGCAGTTCGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCgg                                                                           <  1:116361/146‑1 (MQ=255)
                                                                          aCCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAg                                                                       >  8:308003/1‑149 (MQ=255)
                                                                          aCCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTA‑GCCGCGGTGTAGAGGAAGTTGCCCTAGTGCTTGTCCACGATCTTCGAAACTCGCGTAATGTCCGAGTAg                                                                       >  4:213865/1‑149 (MQ=255)
                                                                            cccAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGtt                                                                     >  4:375056/1‑149 (MQ=255)
                                                                              cagcagACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGc                                                                                                                                                 >  1:130067/1‑72 (MQ=255)
                                                                              cagcagACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGc                                                                                                                                                 <  2:130067/72‑1 (MQ=255)
                                                                              cagcagACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTgg                                                                   >  3:21965/1‑149 (MQ=255)
                                                                                     cGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGcc                                                            >  5:143669/1‑149 (MQ=255)
                                                                                            tGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGGGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGaa                                                      >  1:246347/1‑148 (MQ=255)
                                                                                                    tCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAggt                                             <  6:24158/149‑1 (MQ=255)
                                                                                                      ggCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCGGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAggtg                                            <  6:80438/148‑1 (MQ=255)
                                                                                                       gcgcCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAggtggt                                          >  3:379858/1‑149 (MQ=255)
                                                                                                            ttCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAAt                                                                                 >  5:390155/1‑105 (MQ=255)
                                                                                                            ttCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAAt                                                                                 <  6:390155/105‑1 (MQ=255)
                                                                                                            ttCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTc                                                         >  5:12645/1‑129 (MQ=255)
                                                                                                            ttCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTc                                                         <  6:12645/129‑1 (MQ=255)
                                                                                                            ttCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGgcgc                                     >  4:287714/1‑149 (MQ=255)
                                                                                                              cAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGcc                                    <  2:292281/148‑1 (MQ=255)
                                                                                                               aaCAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGt                                  >  6:51712/1‑149 (MQ=255)
                                                                                                                aCAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGCGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGtt                                 <  2:246347/149‑1 (MQ=255)
                                                                                                                 cAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTg                                <  1:277490/149‑1 (MQ=255)
                                                                                                                 cAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTg                                <  3:375056/149‑1 (MQ=255)
                                                                                                                  aGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTgg                               >  1:28103/1‑149 (MQ=255)
                                                                                                                         tgccATTCGCGCCATGATGGC‑CGAAGGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCgg                        >  8:349404/1‑149 (MQ=38)
                                                                                                                           cccTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGc                      >  1:388050/1‑149 (MQ=255)
                                                                                                                            ccTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGcc                     >  4:61455/1‑149 (MQ=255)
                                                                                                                            cattcgcgCCATGATGGC‑CGAAGGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGcc                     >  5:111754/8‑149 (MQ=38)
                                                                                                                             cTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCg                    <  3:287714/149‑1 (MQ=255)
                                                                                                                                cG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGt                   <  4:339029/147‑1 (MQ=255)
                                                                                                                                cG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTa                  <  5:25071/148‑1 (MQ=255)
                                                                                                                                         atGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGAtata         <  7:308003/149‑1 (MQ=255)
                                                                                                                                         atGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGACACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGAtata         <  7:154613/149‑1 (MQ=255)
                                                                                                                                          agcgccgAAGGGCAGTT‑GGCGCTGTGTAGAGCATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCTCGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAg        <  6:111754/144‑1 (MQ=255)
                                                                                                                                              gcgAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGGGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCtgt    <  2:377580/149‑1 (MQ=255)
                                                                                                                                                gattGGCAGTT‑GGCGCTGTGTAGAGTTGGTTGAGCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTTGTTGGGCACGCCTTCGTACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTgg  <  2:249419/146‑1 (MQ=255)
                                                                                                                                                gAAGGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTgg  <  7:349404/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGATATAGCTGTGG  >  NC_002947/5098036‑5098327

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