Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,098,171 C→G 5.2% M407I (ATG→ATC)  acsA‑I ← acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,098,1710CG 5.2% 237.5 / 3.8 78M407I (ATG→ATC) acsA‑Iacetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base G (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

GTGTCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTA  >  NC_002947/5098028‑5098311
                                                                                                                                                |                                                                                                                                             
gtgtCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGgc                                                                                                                                          >  8:83068/1‑149 (MQ=255)
gtgtCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGgc                                                                                                                                          <  7:79174/149‑1 (MQ=255)
   tCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGCTCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGgcgc                                                                                                                                        <  3:269099/148‑1 (MQ=255)
     gacgaTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGaa                                                                                                                                     >  3:104675/1‑149 (MQ=255)
     gacgaTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAAACGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGaa                                                                                                                                     >  6:302855/1‑149 (MQ=255)
      acgaTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGa                                                                                                                                      >  4:164977/1‑147 (MQ=255)
        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                                                                                                               <  7:262097/148‑1 (MQ=255)
                           ccGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTaga                                                                                                               <  4:104675/148‑1 (MQ=255)
                           ccGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTaga                                                                                                               <  7:83068/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                                                                                                <  4:116326/149‑1 (MQ=255)
                                         gtaCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTtgt                                                                                                >  5:80438/1‑149 (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                                                                                 >  1:377580/1‑149 (MQ=255)
                                                        cTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCt                                                                                 >  4:143333/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                                                                               >  4:372389/1‑131 (MQ=255)
                                                                            ccaccGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTc                                                                               <  3:372389/131‑1 (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                                         <  6:12645/129‑1 (MQ=255)
                                                                                                                    ttCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTc                                         >  5:12645/1‑129 (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                <  3:375056/149‑1 (MQ=255)
                                                                                                                         cAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTg                <  1:277490/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)
                                                                                                                                                |                                                                                                                                             
GTGTCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTT‑GGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTA  >  NC_002947/5098028‑5098311

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: ‑