Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,845,237 C→A 8.0% R50R (CGG→AGG)  coaD → phosphopantetheine adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,845,2370CA 8.0% 227.0 / 8.5 75R50R (CGG→AGG) coaDphosphopantetheine adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base A (4/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

ACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTCGAGTACGAGTT  >  NC_002947/5845094‑5845385
                                                                                                                                               |                                                                                                                                                    
aCCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGtt                                                                                                                                                 <  5:281605/149‑1 (MQ=255)
aCCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGtt                                                                                                                                                 <  4:138889/149‑1 (MQ=255)
aCCGAGTGTTCTACCCGGGGACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGCTCCAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCATTGTTCCCCCTGCAACAACGGGtt                                                                                                                                                 <  5:72067/149‑1 (MQ=255)
  cGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGa                                                                                                                                               <  1:29568/149‑1 (MQ=255)
  cGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAg                                                                                                                                              >  4:250527/1‑149 (MQ=255)
        ttgtACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGc                                                                                                                                         <  7:990/149‑1 (MQ=255)
        ttgtACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGc                                                                                                                                         <  5:93822/149‑1 (MQ=255)
              ccGGGTACTTTCGACCCCATTAGAAAAGGCCATGGCGATCTGGTCGAGCGCGCTGCGCGGCTATTCGACCACGTGAACATCGCGGTGGCGGCCCGCCCCAAGAAAAACCATTTGTACCCCCTGGAACAACGGGTAGAGAGCGCCCga                                                                                                                                     >  7:58154/1‑146 (MQ=255)
              ccGGGTACTTTCGACCCCATTAGAAAAGACAATGGCAATCTGGTCGAGCGCGATTCGCGGATATTAAACAAAGTGATAAAAGAGGTGGCGGCCAGCCCAAAAAAAAACCATTTTTAACAACTGGAACAACGGTTTAAAAACGCCCGt                                                                                                                                     >  7:56680/1‑147 (MQ=255)
                gggTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAgg                                                                                                                                 <  1:222348/149‑1 (MQ=255)
                gggTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAgg                                                                                                                                 >  7:155667/1‑149 (MQ=255)
                 ggTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTc                                                                                                                               <  3:250527/149‑1 (MQ=255)
                          gACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCACCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCa                                                                                                                       >  2:259629/1‑149 (MQ=255)
                             cccATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCt                                                                                                                    >  7:143983/1‑149 (MQ=255)
                               cATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGaacaac                                                                                                                                                      <  7:160994/113‑1 (MQ=255)
                               cATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGaacaac                                                                                                                                                      >  8:160994/1‑113 (MQ=255)
                               cATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGGGGAACAAAggg                                                                                                                                                   <  6:363791/116‑1 (MQ=255)
                               cATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCGTTGTTCCAGGGGGAACAAAggg                                                                                                                                                   >  5:363791/1‑116 (MQ=255)
                                aTTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGcc                                                                                                                 >  7:28797/1‑149 (MQ=255)
                                 ttACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGccc                                                                                                                <  5:287542/149‑1 (MQ=255)
                                    ccAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGa                                                                                                                                   <  5:236852/127‑1 (MQ=255)
                                    ccAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGa                                                                                                                                   >  6:236852/1‑127 (MQ=255)
                                    ccAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAAt                                                                                                             >  6:312525/1‑149 (MQ=255)
                                      aaaGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAGCAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGt                                                                                                           >  6:166488/1‑149 (MQ=255)
                                          gCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAGCAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGaa                                                                                                       <  5:166488/149‑1 (MQ=255)
                                          gCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGaa                                                                                                       <  8:143983/149‑1 (MQ=255)
                                           ccATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGGCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAg                                                                                                      <  7:178888/149‑1 (MQ=255)
                                           ccATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAg                                                                                                      >  2:287077/1‑149 (MQ=255)
                                           ccATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAg                                                                                                      >  2:282309/1‑149 (MQ=255)
                                             aTGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGt                                                                                                     >  4:232825/1‑148 (MQ=255)
                                                 cGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCg                                                                                                >  5:78168/1‑149 (MQ=255)
                                                  gATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGc                                                                                                                                         >  5:86078/1‑107 (MQ=255)
                                                  gATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGc                                                                                                                                         <  6:86078/107‑1 (MQ=255)
                                                       ggTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTtct                                                                                          >  7:253930/1‑149 (MQ=255)
                                                       ggTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTtct                                                                                          >  2:69571/1‑149 (MQ=255)
                                                        ggcGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGGGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTtc                                                                                           <  8:58154/145‑1 (MQ=255)
                                                        gTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTtc                                                                                           <  8:56680/147‑1 (MQ=255)
                                                         tCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTctct                                                                                        >  1:135452/1‑149 (MQ=255)
                                                         tCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTctct                                                                                        <  1:69571/149‑1 (MQ=255)
                                                           gAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTc                                                                                      <  4:98409/149‑1 (MQ=255)
                                                                  cTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTg                                                                               >  8:176710/1‑149 (MQ=255)
                                                                   ttCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTACCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCAACATGGACCTGCCGATTGTCGAAGTCATCGGCTTCGCTTCCCTGTTgg                                                                              >  4:270005/1‑149 (MQ=255)
                                                                         gCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCAtt                                                                        <  8:155667/149‑1 (MQ=255)
                                                                         gCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCAtt                                                                        >  5:87894/1‑149 (MQ=255)
                                                                          cTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGTAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCAttt                                                                       <  1:259629/149‑1 (MQ=255)
                                                                             ttCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCg                                                                     <  5:49035/148‑1 (MQ=255)
                                                                               cGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGAAATGCGcc                                                                   <  2:135452/148‑1 (MQ=255)
                                                                                   cACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGAAATGCGCCaacag                                                              <  1:358078/149‑4 (MQ=255)
                                                                                      gTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGaa                                                             >  2:189153/1‑147 (MQ=255)
                                                                                         atcatcGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAggg                                                        <  1:287077/149‑1 (MQ=255)
                                                                                         atcatcGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAggg                                                        <  1:282309/149‑1 (MQ=255)
                                                                                            atcGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGgcgc                                                     >  6:284736/1‑149 (MQ=255)
                                                                                                cgGTGGCGGGCAGCCCCAAGAAAAACGCTTTGTTCGCCCTTGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAAc                                                 <  4:191822/149‑1 (MQ=255)
                                                                                                 ggTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACg                                                >  8:352775/1‑149 (MQ=255)
                                                                                                  gTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGt                                               <  5:284736/149‑1 (MQ=255)
                                                                                                      cggcCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTc                                           <  7:176710/149‑1 (MQ=255)
                                                                                                       ggcCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGccc                                                                                                                                       >  2:98823/1‑56 (MQ=255)
                                                                                                          cAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGc                                       <  5:312525/149‑1 (MQ=255)
                                                                                                              cccAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTgg                                   >  2:186610/1‑149 (MQ=255)
                                                                                                               ccAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGc                                  >  5:70108/1‑149 (MQ=255)
                                                                                                                 aaGAAAAACCCATTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGGCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCt                                <  6:295560/149‑1 (MQ=255)
                                                                                                                   gctcaaCCCGTTGTTCCAGGGGGAACAAAGGGTTGAGCTCGCCCGTGAGGTCACCGAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCtgc                              >  5:306570/5‑149 (MQ=255)
                                                                                                                   gAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGGCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCtgc                              <  7:352775/149‑1 (MQ=255)
                                                                                                                      caaCCCGTTGTTCCAGGGGGAACAAAGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCtgcgtg                           >  3:95940/2‑149 (MQ=255)
                                                                                                                       aaCCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCtgcgtgc                          >  4:149709/1‑149 (MQ=255)
                                                                                                                       aaCCCGTTGTTCCAGGGGGAACAAAGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCtgcgtgc                          >  7:78340/1‑149 (MQ=255)
                                                                                                                        aCCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAACTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCCCCAACGTCTTCCTGCGTGGCCtgcgtgcg                         <  1:119543/149‑1 (MQ=255)
                                                                                                                          ccTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGGGTGGCCTGCGTGCGgt                       >  7:108860/1‑149 (MQ=255)
                                                                                                                          ccTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGgt                       <  6:87894/149‑1 (MQ=255)
                                                                                                                               gTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCaa                                                                 <  5:144273/102‑1 (MQ=255)
                                                                                                                               gTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCaa                                                                 >  6:144273/1‑102 (MQ=255)
                                                                                                                               gTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTcc                   <  4:39278/148‑1 (MQ=255)
                                                                                                                                   ccccTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACtt              <  8:28797/149‑1 (MQ=255)
                                                                                                                                    cccTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTc             >  3:56698/1‑149 (MQ=255)
                                                                                                                                    cccTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTc             <  6:78168/149‑1 (MQ=255)
                                                                                                                                    agggggAACAAAGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTc             <  4:95940/145‑1 (MQ=255)
                                                                                                                                     ccTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTCg            <  3:232825/149‑1 (MQ=255)
                                                                                                                                          aacaacGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTCGAGTAc       >  5:384121/1‑149 (MQ=255)
                                                                                                                                           acaaaGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTCGAGTAc       <  8:78340/148‑1 (MQ=255)
                                                                                                                                            caacGGTTTGAGCTCGCACGTGAGGTCACCAAGCACCTGCCCAATTTCGAAGTCATCGTCTTATCTTCCCTGTTGGCGCATTTCTCCAAGGAACAGTGCGCCAACGTCTTc                                           <  2:91113/111‑1 (MQ=255)
                                                                                                                                            caacGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCAGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACttcttc                                           >  1:91113/1‑111 (MQ=255)
                                                                                                                                               cGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTCGAGTACGAGtt  >  8:181463/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACCGAGTGTTGTACCCGGGTACTTTCGACCCCATTACCAAAGGCCATGGCGATCTGGTCGAGCGCGCTTCGCGGCTATTCGACCACGTGATCATCGCGGTGGCGGCCAGCCCCAAGAAAAACCCTTTGTTCCCCCTGGAACAACGGGTTGAGCTCGCCCGTGAGGTCACCAAGCACCTGCCCAATGTCGAAGTCATCGGCTTCTCTTCCCTGTTGGCGCATTTCGCCAAGGAACAGGGCGCCAACGTCTTCCTGCGTGGCCTGCGTGCGGTGTCCGACTTCGAGTACGAGTT  >  NC_002947/5845094‑5845385

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: