Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,745,838 C→G 18.9% intergenic (+80/+672) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,745,8380CG18.9% 454.6 / 99.9 200intergenic (+80/+672)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (92/70);  new base G (20/17);  total (112/88)
Fisher's exact test for biased strand distribution p-value = 8.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

Only 100 of 204 total aligned reads displayed.

GCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTG  >  NC_002947/4745691‑4745985
                                                                                                                                                   |                                                                                                                                                   
gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCATGccc                                                                                                                                                   >  1:574627/1‑150 (MQ=255)
gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATAGCGGCCGTGgcc                                                                                                                                                   >  1:2935587/1‑150 (MQ=255)
     ccGCGCCAGCTGTCCACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGTCCGCGaa                                                                                                                                              <  2:4697508/150‑1 (MQ=255)
     ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                              <  1:2887384/150‑1 (MQ=255)
     ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                              <  1:3770147/150‑1 (MQ=255)
     ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                              <  1:2779222/150‑1 (MQ=255)
     ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                              <  1:2500540/150‑1 (MQ=255)
      cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                                                                             <  1:1617942/150‑1 (MQ=255)
      cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                                                                             >  2:472574/1‑150 (MQ=255)
      cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:3864195/150‑1 (MQ=39)
      cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:1035489/150‑1 (MQ=39)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  2:1944711/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  2:2324372/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  2:2696130/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:945738/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:2825216/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:1175202/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  2:4776996/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  2:2852015/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  1:4694013/150‑1 (MQ=255)
       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTATTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                            <  1:984532/150‑1 (MQ=255)
       gcgcCAGCTGTACAACGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                            <  1:2216936/150‑1 (MQ=255)
         gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCATAGCGGCCGGGCCCCCGaatatg                                                                                                                                          >  1:2489199/1‑146 (MQ=255)
            aGCTGTACACCGGCGTGCAGAAGCGCTACATGGTTGCCGTGAAAGACCGCATGTACGCTGCAGTCGCCGAACGCACGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCGGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                                                                       <  2:209750/150‑1 (MQ=255)
            aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                                                                       <  2:3937806/150‑1 (MQ=255)
             gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCgg                                                                                                                                      <  2:3145229/150‑1 (MQ=255)
             gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGGAGAGGCCgg                                                                                                                                      >  2:495394/1‑150 (MQ=255)
                 tACACCGGCGAGCATAAGCGCGACATCGTTGCCCTGAAATACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGACGGCTGCCCTCGGGCAGCCTTTTTTATTTGGCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                                                                  <  2:5608982/150‑1 (MQ=255)
                  acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCGGCGCGGCCGTGCCCGCGCAGAGGCCGGCACAg                                                                                                                                 >  2:37839/1‑150 (MQ=255)
                     ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCCGCCTTTCTTATGTGTCCTGTGCCGGCCGCGTCGCGGCCGGGCCCGCGCCGAGGCCGGCACAGaca                                                                                                                              >  1:4111227/1‑150 (MQ=255)
                       ggCGAGCATACGCGCGACATGGTTGCCCTCAACGACCGCAAGTACGCGGCTGTCCCCTACCGCACGAAGGCTGCCCTGGGGCAGCCTTTTTTATTTGCCCTGTCCCGGCCTCTTGGCGGCCTGGCCCCCGAAGAGGCCGGCACAGAcacc                                                                                                                            <  2:3795113/150‑1 (MQ=255)
                       ggCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAcacc                                                                                                                            >  1:3252884/1‑150 (MQ=255)
                          gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGGCTCATAGCCGCCGTGCCCGCGAAGAGGCCGGCACCAACAAcaca                                                                                                                         >  1:661427/1‑150 (MQ=255)
                          gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                                                                         <  1:3838813/134‑1 (MQ=255)
                              aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                                                                     >  2:5513035/1‑150 (MQ=255)
                              aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                                                                     >  2:1635580/1‑150 (MQ=255)
                               gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCAGACACCACACAAcc                                                                                                                    >  2:920288/1‑150 (MQ=255)
                                     cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                                                                 >  2:5617699/1‑115 (MQ=255)
                                         aTCTTTGCCCTTCACGACCGAAAGTAAGCTGCATTAGCCGAACCCAAGACGGCGGCCATGGGCCAGCCTTTTTTTTTTGACCTGTGCCGGCGTCTTCGGGGCCTTGCCCGCTATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                                                          <  2:4228103/150‑1 (MQ=255)
                                            gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCCCAAAACCTCAGTTTTTCTcc                                                                                                       <  2:878244/150‑1 (MQ=255)
                                                  cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCCCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTaccccccc                                                                                                 >  1:4135084/1‑148 (MQ=255)
                                                   tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCt                                                                                                >  2:3232634/1‑150 (MQ=255)
                                                   tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCCGTTTATCAccgcccgcg                                                                                                >  2:4566644/1‑149 (MQ=255)
                                                          ccGCAAGTAAGCTGCCGTAGCCGAACGCAAGACGGCGGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAg                                                                                         <  2:1697611/150‑1 (MQ=255)
                                                          ccGCAAGTAAGCTGCAGTAGCCGATCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTGTTGGCGGGCTTTGCGGGGAAGAGGCCGGCACAGGGACCACTAAAGCTCAGTTTTTGTCGGCCCGCTCTCGCAg                                                                                         >  2:5348015/1‑150 (MQ=255)
                                                           cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGc                                                                                        >  1:4627298/1‑150 (MQ=255)
                                                            gCAAGTACGCTGCATTAGCCGAACGAAAGACGGCTGCCCTCGGGCAGCATTTTTTTTTTGCCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                                                                       <  1:3230502/150‑1 (MQ=255)
                                                            gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCCGCCGTTCCCGCGCAGATGCCGGCACAGACACCACAAAAAAACAGATTTTATCCGCCCCCGCACGCAGcc                                                                                       >  1:2250490/1‑150 (MQ=255)
                                                                  atgctgctTTAGCCTAACGCATGATGGCTGCCCTCGGGCAGCATTTTTTTTTTGTCCTGGGCGGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                 <  1:1293059/148‑1 (MQ=255)
                                                                  aaGCTGCTGTCGCCGACCGCAAGCAGGCGGCCCCCGGCCAGCCTTTTTTATTGGGCCGGGGCCGCCCCCGGCGCGGCCTTGCCCGCGACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                 <  2:4303956/150‑1 (MQ=255)
                                                                  aaGCTGCAGTAGCCTACCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                 <  2:1049363/150‑1 (MQ=255)
                                                                  aaGCGGCCGTAGCCGACCGCAAGACGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACGTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                                 <  2:4650967/150‑1 (MQ=255)
                                                                   aGCGGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGGCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                                                <  2:4918387/150‑1 (MQ=255)
                                                                    gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCGTGTGCCGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                               <  1:5552656/150‑1 (MQ=255)
                                                                      tGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCCCGCAGCCCGTTCAgggg                                                                             >  2:3407151/1‑149 (MQ=255)
                                                                             gCCGACCGCAAGACGGCTGCCCGCGGGCCGCCTTTTTTTTTGGTCCTGTCCGGGCCTCTACGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                                                                      <  1:2089045/150‑1 (MQ=255)
                                                                                gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACGTCGCGGCCATTCCCGCGAAGAGGCCGGCCCACACACCACAAAACCACAGTTTTTATCCGCCCCCACACGCCCCCCCTGCACGGTATTGAACGgc                                                                   >  1:1778875/1‑150 (MQ=255)
                                                                                   ataaaaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGCCGGGCCCGCGCAGAGGCCGGCCCCGACCCCACCAAAACTCAGTTTTTCTCCCCCCCCGCACGCAGCCCAGTCCGGGGAGTGCACGGCGCa                                                                >  1:1335055/7‑150 (MQ=255)
                                                                                      aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCCGACACCACAAAACCGCAGTTTTTCTCCGCCCGCTCACGCCGCCCGTGCCGGGTATTGAACGGCGCCTcc                                                             >  2:3511305/1‑150 (MQ=255)
                                                                                      aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTCCGCCCGCTCACGCAGCCCGTTCACGGTATTTAACGGCGCATcc                                                             >  1:5152610/1‑150 (MQ=255)
                                                                                       agaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                                                            >  2:5559095/1‑150 (MQ=255)
                                                                                        cccggCTGCCCTCGGGCAGCCTTTTTTATTGGGCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac                                                           <  2:1671516/147‑1 (MQ=255)
                                                                                        aaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTGTCCGCCCGCTCACGCAGCCCGTTCAGGGTATTGAACGGCGCCTccac                                                           >  2:2432397/2‑150 (MQ=255)
                                                                                         atggcggcCCTGGGGCCGACTTTTTTTTTTGTCCTGTGCGGGCATCTGCGCGGCGTTGCCCGCGAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacc                                                          <  1:2706908/148‑1 (MQ=255)
                                                                                         aagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCCCGCAGCCCATTCAGGGGATTGAACGGCGCCTccccc                                                          >  2:2973035/1‑150 (MQ=255)
                                                                                                 ccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCGCCGTTTTTCTCCGCCCGCTCACGCAGCCCGTTCAGGGTATTTAAAGGCGCATCCACCACGAACtt                                                  >  2:522993/1‑150 (MQ=255)
                                                                                                       gCAGCCTTTTTTATTGGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGc                                            <  2:4095944/150‑1 (MQ=255)
                                                                                                          gatttttttttTGGGCGTGGGCCGGCCTCTGCGGGGCCTGGCCCGCGACGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccac                                         <  1:198916/147‑1 (MQ=255)
                                                                                                             ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCGCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAAGGCGCCTCCACCCCGAACGTGTTGGccaccaa                                      >  2:4280439/1‑149 (MQ=255)
                                                                                                             ttttttAGTTGTCCTGTGCCGGCCTCATCGCGGCCATGCCCGCGAAGAGGCCGGCACACACACCACAAAACCACAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAATGAACGGCGCATCCACCACGAACTTGTTGGCCACCAt                                      >  2:3897584/1‑150 (MQ=255)
                                                                                                                  tttttGGCCGGGGCCGGCCTCTGGGGGGCCGCCCCCGCCATGCGCCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCg                                 <  1:495394/148‑1 (MQ=255)
                                                                                                                           tgtgCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                        <  2:687251/150‑1 (MQ=255)
                                                                                                                                 ggCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg                  >  2:4619400/1‑150 (MQ=255)
                                                                                                                                    ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGGCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt               >  2:726464/1‑150 (MQ=39)
                                                                                                                                    ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGGCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt               >  2:1356143/1‑150 (MQ=39)
                                                                                                                                    ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt               >  2:2217861/1‑150 (MQ=255)
                                                                                                                                    ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt               >  1:4622706/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:2509188/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:3434868/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  2:2136340/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:231215/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  2:548807/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  1:268122/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  2:810719/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  2:3356262/1‑150 (MQ=255)
                                                                                                                                     tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc              >  2:1561977/1‑150 (MQ=255)
                                                                                                                                      ctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg             >  2:1464847/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg            >  1:4353484/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg            >  1:5259513/1‑150 (MQ=255)
                                                                                                                                       ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg            >  1:137196/1‑150 (MQ=255)
                                                                                                                                         cgcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCa                              >  1:3010816/1‑130 (MQ=255)
                                                                                                                                          gcggcgTTGCCCGCGAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt         <  1:3142462/150‑1 (MQ=255)
                                                                                                                                          gcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt         >  2:5004252/1‑150 (MQ=255)
                                                                                                                                          gcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt         >  1:2031995/1‑150 (MQ=255)
                                                                                                                                             gCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCa      >  1:3468396/1‑150 (MQ=255)
                                                                                                                                               gtgccccGCCATGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAcc    <  1:4726825/146‑1 (MQ=255)
                                                                                                                                               cTTGCCCGCGAAGAGGCCGGCACAGACACCACAACACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAcc    <  2:2941503/150‑1 (MQ=255)
                                                                                                                                                ttGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCt   >  2:592431/1‑150 (MQ=255)
                                                                                                                                                gggCCCGCGACGCGGCCGGCACAGACACCACAACACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCt   <  1:4396589/148‑1 (MQ=255)
                                                                                                                                                 tGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTg  <  2:5127877/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTG  >  NC_002947/4745691‑4745985

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: